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ABSTRACT 


This  report  describes  technical  findings  concerned  with  single- 
crystal  semiconductors  and  devices  and  with  the  radiation  sensitivity 
of  amorphous  semiconductors .  The  potential  of  the  photomngnetoelectric 
(PME)  and  Dember  effects  in  Cold-doped  silicon  is  suggested  for  photo¬ 
detection  in  the  visible  and  near  infra-red  regions;  measured  response 
times  and  noise  equivalent  powers  are  reported.  For  Cr-doped  gallium 
arsenide  the  development  of  a  theoretical  model  and  correlation  with 
measurements  of  the  PME  and  photoconductivity  (PC)  effects  enables 
description  of  the  recombination  and  trapping  behavior  of  photo-injected 
carriers.  For  n-type  indium  antiminide  a  study  of  the  PME  and  PC  effects 
for  temperatures  between  A.2°K  and  300°K  yields  values  for  the  carrier 
lifetimes  and  the  hole  mobility  and  demonstrates  the  negligible  influence 
of  trapping  processes.  The  presence  of  a  buried  layer  in  a  bipolar 
silicon  transistor  may  cause  the  value  of  the  collector  breakdown  voltages 
to  be  considerably  less  than  existing  theory  predicts;  we  describe  an 
effect  termed  premature  punch-through  that  accounts  for  the  observed 
behavior.  A  method  is  described  for  the  incorporation  of  systematic  model 
selection  into  the  computer  simulation  of  semiconductor  circuit  behavior; 
applicability  to  the  simulation  of  circuit  behavior  in  a  radiation 
environment  is  discussed.  For  a  KPO^-V^O,.  semiconducting  glass  we 
establish  the  fast  neutron  irradiation  damage  threshold  for  atomic 
structural  dependent  electronic  processes  to  be  5  x  10^  nvt  and  interpret 
the  radiation  behavior  by  a  theoretical  model  based  on  cluster  hopping 
processes.  The  fast  neutron  irradiation  sensitivity  of  electronic 
properties  of  heterogeneous  semiconducting  glasses  containing  crystallites 
is  established  at  fluences  from  1  x  10^  to  3  x  10^  nvt;  based  on  these 
results  and  on  softening  point  corre lations ,a  prediction  of  the  radiation 
sensitivity  of  other  amorphous  semiconductor  systems  is  made.  A  computer 
program  is  described  which  manipulates  the  approximately  12,000  data  points 
generated  in  a  combined  dielectric  and  thermogravimetric  absorption 
isotherm  study. 
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SUMMARY 


this  report,  for  the  eighth  scc*i-annual  period  of  contract  support, 
describes  technical  findings  in  two  aa.n  subject  areas:  single-crystal 
semiconductors  and  devices;  and  radiation  sensitivity  of  amorphous  semi¬ 
conductors  . 

SINGLE-CRYSTAL  SEMICONDUCTORS  ANT)  DEVICES; 

lie  suggest  the  potential  of  the  photomagnetoelectric  and  Dember  effects 
in  gold-doped  silicon  for  photodetection  in  the  visible  and  near  I-R  wave¬ 
length  region.  For  a  PME  cell  a  response  time  t  of  0.2  psec  and  a  noise 

-10  r 

equivalent  power  (NEP)  of  6x10  watts  are  demonstrated;  for  a  Dember  cell, 

the  corresponding  data  are  tr<0.03  psec  and  NEP  **  2.5  x  10  ^  watts. 

A  "heoretical  and  experimental  study  is  reported  of  the  recombination  and 

trapping  processes  governing  the  behavior  of  photoinjected  carriers  in  semi- 

insulating  Cr-doped  GaAs.  The  experimental  portion  of  the  study  consists  in 

measuring  the  photomagnetoelectric  (PME)  and  photoconductive  (PC)  effects 

between  80°K  and  300°K.  To  interpret  the  data,  we  develop  a  theory  of  the  PME 

and  PC  effects  that  takes  into  account  the  variation  of  carrier  lifetimes  with 

injected  carrier  density  and  the  trapping  of  holes  in  the  chromium  levels.  The 

theory  predicts  a  power-law  relationship  between  the  PME  short  circuit  current 

IpME  an(*  t^ie  photoconductance  AG, with  different  exponents  in  the  power  law 

applying  for  different  injection  ranges.  Two  well-defined  ranges  of  injection 

were  observed:  For  high  injection,  I  is  directly  proportional  to  AG,  and 

mh  ^2 

for  intermediate  injection,  IpME  varies  in  proportion  as  AG  .  Numerical 
values  of  the  lifetimes  Xr  and  x^  are  given  as  functions  of  AG.  The  results 
are  in  good  agreement  with  the  theoretical  predictions. 

A  study  of  the  photomagnetoelectric  (PME)  and  photoconductive  (PC)  effects 
in  N-type  InAs  single  crystals  is  described  for  the  temperature  range  between 
4.2°K  and  300°K  under  low  injection  conditions.  For  T  >  77°K  the  PME  open- 
circuit  voltage  varies  in  proportion  to  the  magnetic  induction  and  for  T  <  21°K 
the  PME  open-circuit  voltage  saturates  at  large  magnetic  induction.  These 
observations  agree  with  the  predictions  of  large  Hall-angle  PME  theory.  For 
temperatures  between  4.2 °K.  and  300°K,  c'.rrier  lifetimes  are  computed  from  the 
PME  and  PC  data  between  4.2°K  and  300°K;  the  hole  mobility  is  estimated  from 
the  PME  data  at  21°K  and  4.2°K.  The  observed  linear  relationship  between  the 
PME  open-circuit  voltage  and  photoconductance  at  low  light  levels  demonstrates 
that  trapping  contributes  negligible  effect  over  the  entire  temperature  range 
from  4. 2°  K  to  300°  K. 

ii 


Ths  presence  of  a  buried  layer  in  a  bipolar  silicon  transistor  nay  cause 
the  value  of  the  collector  breakdown  voltages  BVCES  and  BVCB0  to  be  considerably 
less  than  existing  theory  predicts.  We  describe  an  effect  termed  premature 
punch-through  that  accounts  for  the  observed  behavior.  Expressions  are  derived 
for  the  premature  punch-through  voltage  BVppT  and  the  thickness  (wEpi)opT  of 
the  epicaxial  collector  that  maximizes  this  voltage  and,  hence,  the  usable 
region  of  device  operation.  The  expression  for  BVpp>p  is  shown  to  be  consistent 
with  measurements  taken  on  transistors  fabricated  in  our  laboratory. 

*e  describe  how  systematic  selection  of  transistor  models  can  be  incorporated 
into  the  computer  simulation  of  semiconductor  circuit  behavior.  The  selection 
is  so  done  that  each  model  chosen  provides  the  best  accuracy-simplicity  compromise 
for  the  circuit  environment  in  which  each  transistor  resides.  Compared  with 
other  procedures  now  available,  the  method  described  here  simplifies  computation, 
saves  computational  cost,  and  enables  simulation  of  larger  circuits.  Applicability 
to  the  simulation  of  circuit  behavior  in  a  radiation  environment  is  discussed. 

RADIATION  SENSITIVITY  OF  AMORPHOUS  SEMICONDUCTORS: 

A  fast  neutron  irradiation  damage  threshold  for  atomic  structural  dependent 
electronic  processes  in  a  KPO^V^  semiconducting  glass  is  established  as  5  x  1016 
nvt.  This  radiation  threshold  depends  neither  upon  thermal  processing  nor  upon 
microstructure  of  the  glasses.  The  radiation  behavior  is  interpreted  in  terms 
of  a  cluster  hopping  model. 

Fast  r eutron  irradiation  sensitivity  of  certain  electronic  properties  of 
heterogeneous  semiconducting  glasses  containing  crystallites  is  established  at 
fluences  from  1-3x10  ^  nvt.  Relaxation,  loss  spectra  are  appreciably  changed  by 
the  irradiation.  Presumably,  this  occurs  because  of  disordering  of  the  phase 
barriers  which  are  responsible  for  interfacial  polarization  in  the  material. 

A  prediction  of  the  radiation  sensitivity  of  other  amorphous  semiconductor 
systems  is  made  based  on  these  results  and  on  softening  point  correlations. 

A  computer  program  is  described  which  manipulates  the  approximately  12,000 
data  points  generated  in  a  combined  dielectric  and  thermogravimetric  adsorption 
isotherm  study.  Applications  of  this  method  to  characterizing  the  surface 
dependent  behavior  of  electronic  materials  are  suggested. 
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I.  Introduction 


The  original  general  objective  of  this  research  program  was  to  establish 
at  the  University'  of  Florida  a  "Center  of  Competence  in  Solid  State  Material? 
and  Devices."  From  the  efforts  expended  In  developing  this  center  of  compe¬ 
tence  have  evolved  technical  findings:  technical  findings  concerning  such 
materials  as  glass  ceramics,  semiconducting  glasses,  magnetic  films,  degenerate 
materials  and  degenerate  semiconductors;  concerning  devices  made  from  these 
materials;  concerning  measurement  techniques;  and  concerning  methods  of  fabri¬ 
cation.  The  first  of  the  findings  In  these  various  areas  are  deser  bed  In 
seven  previous  scientific  reports.*  ^  To  Scientific  Report  No.  1  the  reader 
is  referred  for  a  more  detailed  statement  of  the  research  objectives  than  given 
here  and  for  a  discussion  of  the  means  to  be  used  in  achieving  these  objectives. 

In  accordance  with  an  agreement  made  with  the  contact  monitor  In  January 
1971,  vc  have  focussed  the  objectives  of  this  research  program  to  give  greatest 
emphasis  to  the  following  problem  areas: 

1.  The  development  of  new  types  of  detectors  with  stress  given 
to  utilizing  the  Dcmbcr  and  photosuignetoelectrlc  effects  In  semi¬ 
conductors  and  to  utilizing  new  mask  geometries  compatible  with 
conventional  silicon  technology. 

2.  Full  static  and  dynamic  characterization  of  bipolar  and  field- 
effect  devices  aiming  both  toward  Improved  device  design  and  toward 
Improved  circuit  Implementation.  The  characterization  will  deni 
with  operation  over  a  wide  range  of  temperatures,  with  emphasis  on 

the  range  between  ordinary  room  temperature  and  cryogenic  temperatures. 

3.  The  Incorporation  of  the  effects  of  Irradiation  into  the  charac¬ 
terization  of  device  behavior;  Junction  field-effect  devices  will 
receive  special  attention.  Faphasis  will  be  placed  on  the  develop¬ 
ment  of  characterizations  enabling  computer  simulation  and  design 

of  circuits. 

4.  Noise  studies  proceeding  collaterally  with  each  of  the  efforts 
listed  above.  Moreover,  studies  will  continue  dealing  with  the 
noise  associated  with  Zener  and  avalanche  phenomena  in  field-effect 
and  bipolar  devices. 

5.  Radiation  studies  of  amorphous  semiconductors. 

To  a  large  degree,  the  content  of  the  present  report  reflects  this  change 
In  emphasis.  In  the  presentation  to  follow.  Section  II  describes  the  results 
of  research  concerning  single-crystal  semiconductors  and  semiconductor  devices, 
and  Section  III  reports  findings  concerned  with  radiation  effects  in  semi¬ 
conducting  glasses. 
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II.  Single  Crystal  Semiconductors  and  Devices  (A.  J.  Brodersen,  E.  R.  Chenette, 
S .  S.  Lt  and  F.  A.  Lindholm) 

A.  PHOTODETECTION  USING  PHOTOMAGNETO ELECTRIC  AND  DEMBER  EFFECTS  IN  GOLD-DOPED 
SILICON  (H.  F.  Tseng  and  S.  S.  Li) 

The  first  report  on  utilizing  the  phot omagnetoelectric  (PME)  effect  for 
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light  detection  using  Ge  cell  was  made  by  Boatright  and  Mette.  Zitter  has 
subsequently  developed  an  InSb  PME  detector  capable  of  demodulating  optical 
masers  at  wavelength  up  to  7  ym. 

In  this  letter  we  report  our  study  of  the  PME  and  Dember  effect  detectors 

made  of  gold-doped  n-  and  p-type  silicon  wafers,  operating  between  -20°C  and 

+50°C.  It  is  found  that  the  detector's  detectivity  is  limited  by  Johnson  noise 

and  the  response  speed  is  by  RC  time  constant  of  the  detector  system. 

In  our  recent  work  we  have  developed  a  general  theory  to  account  for  the 

PME  and  Dember  effects  in  gold  over-compensated  silicon  for  the  small  signal 

3 

case.  The  results  are  that  the  PME  open  circuit  voltage  is  given  by 
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and  the  Dember  voltage  is 
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where 

VPME  iS  t,1G  PME  °pen  circuit  volta8e  Per  unit  sample  length 

Vp  is  the  open  circuit  Dember  voltage,  d  is  the  sample  thickness 

B  denotes  the  magnetic  flux  density 

U  ,  T  ,  D  ,  n  ,  p  and  b  (=u  /u  )  has  the  conventional  meaning. 

pn  n  n  o  ro  n  p 

T  is  the  ratio  of  excess  electron  and  hole  density  (F  =  An/Ap). 

AG  and  Gq  are  the  photoconductance  and  dark  conductance  per  unit 
sample  width  to  length  ratio,  respectively. 

The  above  two  equations  indicate  that  the  and  are  directly 

proportional  to  the  photoconductance  AG  which  in  turn  varies  linearly  with 
the  li'.ht  intensity  for  small  signal  cas,  .  This  is  experimentally  observed  for 
all  the  cells  reported  in  this  paper. 

The  reasons  for  using  gold-doped  silicon  as  detector  material  are  two-fold: 
(1)  gold-doped  silicon  has  extremely  short  excess  carrier  lifetime  and  therefore 
fast  response  cell  could  be  built.  (2)  High  resistivity  silicon  could  be 
achieved  by  gold  doping  and  hence  increases  the  sensitivity  of  the  PME  and 
Dember  effects  in  silicon. 
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Gold-doped  n-  and  p-type  silicon  PMF,  and  Dember  cells  with  different 

gold  density  and  shallow  impurity  concentrations  were  fabricated  and  tested. 

3 

The  typical  cell's  dimension  were  5x2xJ.l  mm  .  The  experimental  data  and 

calculated  parameters  for  several  cells  are  summarized  in  Table  I  for 

comparison.  The  spectral  response  of  the  PME  and  Dember  cells  are  similar 

to  the  conventional  silicon  photodetectors  with  the  peak  response  occurring 

near  0.6  ym,  and  the  response  is  relatively  flat  ranging  from  0.45  ym  to 

0.7  ym.  The  normalized  PME  open  circuit  voltage  V  (with  V-.,,,  at  25°C  as 

rML  rMh 

reference)  as  function  of  temperature  between  -20°C  and  +50°C  for  samples 
N-l-d,  N-2-d,  and  N-4-d  is  shown  in  Fig.  1. 

The  result  shows  that  for  each  cell  there  exists  a  peak  sensitivity 
temperature,  and  the  position  of  this  peak  shifts  to  the  lower  temperature 
side  as  the  cell's  resistivity  decreases.  The  result  also  indicates  that 
the  temperature  sensitivity  of  each  cell  around  room  temperature  is  indepen¬ 
dent  of  its  resistivity.  The  change  in  PME  response  (V’  )  for  these  cells 

is  two  percent  per  °C  around  25°  ±10°C. 

The  noise  measurement  on  PME  detectors  indicates  that  the  y  noise 
dominates  for  frequencies  below  1  KHz  and  the  thermal  noise  is  the  only 
component  for  f>l  KHz.  Further  experiments  indicate  that  the  noise  level 
does  not  change  when  the  cell  is  exposed  or  removed  from  the  light.  This 
fact  shows  that  thermal  noise  rather  than  g-r  noise  will  constitute  the  lower 
limit  for  the  detection  of  weak  signal. 

Our  PME  cell's  response  speed  measurements  indicate  that  the  detector 
response  time  is  controlled  by  the  cell's  RC  time  constant  rather  than  by 
the  excess  carrier  lifetimes.  It  is  noted  from  Table  I  that  for  cells  with 
high  internal  resistance  the  responsivity  is  high,  while  cells  with  low 
internal  resistance  have  relatively  low  responsivity  and  short  RC  time  constant. 
To  fabricate  a  PME  light  detector  of  considerable  sensitivity  and  fast  response 
speed,  it  requires  to  compromise  the  doping  concentration  of  gold  and  shallow 
impurities  in  silicon. 

The  general  features  of  a  Dember  cell  are  analog  to  that  of  a  PME 
detector.  Table  I  summarized  results  of  Au-doped  n-  and  p-  type  cells. 

The  measured  sensitivity  in  both  cells  is  two  orders  of  magnitude  lower  than 
the  value  predicted  by  Eq.  (2).  This  might  be  due  to  the  fact  that  the  photo- 
injected  carriers  have  a  very  short  diffusion  length  (of  the  order  of  few  ym) 
and  the  large  ohmic  contact  area  which  shortens  a  large  part  of  the  Dember 
voltage. 
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In  sunnary,  as  a  photodetector,  a  Derober  cell  has  a  compact  structure 
with  a  faster  response  speed  and  lower  sensitivity  than  a  PME  cell.  However, 
a  PME  cell,  aside  from  using  it  as  a  photodetector,  can  also  be  used  as  a 
Gaussmeter  by  incorporating  it  with  a  small  GaAs  light  emitting  diode. 
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TABLE  I  Performance  Characteristics  of  Au-doped  Si  PME/Dember  Detectors  at  300° K 


B.  INVESTIGATION  OF  THE  RECOMBINATION  AND  TRAPPING  PROCESSES  OF  PHOTO- INJECTED 
CARRIERS  IN  SEMI- INSULATING  Cr-DOPED  GaAs  USING  PME  AND  PC  METHODS  (S.  S.  Li 

anc  C.  I.  Huang) 

I.  Introduction 

Photoconductivity  (PC)  and  photomagnetoelectric  (PME)  effects  in  undoped 

semi-insulating  and  semiconducting  GaAs  single  crystals  were  first  reported 

1  2 

by  Holeman  and  Hilsum.  Hurd  reported  the  experimental  results  of  the  Deuber 

and  PME  effects  in  02~doped  GaAs  between  140°K  and  300°K.  The  photoe lectronic 

3 

analysis  of  high-resif tivity  n-type  GaAs  was  given  by  Bube,  based  on  the 

4 

recombination  theory  of  Rose.  A  more  recent  study  on  the  photoconductivity 
of  high  resistivity  GaAs  and  Cr-doped  GaAs  was  reported  separately  by  Broom, ^ 
and  Milner-Brown  and  Fortin. ^  These  authors  analyzed  the  observed  photocon¬ 
ductivity  data  mainly  for  the  small  injection  conditions.  In  this  paper  we 
report  the  results  of  our  measurement  of  the  PME  and  PC  effects  in  semi- 

Q 

insulating  chromium-doped  GaAs  (p  ~  10  fi-cm)  at  T  =  80°K  and  300°K.  A  gen¬ 
eralized  theory  to  account  for  the  observed  PME  and  PC  effects  in  Cr-doped 
GaAs  is  developed  by  considering  the  carrier  lifetimes  as  a  function  of 
injection  and  the  effect  of  the  trapping.  Based  on  the  theory,  we  are  able 
to  define  different  ranges  of  injections  in  which  a  power  law  relationship 
is  obtained  between  the  PME  short  circuit  current,  Ipj^E  and  the  photoconduc¬ 
tance,  ffi.  This  is  verified  by  our  PME  and  PC  experiment,  in  which  we 
observed  two  distinct  and  well-defined  injection  ranges. 

I I .  Trapping  and  Recombination  Processes 

In  steady  state  injection,  the  emission  and  capture  of  carriers  by 
the  impurity  centers  cause  the  charge  in  such  centers  to  change  from  its 
equilibrium  values.  Therefore,  the  density  of  carriers  trapped  in  the  impu¬ 
rity  states  is  a  function  of  injected  carrier  densities.  If  the  density  of 
these  traps  is  larger  than  the  excess  carrier  density,  the  charge  in  them  will 
play  an  important  role  in  preserving  the  charge  neutrality  under  steady  state 
injection.  This  has  been  shown  by  Agraz  and  Li^  in  gold-doped  silicon,  where 
they  observed  that  the  charged  impurity  centers  completely  dominate  in  charge 
neutrality  except  at  high  injection.  This  theory  has  been  verified  by  them 
from  the  results  of  PME  and  PC  measurements  in  gold-doped  silicon  under  small 

8_9 

and  intermediate  injection  conditions. 

The  recombination  model  for  the  semi-insulating  Cr-doped  GaAs  resembles 
that  of  gold-doped  silicon  due  to  the  fact  that  chromium  atoms  introduce  one 
single  deep  acceptor  level  in  GaAs.^  This  level  is  located  near  the  middle 
of  the  band  gap,  and  is  an  effective  recombination  center  for  the  excess 
carriers  under  illumination  conditions.  A  very  high  resistivity  n-type  GaAs 
can  be  obtained  by  doping  semiconducting  n-type  GaAs  with  chromium.  This  is 
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due  to  the  auto-compensation  effect  of  chromium  impurity  with  the  shallow  donor 
impurities  in  GaAs.^ 

In  order  to  understand  the  charge  states  in  the  chromium  levels  under  dark 
and  illumination  conditions,  an  energy  band  diagram  for  the  Cr-doped  GaAs  is 
presented  in  Fig.  1.  In  thermal  equilibrium,  the  chromium  levels  are  occupied 
by  the  conduction  electrons  (see  Fig.  1(a)),  and  are  in  negatively-charged 
states.  In  steady  state  injection  (e.g.,  by  means  of  photoexcitation),  the 
free  holes  in  the  valence  band  are  captured  by  the  negatively— charged  chromium 
centers  due  to  coulomb  attractive  potential.  As  a  result,  these  trapped  holes 
will  tend  to  destroy  the  charge  balance  in  the  localized  and  band  states. 

This  is  true  for  the  case  when  the  chromium  density  is  much  higher  than  the 
thermal  equilibrium  electron  concentration.  If  the  injection  is  increased 
subsequently,  more  and  more  negatively-charged  chromium  centers  will  become 
neutral  as  a  result  of  hole  trapping  at  that  level.  This  is  illustrated  in 
Fig.  1(b),  which  corresponds  to  the  situation  of  the  intermediate  injection 
case.  In  this  range  of  injection,  the  excess  carrier  densities  An  and  Ap  can 
be  correlated  by  a  trapping  constant  T  (where  F  =  ^  **  ^  < 1) .  In  addition  to 
the  trapping  effect,  the  carrier  lifetimes  may  also  depend  on  the  injected 
carrier  density  (or  incident  photon  flux  density).  This  dependence  can,  in 

general,  be  expressed  in  terms  of  a  power  law  relationship  between  T  and  An, 

S  n 

namely  =  <An  (values  of  k  and  0  will  be  discussed  later).  In  the  high 

injection  region,  the  injected  carrier  density  exceeds  the  chromium  concen¬ 
tration  (i.e..  An  =  Ap  >  N  )  and  the  chromium  acceptor  levels  are  occupied  by 
holes,  being  in  neutral  charge  states.  As  a  result,  the  charge  states  in 
localized  and  band  states  are  essentially  controlled  by  injected  carrier 
densities  (i.e..  An  =  Ap),  and  the  electron  and  hole  lifetimes  are  equal  and 
independent  of  injection.  The  trapping  constant  T,  for  this  case,  is  equal 
to  unity  and  0=0.  The  corresponding  energy  band  diagram  for  this  range  of 
injection  is  shown  in  Fig.  1(c). 

The  above  simple  physical  model  can  be  used  in  the  development  of  the 
generalized  theory  for  PME  and  PC  effects  in  chromium-doped  GaAs  under  large 
injection  conditions. 

III.  Photoconductive  and  Photomagnetoelectric  Effects 

Since  the  experimental  conditions  for  the  PME  and  PC  measurements  on  Cr- 

doped  GaAs  are  such  that  the  excess  electron  density,  An,  is  always  greater 

than  n^,  the  derivation  of  PME  and  PC  theory  will  be  made  under  the  assumption 

of  large  injection  conditions  (i.e.,  An  >>  n  and  Ap  »  p  ). 

o  o 
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The  effect  of  trapping  can  be  expressed  by 

Ap  =  TAn 

where  T  ■  Tp^Tn  <  ^  f°r  hole  trapping. 

And  the  dependence  of  carrier  lifetimes  on  the  injected  carrier  density  can 
be  written  as 

T  =  icAn^  (2) 

n 

With  the  help  of  Eqs.  (1)  and  (2),  the  generalized  expressions  for  the 
PME  short  circuit  current  and  the  photoconductance  can  be  obtained  as  follows: 
The  photoconductance  per  unit  sample  length-to-width  ratio  is  given  by:^ 


where 


and 
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is  the  effective  diffusion  constant  (3-a) 


■  20  r  for  b  »  r 
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R  *«  —  =  •  =  —  An^  ^  is  the  rate  of  recombination 

T  ,  l*  K  /*jt\ 

n  K  Lin  (3-b) 


Since  the  electron-hole  mobility  ratio,  b,  in  GaAs  is  much  greater  than 
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unity  (we  use  b  =  20  in  this  paper)  and  the  trapping  constant  r  is  smaller 
than  unity  (we  found  that  in  the  intermediate  injection  range,  T  ~  0.18  at 
300  K  and  ^  -  0.084  at  80°K),  it  is  therefore  adequate  to  assume  that  b  >>  F 
in  the  derivations  of  AG  and 

PME 

Solving  Eqs.  (3),  (3-a)  and  (3-b),  we  obtain  the  general  expression  of 
photoconductance,  AG 

A  2qUn  I-  2+S 

AG  =  (Tmr)  [ (2-p;l)  r<]'2  An  2  (4) 

P  P  O 

Eq.  (4)  is  valid  for  different  injection  ranges  so  long  as  the  conditions 
given  by  Eqs.  (1)  and  (2)  are  satisfied.  The  values  of  6  in  Eq.  (4)  may  change 
from  one  injection  range  to  the  other  and  can  be  determined  from  measurement 
of  the  slope  of  the  I  versus  AG  plot.  An  in  Eq.  (4)  denotes  the  excess 

rMfc  O 

electron  densities  at  the  illuminated  surface,  which  is  related  to  the  incident 
photon  flux  density. 
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The  PME  short  circuit  current  per  unit  width  of  the  sample,  I_,w„,  is 
y  PME 

given  by 

An 


JpHE  =  q  u  (1+b)  B  /  D  dAn 
p  o 


(5) 


Solving  Eqs.  (3-a)  and  (5),  we  obtain  the  general  expression  for  the 
PME  short  circuit  current: 
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PME 
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Tlie  relationship  between  I 


and  (6)  by  eliminating  AnQ  from  these  two  equations 


and  AG  can  now  be  obtained  from  Eqs.  (4) 

The  result  yields 
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The  general  expression  of  IpM£  versus  AG  given  in  Eq.  (7)  is  rather 
important,  since  it  provides  a  direct  correlation  between  the  two  measurable 
parameters  Ip^  and  Ag,  and  allows  us  to  make  a  direct  comparison  of  the 
theoretical  prediction  with  the  experimental  results. 

The  general  expressions  for  and  AG,  given  above,  can  be  reduced  to 

a  very  simple  form  under  certain  injection  conditions  such  as  in  the  high 
injection  region  in  which  these  parameters  are  linearly  related  to  one  another. 
To  be  more  specific,  we  have  observed  two  well-defined  injection  regions  in 
the  PME  and  PC  measurements  for  the  Cr-doped  GaAs .  This  will  be  discussed  in 
detail  as  follows: 

(1)  High  injection  region  (An  =  Ap  >>  n^ ,  N  ) 

In  this  injection  region  (see  Fig.  1-c),  the  excess  electron  and  hole 
densities  are  much  higher  than  the  thermal  equilibrium  carrier  concentration 
n^  and  the  electrically  active  chromium  density  N  .  The  charge  neutrality 
condition  is  essentially  controlled  by  the  excess  electron  and  hole  densities. 
The  carrier  lifetime  is  independent  of  light  intensity  (or  injection)^  and 
the  effect  of  trapping  is  completely  negligible.  As  a  result,  the  conditions 
given  by  Eqs.  (1)  and  (2)  can  be  reduced  to  the  forms  that 


and 


An  =  Ap  (i.e. ,  T  =  1)  (1-a) 

Tn  =  TP  =  T»  (i-e”  6  =  0)  <2-a> 
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where  Tm  Is  the  carrier  lifetime  at  high  injection  as  defined  by  Shockley  and 
Read.11 

With  the  help  of  Eqs.  (1-a)  and  (2-a),  the  photoconductance  for  this  case 
is  obtained  from  Eq .  (A)  by  setting  8  -  0  and  T  -  1,  which  yields 

fr'  c  qw  (lib)  (!)  ;  )lj  An 
I'  n  n  o 


where 


!)a  °  Hb  ls  t*lc  polar  diffusion  constant 


The  l’ME  short  circuit  current  is  obtained  from  |> .  (6)  for  J'  -  1 

Vir  “  q’.i  (1+b)  D  1!  An 
i.n.  p  a  o 

The  relationship  between  1^  and  AG  for  this  case  is  obtained  from 
Eqs.  (8)  and  (9).  The  result  yields 

D  h 

I!’MK  "  B  tT’ 

n 

Furthermore,  the  PME  open  circuit  voltage  per  unit  length  is  obtained 
from  Eq.  (10)  using  the  relation:1^ 

^ I’M F 

V  *  e  r_i\  n 

PM  I'  AG  S  1 
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The  above  results  are  consistent  with  the  PME  theory  derived  by  Van 
Roosbroeck  for  the  high  injection  case.  In  this  range  of  injection,  the 

IPME  varies  linearly  with  AG,  and  the  VpM£  is  independent  of  light  intensity 
(or  injection). 

(2)  Intermediate  injection  region  (An  >  n^.  An  <  N  ,  An  4  Ap) 

In  this  range  of  injection,  the  charge  states  in  the  chromium  levels 
contribute  a  significant  part  in  the  charge  neutrality  condition.  The  free 

holes  are  trapped  in  the  chromium  levels.  As  a  result,  the  relationship  between 
An  and  Ap  is  given  by  Eq.  (1). 

In  addition,  the  carrier  lifetime  in  this  range  of  injection  may  also 
depend  on  the  injection  (cr  light  intensity).  The  value  of  8  in  Eq.  (2), 
depending  on  the  range  of  injection,  can  only  be  determined  by  the  experimental 
observation  such  as  from  the  IpM£  versus  AG  plot  (see  Fig.  3).  In  the  Cr-doped 
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In  iddltlon,  the  carrier  llfctlae  in  th.'s  range  of  Injection  may  also 
depend  on  the  injection  (or  light  intensity).  The  value  of  8  in  Eq.  (2), 
depending  on  the  range  of  injection,  can  only  be  detominod  by  the  experi- 
ocntal  observation  such  as  from  the  I  versus  tC  plot  (see  Fig.  3).  In  the 
Cr-doped  CaAs  samples,  we  found  that  in  the  intermediate  injection  range, 
the  value  of  6  is  equal  to  -  and  therefore  the  electron  lifetime,  tn»  for 
this  case  is  given  by: 

T  -  .:An“1/3  (2-b) 

n 


With  the  help  of  Eqs.  (1)  and  (2-b),  the  photoconductancc  for  this  case 
can  be  deduced  from  Eq.  (4)  by  setting  8  "  “  ■j*  which  yields 


AC  4  <£,,,„)  |i»p  I  k-J*1  A"05/6 


(12) 


The  expression  of  I_.c  for  this  case  remains  the  same  as  that  of  Eq.  (6). 

I  rl  t  ^ 

The  IpME  in  terms  of  AC  is  obtained  by  substituting  8  ■  ~  y  Into  Eq.  (7). 

The  result  yields: 

27D2r2  ±  * 

W  “  2*" n(Hb)  n  l—1 H  ^  <13> 

ME  p  343  ^3  6q».n 

1  2 

Eq.  (13)  predicts  that  IpME  varies  as  AG  ,  if  other  parameters  on  the 
right-hand  side  are  independent  of  injection.  This  is  observed  for  sample  S-l 
at  80*K  and  sample  S-2  at  80*K  and  300°iC  In  the  intermediate  injection  range. 

IV.  Experimental  Details 

TWo  slices  of  rectangular  bar  were  cut  from  a  Cr-doped  CaAs  single  crystal 
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wafer  (total  chromium  density,  1-3  x  10  cm  )  made  by  zone-melting  method. 

3 

The  sample  dimensions  are  0.5  x  1  x  0.1  cm  for  sample  #  (S-l)  and  0.4  x  0.8  x 
o 

0.1  cm  for  sample  t  (S-2).  The  samples  were  mechanically  lapped  on  all  six 
faces  and  chemically  etched  on  the  illuminated  surface  by  using  a  solution 
mixed  with  3H2SO^:H202:H20.  Ohmic  contact  was  made  by  indium  alloying. 
Resistivity,  Hall  effect  and  Photo-Hall  effect  measurements  were  made  on  both 
samples  at  80#K  and  300#K. 

8  8 

The  dark  resistivity  for  sample  I  S-l  is  1.2  x  10  ohm-cm  and  3  x  10  ohra-eo 

2 

f*>r  sample  I  S-2  at  300*K.  The  Hall  mobility  was  found  to  be  about  250  cm  /V-soc 
at  30<f  K.  The  low  Hall  mobility  is  attributed  to  the  mixed  conduction  in  Cr-doped 
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GaAs,  occurring  at  took  temperature,  as  pointed  out  by  Inoue  and  Ohyama.1 
The  electron  mobility  can  be  calculated  from  the  conductivity  expression  given 
in  reference  13,  and  the  result  for  the  present  case  is  found  to  be 
Dn  *  1000  co  /V-sec  at  300*K.  The  slight  difference  in  reslstlrit  ;*■  tJ&tvecn 
#  S-l  and  #  S-2  implies  that  the  chromium  densities  are  not  equal  In  both 
samples.  Presumably,  sample  #S“1  contains  slightly  lower  chromium  Ir^iurity 
density  than  that  of  #  S-2 .  The  Photo-llall  effect  measurements  over  two  samples 
indicate  that  the  Hall  mobility  is  a  slow  varying  function  of  light  intensity 
in  measured  injection  ranges  of  our  PME  and  PC  experiments.  For  sample  #  2, 
it  is  found  mh  ^  150  -  30  cm2 /V-sec  at  300*K  and  UtJ  $  450  ♦  50  cm2/V~sec  at  80*K, 
in  the  intermediate  injection  range. 

The  technique  for  measurements  of  the  PME  and  PC  response  is  described  in 
14 

detail  by  Li  and  will  not  be  repeated  here.  A  high  intensity  tungsten  light 
source  is  used  for  the  PME  and  PC  measurements. 

V.  Results  and  Analysis 

In  Section  IV  wc  have  reported  the  results  of  resistivity.  Hall  effect  and 

photo-Hall  effect  measurements  on  two  samples.  In  this  section,  we  will  present 

the  results  of  PME  and  PC  measurements  at  T  ■  300*K  and  80*K. 

The  PME  open  circuit  voltage  (V^) ,  the  I  and  AC  are  measured  as 

functions  of  incident  radiant  power  densitv,  1  ,  for  T  •  80®K  and  300*K.  The 

o 

results  arc  interpreted  in  :erms  of  the  theory  developed  in  Section  ill. 

Comparisons  between  the  theory  ard  the  experimental  data  are  made  for  both 

samples.  Important  parameters  such  as  t  and  T  are  deduced  from  the  present 

ty  p  r 

results.  It  is  noted  that  the  results  Indicate  that  the  effect  of  hole  trapping 

in  sample  #  S-l  is  less  severe  than  In  sample  #  S-2.  This  Is  consistent  with 

the  data  of  resistivity  and  Hall  effect  measurements  in  which  we  have  found  that 
the  chromium  density  is  slightly  lower  in  sample  #  S-l  than  In  sample  I  S-2. 

Figs.  2  through  4  show  the  results  of  the  PME  and  PC  measurements,  which  will  be 
discussed  subsequently,  as  follows: 

(a)  The  PME  open  circuit  voltage 

The  PME  open  circuit  voltage,  as  a  function  of  magnetic  field,  B, 

for  sample  #  S-l  at  T  »  80*K  and  300*K,  for  different  light  Intensity,  1  is 
shown  in  Fig.  2.  The  linear  dependence  of  V  on  B  at  300*K  for  different 

I  *1  u 

light  Intensity  is  in  accord  with  the  theoretical  prediction  of  Eq.  (13). 
Deviation  from  linearity  at  80r K  is  due  to  the  increase  in  electron  mobility 
with  decreasing  temperatures.  Saturation  of  with  increasing  light  inten- 
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aity  as  shown  in  Fig.  2  provides  Che  evidence  chat  experiments  are  indeed 
perforated  at  high  lnjeccion  regions. 

(b)  The  PME  short  circuit  current  vorsu9  photoconductance  (1^  versus  AG) 

In  order  to  correlate  the  experimental  results  with  the  theory  developed 
in  the  previous  section,  we  plot  1^,  versus  Ac  for  T  ■  80°K  and  300  °K  for 
both  samples,  and  the  result  is  shown  in  Fig.  3.  It  is  interesting  to  note 
that  two-well-defined  regions  arc  obtained  in  the  Ip^g  versus  AG  pi _ts.  In 
region  I  (high  injection),  I_u_  is  directly  proportional  to  AG;  this  is 
observed  for  sample  #  S-l  at  T  ■  300*K.  In  region  II  (intermediate  injection), 
I  varies  with  this  is  observed  tor  both  samples  at  80°K  and  for  sample 

s  *1  Li 

#  S-2  at  300*K.  To  interpret  this  result.,  the  generalized  theory  enveloped 

in  the  previous  section  is  used  for  this  purpose. 

The  high  injection  linear  region  observed  for  sample  #  S-l  at  300°K  is 

in  accord  with  the  theoretical  prediction  of  Eq.  (11).  The  electron  and  hole 

lifetimes  are  equal  and  independent  of  light  intensity.  In  the  intermediate 

injection  region  (region  II),  the  effect  of  hole  trapping  at  the  chromium 

-1/3 

acceptor  levels  and  the  dependence  of  T  on  An  leads  to  the  observed 

i  IC  ^ 

dependence  of  Ipjjg  on  Ac  for  both  samples  at  80°K  and  for  sample  t  S-2  at 

300*K.  This  is  predicted  by  Eq.  (13),  derived  in  the  previous  section  under 

-1/3 

the  assumption  that  T  »  KAn  and  Ap  *  TAn.  The  dependence  of  carrier 

0  1 
lifetimes  on  the  injected  carrier  densities  has  been  cited  by  several  authors, 

although  no  explicit  form  has  so  far  been  proposed  in  the  literature.  In  fact, 

g 

the  powerlav  dependence  of  in  on  An  can  be  established  for  different  injection 
ranges.  The  values  of  8  may  vary  from  material  to  material  and  the  ranges  of 
injection  as  well.  (Li1  has  observed  varies  with  An  ^  in  silicon  at  low 
temperatures) . 

(d)  Carrier  lifetimes  versus  photoconductancc 

We  have  shown  in  the  previous  section  that  the  lifetimes  of  holes  and 
electrons  are  equal  and  independent  of  light  intensity  under  high  injection 
conditions.  This  is  observed  for  sample  #  S-l.  However,  in  the  intermediate 
injection,  the  electron  and  hole  lifetimes  are,  in  general,  functions  of  injec¬ 
tion.  In  order  to  determine  the  electron  and  hole  lifetimes  as  functions  of 
injection,  it  is  necessary  to  define  an  apparent  lifetime  T^,  in  the  IpME 
expression.  This  can  be  achieved  by  rewriting  Eq.  (13)  in  the  form  that 


11 


16 


where 


T 

a 


PME 


(*£) 

V25; 


Da  ** 

(-“)  B  AG 
a  a 
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(15) 


is  the  PME  apparent  lifetime.  Eq.  (14)  resembles  Eq.  (10)  for  the  high 

injection  case,  with  t  replaced  by  x  . 

n  ■'a 

Eq.  (15)  predicts  that  the  apparent  lifetime  x  is  proportional  to 
the  ZiG  ^  in  the  intermediate  injection  range. 

With  the  help  of  Eq.  (14)  and  Fig.  3,  the  apparent  lifetime  x  is 

3 

computed  as  a  function  of  AG  and  the  results  are  shown  in  Fig.  4.  The 

result  of  Xfl  versus  tC  plot  yields  a  slope  of  -  ”,  which  is  in  accord  with 

the  theoretical  prediction  of  Eq.  (15).  Note  that  x  is  equal  to  x  and  T 

a  n  p 

in  the  high  injection  region. 

In  the  intermediate  injection  region  where  trapping  of  holes  at  the 
chromium  acceptor  levels  is  essential,  the  electron  lifetime  is  expected  to 
be  longer  than  the  hole  lifetime  (i.c.,  T  <  1).  This  is  because  the 
negatively-charged  chromium  acceptor  levels  have  a  larger  capture  cross 
section  for  holes  than  for  electrons.  To  deduce  the  electron  and  hole  life¬ 
times  from  the  PME  apparent  lifetimes,  we  make  use  of  the  following  relation 
in  the  presence  of  trapping 


T 

a 


Apx  +  Anx 

J .  n„ .  -  . P 

An  +  Ap 


+  Tx 

n 

i  +  r 


A 

VS 


(16) 


where  wc  assume  F  =  x  h  <  1. 

p  n 


In  essence,  the  PME  apparent  lifetime,  X  is  a  parallel  sum  of  x  and  x 

a  n  p 

and  is  controlled  by  the  shorter  of  the  two  lifetimes,  namely,  x  .  With  the 
help  of  Eq.  (15),  x^  is  deduced  from  Eq.  (16)  and  the  results  are  shown  in 
Fig.  4. 


The  remainder  of  the  parameters  x  ,  T  and  <  can  be  deduced  subsequently 

by  using  Eqs.  (15),  (2-b)  and  Fig.  3,  and  the  relation  that  An  =  AG/qu  d 

n 

(where  d  is  the  thickness  of  the  sample).  It  is  worth  noting  that  the 
computed  values  of  T  and  <  are  fairly  constant  over  the  entire  intermediate 
injection  range.  The  value  of  T  is  found  equal  to  0.18  at  300°K  and  0.084 
at  80° K  for  sample  8  S-2,  which  shows  that  the  ratio  of  x  /x  is  about  5.5 

n  p 

at  300  K  and  is  about  12  at  80° K  for  sample  8  S-2.  The  results  are  in  good 
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agreement  with  the  statement  by  Holeman  and  Hilsum  in  that  they  predicted  a 
ratio  of  <  10  for  semi-insulating  GaAs.  The  present  results  show  that 

the  effect  of  hole  trapping  in  the  chromium  levels  is  not  too  severe  even  at 
80  K  in  the  intermediate  injection  region,  and  can  be  ignored  at  high  injection. 

The  decrease  of  electron  and  hole  lifetimes  (see  Fig.  4)  with  increasing 
injection  indicates  that  more  and  more  chromium  acceptor  centers  are  converted 
from  the  originally  hole-trapped  centers  to  effective  recombination  centers  for 
electrons  as  injection  is  increased.  This  is  consistent  with  the  recombination 
model  proposed  in  Fig.  1. 

VI.  Conclusions 

Photomagnetoelectric  (PME)  and  photoconductive  (PC)  effects  in  semi- 

insulating  chromium-doped  GaAs  have  been  investigated  at  T  =  80°K  and  300°K. 

Using  a  simple  recombination  and  trapping  model,  a  generalized  theory  has  been 

developed  for  the  PME  and  PC  effects,  taking  into  account  the  effect  of  trapping 

and  the  dependence  of  the  carrier  lifetimes  on  the  injection,  under  large 

injection  conditions.  Good  agreement  between  theory  and  experiment  is  obtained 

for  two  well-defined  injection  ranges;  one  region  is  characterized  by  a  linear 

dependence  of  I_.._  on  AG  and  constant  lifetime  (t  =  t  =  2.5  x  10  ^  sec.  at 
PME  n  p 

300°K  for  sample  //  3-1),  and  the  other  region  is  characterized  bv  a  nonlinear 

1  2 

dependence  of  IpME  on  AG  (i.e.,  Ip^E  varies  as  AG  ’  ).  Numerical  values  of  xn» 
Tp,  AG  and  Ip^E  as  functions  of  injection  are  given  for  both  samples. 
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Fig.  1 


Energy  band  diagram  for  Cr-dopcd  CnAs  in  thermal  equilibrium  and 
under  steady  state  illumination 


Fig.  2  The  PME  open  circuit  voltage,  V  versus  magnetic  flux  density,  b,  for 
different  light  intensities 

a)  IQ  =  100  mW/cm  ,  b)  40  mU'/cm^,  c)  22  mW/cm^,  d)  11  nV/cm^, 

2  n 

c)  4  mW/cm  ,  f)  2  inW/cm- 


a')  I  =  256  mW/cm 


b ')  I  =  100  tnW/cm^,  c1)  I  =  40  r.V'/cra^, 
o  o 


Fig.  3  The  PME  short,  circuit  current  per  unit  width  of  sample,  I  „,  versus 

photoconductance,  AG,  for  samples  US-1  and  fS-2,  at  S0°K  arid  300°K.  In 
region  (1),  I_  is  proportional  to  AG,  and  in  region  (II),  In>,,.  varies 

"ith“  |B=0  1wb/^).  PL 


The  PME  apparent  lifetime  ,  the  electron  lifetime,  t  and  the  hole 
lifetime  T  versus  AG  in  the  intermediate  injection  region  (region  II) 
for  sample  v‘S-2  at  80  k  and  300°K.  The  slope  of  t  versus  AG  is  -  2/5. 


Fig.  4 


C.  LOW  TEMPERATURE  PHOTOMAGNETOELECTRIC  AND  PHOTOCONDUCTIVE  EFFECTS  IN 
N-TYPE  InAs  (S.S.  Li  and  C.I.  Huang) 

Introduction 

The  room  temperature  photomagnetoelectric  (PME)  effect  in  both  n-  and 

p-type  InAs  was  reported  previously  by  Dixon.1  The  results  were  interpreted 

in  terms  oi  the  large  Hall  angle  PME  theory  for  p-type  samples  and  of  the 

3 

small  Hall  angle  PME  theory  for  n-type  samples. 

In  the  present  work  we  extended  the  experimental  study  of  the  PME  and  PC 
effects  in  single  crystal  n-type  InAs  from  room  temperatures  down  to  4.2°K. 
Incorporating  the  two  measurements  with  resistivity  and  Hall  effect  experiments 
would  allow  us  to  determine  the  carrier  lifetimes  as  well  as  electron  and  hole 
mobilities  at  low  temperatures.  The  carrier  lifetimes  were  computed  from  the 
PME  and  PC  data,  and  the  electron  and  hole  mobilities  were  deduced  from  the 
Hall  effect  and  the  PME  data  between  4.2°K  and  300 °K. 

Experimental 

Undoped  n-type  InAs  grown  by  the  Czochralski  method  with  electron  con- 
16  -3 

centration  of  2.4  x  10  cm  was  used  in  the  present  work.  The  sample  was  cut 
into  rectangular  form  by  using  a  diamond  saw.  The  dimensions  of  the  sample, 
after  being  mechanically  lapped  by  SiC  powder  and  chemically  etched  by  CP^, 
were  1.3  x  0.4  x  0.04  cm^.  The  sample  was  then  mounted  onto  the  cryotip  of 
the  AC-310  refrigeration  system.  The  resistivity  and  Hall  effect  measurements 
were  performed  by  using  the  a-c  method,  similar  to  that  described  in  Reference  1. 
The  tungsten  lamp  light  source  was  chopped  at  400  cps,  and  the  PME  and  PC 
signals  picked  up  from  the  sample  were  passed  through  amplifiers  and  a  wave 
analyzer.  The  system  provides  a  gain  of  10^,  which  would  allow  uj  to  measure 
a  relatively  small  signal. 

Results  and  Discussion 

The  resistivity  and  Hall  effect  were  measured  between  4.2°K  and  300°K. 

The  results  showed  that  the  electron  concentration  was  nearly  constant  over  the 
entire  temperature  range,  and  the  Hall  mobility  (assuming  equal  to  electron 
mobility)  reached  a  peak  value  at  77°K.  The  results  are  summarized  in  Table  I. 

The  PME  open-circuit  voltage  and  the  photoconductance  as  functions  of 
magnetic  induction  and  light  intensity  were  measured  between  300°K  and  4.2°K. 

The  results  for  T  >  77°K  were  in  accord  with  the  small  Hall  angle  PME  theory, 
and  for  T  <  21°K,  the  large  Hall  angle  PME  theory  predicted  the  observed  PME 
data.  This  will  be  discussed  separately  in  the  following  subsections: 

A.  Results  for  77°K  <  T  <  300°K. 

In  this  temperature  range,  the  small-signal  PME  and  PC  data  are  shown  in 

Figs.  1  and  4.  Fig.  1  shows  the  PME  open-circuit  voltage  V_._  versus  the 

PME 
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as  a 


magnetic  induction,  B,  for  T  =  300°K  and  77 °K,  with  light  intensity,  1 
parameter.  The  results  are  in  accord  with  the  small  Hall  angle  PME  theory 

3 

given  by  Van  Roosbroeck,  which  reads 


PME 


»  £ 
P  0 


(1) 


where  D  is  the  hole  diffusion  coefficient,  x  is  the  hole  lifetime,  G  is 
P  p  o 

the  dark  conductance  and  AG  is  the  photoconductance.  The  linear  relationship 

between  Vp^E  and  B,  observed  at  77 °K  and  300°K,  indicates  that  the  hole 

mobility  is  low  in  this  temperature  range.  Assuming  that  the  electron  mobility 

equals  the  Hall  mobility,  and  electron-hole  mobility  ratio  of  100,  the  holi 

mobility  was  estimated  at  300°K  and  77 °K  (see  Table  I).  The  hole  lifetime 

(which  is  equal  to  electron  lifetime  for  the  present  case)  is  deduced  from 

Eq.  (1)  and  the  PME  data  shown  in  Fig.  4.  The  carrier  lifetime  estimated  at 

-8 

300°K  (t  =  3.7  x  10  sec)  is  in  good  agreement  with  that  reported  previously 
by  Dixon.  The  substantial  increase  in  the  photoconductance  at  77°K  (as 
compared  with  value  at  300°K)  is  due  to  the  fact  that  both  the  carrier  mobility 
and  the  carrier  lifetime  are  Increased  substantially  at  77 °K. 

B.  Results  for  4.2°K  <  T  <  21°K 

In  this  temperature  range,  the  VpME  saturates  at  high  magnetic  field 
strength,  indicating  that  the  hole  mobility  is  greatly  increased  at  low  temper¬ 
atures.  To  interpret  the  observed  VpM£  versus  B  plot  as  shown  in  Fig.  2,  the 
large  Hall  angle  PME  theory  is  used.  The  theory  predicts  that^ 


PME 


D  H 
=  (-*) 


(1 


+  W2bV 


AG 

G 


(2) 


where  y  «*  (y^y^)  Is  the  effective  carrier  mobility.  Eq.  (2)  predicts  that 

the  Vp^E  is  proportional  to  B  at  low  magnetic  field  and  is  independent  of  B 

at  high  magnetic  field.  The  results  shown  in  Fig.  2  for  T  =  21°K  and  4.2°K 

are  in  accord  with  the  prediction  given  by  Eq.  (2).  To  estimate  effective 

mobility  y  and  hole  lifetime  T  from  Eq.  (2)  and  Fig.  2,  a  plot  of  (B/V  M„)^ 

2  P  rWCi 

versus  B  is  shown  in  Fig.  3.  The  intercept  of  the  straight  lines  snown  in 

Fig.  3  with  the  ordinate  yields  the  hole  lifetime  and  the  slope  of  these  lines 

yields  the  effective  carrier  mobility.  The  computed  values  of  y  and  tp  at 

21° K  and  4.2°K  from  Fig.  3,  are  summarized  in  Table  I.  The  electron  mobility 

is  deduced  from  resistivity  and  Hall  effective  measurements,  and  the  hole 
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mobility  is  computed  from  Mp  -p2/pn.  The  results  show  that  the  electron  and 
hole  mobility  ratio  changes  drastically  from  assumed  value  of  1001  at  300°K 
to  the  value  of  b  -  5  (computed  value)  at  4.2°K.  This  is  indeed  consistent 
with  the  observed  VpME  versus  B  data  for  T  at  300°K  and  4.2°K.  The  change  in 
mobility  ratio  with  temperature  has  also  been  observed  in  InSb.^ 

Fig.  4  illustrates  the  plot  of  versus  AG/Gq  for  T  »  300°K,  77°K, 

21°K  and  4.2°K.  The  results  (slope  of  this  plot  is  one)  at  all  temperatures, 
show  that  no  trapping  effect  is  involved  in  the  entire  temperature  range.6,7 
It  is  noted  in  Fig.  3  that  the  deviation  from  linearity  (as  is  predicted  by 
Eq.  5)  at  higher  magnetic  field  is  due  to  the  fact  that  the  photoconductance 
becomes  magnetic  field  dependence  as  B  is  increased. 

Conclusions 

Experimental  study  of  the  PME  and  PC  effects  in  n-type  InAs  single  crystals 
has  been  extended  from  room  temperature  down  to  4.2°K.  The  results  in  good 
agreement  with  the  theory  developed  previously.  Carrier  lifetimes  and  mobilities 
are  deduced  from  the  present  results.  Measurements  of  the  PME  and  Hall  effect 
at  low  temperatures  enable  us  to  determine  both  the  majority  and  minority  carrier 
mobilities  in  n-type  InAs.  Study  of  the  PME  and  PC  effects  show  that  no  trapping 
is  involved  in  the  entire  temperature  range  for  the  present  samples. 
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Table  I  Transport  and  Recombination  Parameters  for  n-type 
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Figure 


Fig.  1 


Fig.  2 


Fig.  3. 


Fig.  4. 


Ions 

PML  open  circuit  voltnge  per  unit  sample  length, 

VPHE *  vcrsus  “agnetic  Induction,  B,  with  light  intensity  T  ns 
a  parameter,  for  T  -  300*K  and  77°K. 

Curve  a)  Iq  -  430  mW/ca2  b)  270  mK/cm2  c)  81  mW/cm2 

d)  430  mW/cm2  e)  224  mW/cm2  f)  129  mW/cn2 

PME  open  circuit  voltage  per  unit  sample  length,  V  , versus 

rMc 

esagne  :1c  induction,  B,  with  light  intensity  I  as  a  parameter, 
for  T  «  21*k  and  4.2®K. 

Curve  a)  I  -  430  mW/cm2  b)  224  mW/cn2  c)  81  mW/cm2 

°  2  1  o 

d)  350  eW/cm  e)  156  mW/cm  f)  106  mW/cm 

2  2 

versus  B  with  light  intensity  Iq  as  a  parameter  for 

'  *  21*K  a^d  4.2®K.  The  intercept  of  this  plot  with  ordinate  is 

equal  to  (r^)  and  the  slope  of  this  plot  yields  (t^J*2. 

.  u  _ 

Curve  a)  Ip  -  106  mW/cm  b)  224  mW/crn  c)P  430  mW/cm2 

d)  106  mW/cm2  e)  224  mW/cm2 

PME  open  circuit  voltage,  Vp^,  versus  photoconductance  ratio, 

For  T  ■  300  K,  77  K,  21®K  and  4.2°K.  The  slope  of  these  plots  ° 
is  unity. 
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D.  THE  EFFECT  OF  A  BURIED  LAYER  ON  THE  COLLECTOR  BREAKDOWN  VOLTAGES  OF 
BIPOLAR  JUNCTION  TRANSISTORS  (F.W.  Hewlett,  Jr.,  F. A.  Lindholm  and 
A.  J.  Brodersen) 


I.  Introduction 

Bipolar  transistors  designed  for  high-frequency  or  high-speed 
applications  are  very  shallow  structures.  For  example,  devices  used  for 
small-signal  microwave  amplification  may  have  emitter  and  collector 
junctions  that  lie  withia  5000  Angstrom  units  from  the  surface;  typically, 
the  devices  are  fabricated  by  successive  diffusions  into  N/N+  epitaxially 
grown  silicon  [1] .  The  N-type  epitaxial  material  is  made  thin,  less  than 
10  microns,  to  reduce  collector  series  resistance.  Similar  devices  are 
used  commonly  in  high-speed  digital  integrated  circuits. 

In  this  paper,  we  describe  a  premature  punch-through  effect  occurring 
in  this  type  of  transistor  that  can  lower  the  collector  breakdown  voltages 
below  the  values  predicted  by  existing  theory.  We  express  the  premature 
punch-thiough  voltage  BVppTi  which  characterizes  the  effect,  in  terms  of 
parameters  of  the  fabrication  processing,  and  from  this  derive  an  expression 
for  the  thickness  of  the  epitaxial  collector  that  maximizes  BVppT  and,  hence, 
the  usable  region  of  device  operation.  Measurements  are  presented  that 
support  the  theoretical  predictions. 

Punch-through  occurs  when  the  space  charge  layer  of  a  reverse  biased 
PN  junction  contacts  the  space  charge  layer  of  a  nearby  PN  junction.  Reach- 
through  occurs  when  the  depletion  region  of  a  reverse  biased  junction  contacts 
the  boundary  between  the  P-type  or  N-type  material  and  that  of  a  nearby  layer 
of  the  same  but  higher  conductivity  type.  The  punch- through  condition  may 
occur  in  shallow  bipolar  transistors  whose  doping  in  the  collector  is 
comparable  with  that  in  the  base.  The  reach-through  condition  may  occur 


/ 


in  PN  junctions  fabricated  on  high  conductivity  buried  layers.  In  this  case, 
premature  avalanche  breakdown  due  to  reach-through  takes  place. 

In  the  premature  punch-through  effect,  a  sufficiently  large  reverse 
Dias  placed  on  the  collector  base  junction  can  cause  the  collector  base 
depletion  region  to  reach  through  the  lightly-doped  collector  to  the  N+ 
buried  layer.  Further  increase  of  the  collector  base  voltage  acts  to  widen 
the  collector  depletion  region  mainly  on  the  base  side,  forcing  it  to  punch 
through  the  base  into  the  emitter  space  charge  region.  This  occurrence 
results  in  a  large  current  flow,  or  breakdown,  if  che  emitter  and  base  are 
short-circuited. 

In  Section  II,  we  summarize  the  existing  theory  for  punch-through  [2] 
and  reach— through  [3]  and  present  a  theory  which  is  new  for  premature  punch- 
through.  Section  III  contains  experimental  evidence  which  substantiates 
the  theory  presented  in  Section  II. 

II.  Theory 

A.  PUNCH-THROUGH  IN  TRANSISTORS 

In  this  part,  we  present  a  method,  discussed  in  Phillips  [2],  which 
enables  one  to  distinguish  between  breakdown  resulting  from  punch-through 
in  the  base  and  that  resulting  from  avalanche  multiplication  in  the  collector. 
The  discussion  herein  is  qualitative  and  yields  relations  between  the  break¬ 
down  voltages  of  a  bipolar  transistor  valid  for  both  punch-through  and 
premature  punch-through. 

If  breakdown  is  due  to  avalanche  multiplication  in  the  collector  space 
charge  region,  we  may  determine  this  as  follows.  One  measures  the  breakdown 
voltage  from  emitter  to  collector  with  the  base  shorted  to  the  emitter, 

BVCES’  ant*  tlie  breakdown  voltage  from  collector  to  base  with  the  emitter 
open  If  avalanche  multiplication  in  the  collector  is  causing  break¬ 

down,  then 


BV 


CES 


BV 


CBO ' 


(1) 


Thus,  avalanche  breakdown  and  the  ensuing  heavy  current  do  not  depend  on 
the  constraint  applied  at  the  emitter  lead.  Large  currents  flow  at  the 
same  value  of  base  collector  voltage  both  when  the  emitter  is  shorted  to 
the  base  and  when  the  emitter  is  open  circuited. 

Next,  consider  a  device  in  which  the  breakdown  voltage  is  punch-through 
limited.  Here,  the  breakdown  voltage,  BVCB0,  is  related  to  the  punch-through 
voltage  BVpT,  by 
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(2) 


where  BV_,Dn  is  the  breakdown  voltage  applied  between  emitter  and  base 

CiDU 

with  the  collector  open  circuited.  This  relation  may  be  explained  as  follows. 

As  the  collector  base  reverse  bias  is  increased;,  the  collector  transition 

region  punches  through  into  the  emitter  space  charge  region.  If  the  emitter 

is  open  circuited,  this  condition  will  not  result  in  a  large  current  flow. 

A  further  increase  in  voltage  is  necessary  to  cause  avalanche  breakdown  of 

the  emitter-base  sidewall,  which  then  results  in  the  flow  of  heavy  current 

from  the  base  lead.  Thus,  for  a  device  limited  by  punch-through  is 

equal  to  the  voltage  necessary  to  cause  the  collector  transition  region  to 

punch  through  the  base  plus  the  breakdown  voltage  of  the  emitter  base  sidewall. 

In  the  measurement  of  one  shorts  the  emitter  to  the  base.  In  this 

case,  when  punch  through  occurs,  we  have  two  layers  of  N-type  material 

separated  by  a  layer  of  space  charge.  It  is  not  necessary  to  breakdown  the 

emitter  side  wall  since  the  shorted  emitter  lead  provides  a  current  path. 

The  energy  band  diagrams  for  three  values  of  Vcfi  are  shown  in  Figure  1. 

Figure  1(a)  represents  a  thermal  equilibrium  situation.  The  energy  band 

of  Figure  1(b)  illustrates  a  reverse  biased  collector  junction.  Figure 

lCc)  shows  the  punch- through  condition.  For  VCfi  -  BVpT  the  potential  barrier 

has  been  reduced  so  that  electrons,  supplied  by  the  shorted  emitter  lead, 

can  overcome  the  barrier  and  slide  down  the  potential  "hill."  As  V  is 

CB 

increased  above  BVpip,  space-charge-limited  flow  ensues  [4], 

Therefore, 


bvces  =  Bvpr 


(3) 


for  a  transistor  in  which  punch  through  determines  the  breakdown  characteristics; 
and,  from  equations  (2)  and  (3), 


bvces  +  BW 


Phillips  [2]  gives  a  theoretical  expression  for  BVpT  in  terms  of 
fabrication  parameters  applicable  to  a, uniform  base  alloy  transistor: 


(5) 


where  is  the  acceptor  concentration  in  the  base  region,  and  Wfi  is  the 
width  of  the  base  region.  Figures  2(b)  and  2(c)  illustrate  the  electronic 
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(a)  Thermal  equilibrium. 


E 


Figure 


(b)  The  collector  reverse  biased. 


(c)  Onset  of  space  charge  limited  flow. 


1  Energy  band  diagram  for  an  NPN  transistor  with 
various  applied  collector  voltages.  Only  the 
conduction  band  edge  is  shown. 


charge  and  electrostatic  field  distributions  at  the  punch-through  condition. 
Equation  (5)  assumes  that  the  collector  base  junction  is  asymmetrical,  the 
collector  being  much  more  heavily  doped  than  the  base.  Therefore, 
approximately  all  the  voltage  difference  appears  on  the  base  side  of  the 
collector  transition  region.  The  collector  base  voltage  at  punch-through 
is  approximated  by  the  area  of  the  shaded  triangle  in  Figure  2(c). 

For  a  double  diffused  structure,  one  may  obtain  an  approximate  value 
for  BVp^,  from  the  Lawrence  and  Warner  curves  [5]  ,  using  an  iterative 
procedure.  However,  when  space  charge  region  widths  become  comparable 
with  the  quasi-neutral  base  width,  a  computer  solution  may  be  necessary. 

B.  REACH-THROUGH  IN  DIODES 


Reach-through  is  an  effect  which  can  limit  the  breakdown  voltage  of 
P+N  junction  diodes  fabricated  on  N/N+  epitarial  material,  of  the  type 
shown  in  Figure  3(a).  An  increase  of  the  reverse  voltage  widens  the 
transition  or  space-charge  region,  which  extends  primarily  into  the  lightly 
doped  N  region.  For  a  particular  value  of  reverse  bias,  the  space  charge 
region  extends  to  the  buried  layer.  Further  increase  in  reverse  voltage 
causes  uncovering  of  immobile  ion  charge  in  the  N+  layer,  as  shown  in 
Figure  3(b).  This  produces  a  sensitive  dependence  on  reverse  voltage  of 
the  maximum  electric  field  E^^,  which  grows  until  It  reaches  the  critical 
field  for  avalanche  multiplication  This  critical  field  supplies 

conduction  electrons  with  sufficient  energy  to  impact  ionize  valance 
electrons,  resulting  in  heavy  current  flow.  The  charge  and  electrostatic 
field  distributions  at  breakdown  are  illustrated  in  Figures  3(b)  and  3(c). 

Grove's  expression  for  the  breakdown  voltage  due  to  reach-through 
BVrt  in  terms  of  the  critical  field  and  fabrication  parameters  is  [3] 


RV  =  F  W 

RT  CRIT  EPI 


qNp  wepi 

2Ks  eo 


(6) 


Here  is  the  donor  concentration  in  the  epitaxial  silicon  and  i-s 

the  width  in  epitaxial  material  measured  from  the  metallurgical  P  N  step 
junction  to  the  buried  layer.  Equation  (6)  is  an  approximation  for  the 
shaded  area  indicated  in  Figure  3(c).  It  is  a  good  approximation  for  a 
single  diffused  p"*"nN+  diode. 

This  development  assumes  a  step  function  change  in  collector  doping 
as  one  moves  from  the  lightly  doped  collector  toward  the  buried  layer. 
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N 


(a)  Intrinsic  structure. 


(b)  Electronic  charge  distribution 


E 


EMAX  "  ecrit 


(c)  Electrostatic  field  distribution. 

Figure  3  Illustration  of  a  P+NN+  diode  in  the 
reach-through  condition. 
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In  reality,  the  out-diffusion  of  the  substraight  impurity  during  epitaxial 

growth  and  subsequent  heat  treatments  may  cause  BV  to  be  less  than  the 

KT 

value  predicted  by  equation  (6).  A  similar  reduction  of  the  breakdown 
voltage  resulting  from  the  out-diffusion  of  the  substrate  impurity  also 
holds  for  the  premature  punch-through  effect  in  Junction  transistors, 
which  we  are  now  prepared  to  examine. 

C.  PREMATURE  PUNCH-THROUGH  IN  TRANSISTORS 

The  breakdown  of  the  P^NN*  diode,  illustrated  in  Figure  3  results 
from  avalanche  multiplication  in  the  P+N  junction  transition  region.  The 
maximum  field,  attains  the  value  of  the  critical  field  prematurely 

because  of  the  presence  of  the  buried  layer.  The  presence  of  another 
transition  region  creates  an  alternative  mechanism  for  the  occurrence  of 
large  current  flow.  This  possibility  exists  in  NPNN+  bipolar  structures. 

For  simplicity  of  discussion,  we  deal  first  with  a  one  dimensional  structure 
with  homogeneous  base  doping  which  is  illustrated  in  Figure  4(a). 

In  this  type  of  shallow  device,  the  emitter  and  collector  transition 
regions  are  separated  by  a  narrow  distance,  typically  less  than  om  micron. 
Thus,  the  collector  space  charge  region,  broadening  with  incrca«dr$, .reverse 
voltage,  may  punch  through  the  quasi-neutral  base  into  the  emitter  space 
charge  region  before  the  critical  field  for  avalanche  multiplication  in 
the  collector  is  reached.  The  charge  and  electric  field  distributions  are 
Illustrated  in  Figures  4(b)  and  4(c),  respectively,  for  a  transistor  with 
uniform  doping  in  the  collector  and  base.  The  breakdown  voltages  of  a 
structure  experiencing  premature  punch-through  obey  equation  (4). 

To  obtain  an  expression  for  the  breakdown  voltage  due  to  premature 
punch-through,  BVppj,  we  must  Include  the  voltage  drops  in  both  the  base 
and  the  collector  as  represented  by  the  shaded  area  of  Figure  4(c).  This 
area  may  be  expressed  as  the  sum  of  the  areas  of  Figure  2(c)  and  3(c), 
which  corresponds  to  replacing  ECRIT  by  E^  in  equations  (5)  and  (6). 

Thus,  an  approximate  expression  for  the  breakdown  voltage  due  to  premature 
punch-through  may  be  written  as 


where 


bvppt  "  FMAX  wepi 


W2 


EPI 

0 


+ 


w‘ 


(7) 


EMAX 


<na  wb 

Ks  co 
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EMIT 


(a)  Intrinsic  structure. 


(b)  Electronic  charge  distribution 


(c)  Electrostatic  field  distribution. 


Figure 


4  Illustration  of  an  N+PNN+  transistor  in  the 
premature  punch-through  condition. 


i 
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when  the  punch-through  condition  occurs. 

This  development  applies  for  a  homogeneous  base  transistor.  However, 
one  may  approximate  the  value  of  BVpp,p  for  a  double  diffused  device  by 
combining  equations  (7)  and  (8), 


BV 


PPT 


KS  eo 


(l\WBWEPI 


2  2 
N  W  N  W 

D  EPI  +  NAWB  ) 


(9) 


and  by  noting  thac  the  product  of  N.  and  WD  v:  *<  ,i  appears  in  equation  (9) 

A  b 

is  Gummel's  number  N„  for  a  homogeneous  base  transistor.  Therefore,  as  an 

b 

approximation  to  BVppT  for  the  double  diffused  structure  we  write 


BV 


PPT 


KS  £0 


(nbwepi 


N  N  W 

D  EPI  B  B 

2  +  2  J 


(10) 


where  N  may  be  determined  via  terminal  measurements  [6]  and  W  may  be 
b  b 

determined  by  the  groove  and  stain  methc  1  (7). 

Equation  (10) indicates  that,  for  a  device  whose  breakdown  voltages 

are  limited  by  premature  punch-through,  there  is  a  value  of  W  T  which 

hr  I 

gives  a  maximum  value  of  BVppT«  Differentiating  equation  (10),  we  find 

the  optimum  epitaxial  thickness,  (W  )  ,  is 

EPI  OPT 


<WEPX> 


OPT 


V 


(11) 


Thus,  ®VppT  will  be  the  largest  for  a  particular  fabrication  process  when 

the  epitaxial  thickness  is  equal  to  Gummel's  number  (in  atoms/cm^)  divided 

3 

by  the  doping  concentration  of  the  lightly  doped  collector  (in  atoms/cm  ). 

The  next  section  contains  data  on  bipolar  structures  that  experience 
premature  breakdown  resulting  om  reach-through  and  from  punch-through 
induced  by  reach-through.  For  comparison,  it  also  shows  the  breakdown 
voltage,  of  a  related  device  which  experiences  neither  effect. 

III.  Experimental 

We  have  designed  and  fabricated  two  runs  of  bipolar  transistors 
whose  fabrication  parameters  and  profiles  were  similar  but  whose  values 
of  and  BV^^  were  considerably  diff^'  '*nt.  The  feature  which 

distinguished  the  two  runs  was  that  one  of  them  employed  a  buried  layer. 

These  devices  were  designed  to  avoid  punch  through  of  the  base  in  the 
usual  sense.  The  purpose  of  the  experiment  was  to  investigate  the  effect 
of  a  buried  layer  on  the  breakdown  characteristics  of  these  NPNN+  transistors. 
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Finnic  5(a)  Photograph  of  HV 


for  test  unit 


Figure  6(a)  Photograph  of  BV  for  test  unit 
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Figure  5(b)  Photograph  of  BV  for  test  uiit  ti  1  -<1I 


Figure  6(b)  Photograph  of  BV 


for  test  unit  #3-41 
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Data  taken  from  the  bipolar  transistors  fabricated  are  contained  in 

Table  1.  Test  unit  #1-4L  was  fabricated  on  N-type  material  without  a 

buried  layer.  For  this  device  BV„_.  is  equal  to  BV„,,  (equation  (1)) 

which  indicates  that  avalanche  multiplication  in  t’  e  collector  transition 

region  causes  the  large  current  flow.  BV___  and  BV„_.  for  this  device 

are  pictured  in  Figure  5  (a)  and  (b),  respectively.  The  other  test  units 

in  Table  1  were  fabricated  on  N/N+  epitaxial  material.  Test  units  #3-2L 

and  It 3-3L  exhibit  breakdown  due  to  premature  avalanche  multiplication  in 

the  collector  transition  region  caused  by  reach-through.  BV-,,,.,  is  equal 

to  for  these  test  units.  Test  units  It 3-4L  and  It 3-6L  are  limited  by 

premature  punch-through.  Figure  6(a)  shows  a  photograph  of  BV,„,0  for 

CEb 

test  unit  //3-4L.  Figure  6(b)  shows  a  photograph  of  BV  for  this  test 

CBO 

unit  which  is  greater  than  BVCES  by  the  value  of  BV£B0  obeying  equation  (4). 
The  values  of  and  IEg  were  determined  by  measurement  using  the  method 
described  by  Gummel  [6] . 

Test  unit  It 3-3L,  which  has  a  higher  value  of  N  and  hence  more 

+12  -2 

impurities  in  the  base  layer  per  unit  area,  5.6  x  10  cm  ,  than  test 

+12  -2 

unit  i!‘3-2L,  4.7  x  10  cm  ,  breaks  down  in  the  collector  at  lower 

voltages.  However,  It 3-6L,  which  has  a  larger  value  of  N  than  It 3-4L, 

breaks  down  at  correspondingly  higher  values  of  BVCEg  and  BVCBQ.  Thus,  the 

addition  of  more  impurities  in  the  base  increases  BV  ■  BV,,,,™  for  test 

CEb  PPi 

units  It 3-4L  and  It 3-6L,  in  which  premature  punch  through  dominates,  and 
decreases  BV^g  for  units  It 3-3L  and  It 3-2L,  in  which  the  dominant  mechanism 
is  avalanche  multiplication  induced  by  reach-through. 

IV.  Summary  and  Conclusions 

We  have  combined  existing  theory  for  punch-through  and  reach-through 
to  obtain  an  expression  for  the  premature  punch-through  breakdown  voltage, 
equation  (10).  Device  characteristics  which  show  premature  punch  through 
are  presented  and  compared  with  those  which  show  the  usual  avalanche 
breakdown  in  the  collector. 

Equation  (11)  gives  in  terms  of  fabrication  parameters  the  thickness 

of  the  epitaxial  layer  which  maximizes  the  breakdown  voltages  (both  BV„_,C 

CES 

and  BVebq)  of  a  device  limited  by  premature  punch- through. 

Experimental  results  indicate  that  large  variations  in  the  breakdown 
voltages  of  devices  fabricated  on  the  same  thin  epitaxial  layer  may  be 
expected  (37  to  50  volts).  These  valuations  likely  result  from  two  causes; 
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TABLE  1 


BREAKDOWN  VOLTAGES 


Test  Unit 

bvces 

(volts) 

DVCBO 

(volts) 

BV 

EBO 

(volts) 

nb 

cm-2 

#3-2L 

48. 

48. 

4.7  X  10+12 

#3-3L 

46. 

46. 

5.6  X  10+12 

A3-4L 

37. 

45. 

8. 

2.3  x  10+12 

#3-GL 

42. 

50. 

8. 

3.1  x  10+12 

#1-4L 

75. 

75. 

2.9  x  10+12 

48 


(1)  because  of  nonuniformities  in  the  diffusion  process  and  the  fact  that 
the  base  widths  were  determined  by  the  tails  of  two  diffusions  result  in 
variations  in  the  number  of  impurities  in  the  base,  and,  therefore,  in  the 
breakdown  voltages;  (2)  because  of  nonuniformities  in  the  epitaxial  layer 
thickness,  Wgpp.  More  important,  the  fact  that  both  premature,  punch- 
through  and  reach-through  occur  in  different  devices  taken  from  the  same 
wafer  causes  wide  variations  in  the  values  of 

Gummel's  number  is  interpreted  in  equation  t  lO)  as  the  net  number 
of  impurities  per  unit  area  in  the  base  layer In  certain  shallow 
transistors,  this  interpretation  is  no  longer  valid  [8],  In  these  devices, 
Ng  is  the  number  of  holes  per  unit  area  in  the  transistor  [8,9].  Therefore 
we  would  not  expect  equation  (10)  to  accurately  reflect  ®Vpp.j.  for  those 
devices  for  which  the  conventional  interpretation  of  N  is  not  valid. 

For  such  devices  a  computer  solution  may  be  required. 
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E.  DEVICE  CHARACTERIZATION  FOR  COMPUTER  ANALYSIS  OF  LARGE  SEMICONDUCTOR 
CIRCUITS  (F.  A.  Lindholm) 


I.  Introduction 

trom  a  general  viewpoint,  the  starting  point  in  complex  microcircuit 
design  is  the  processing  (diffusion  of  impurities,  mask  generation,  -_tc.) 
used  in  fabrication.  The  end  point  is  a  circuit  that  meets  specifications 
in  system  performance  and  cost.  In  between  are  four  modeling  problems, 
each  involving  abundant  use  of  computer  calculation  as  the  flow  chart  of 
Fig.  1  illustrates. 

The  first  of  these  problems  seeks  to  describe  the  relationship  between 
processing  and  resulting  physical  make-up  (spatial  distribution  of  impurities 
and  defects,  corresponding  distribution  of  the  pertinent  material  parameters, 
etc.)  of  each  device  in  the  circuit.  In  the  second  modeling  problem  one 
solves  the  relevant  partial  differential  equations,  often  by  numerical  methods, 
to  relate  the  physical  make-up  to  the  spatial  and  temporal  dependence  of  such 
internal  variables  as  free  carrier  densities,  quasi-Fermi  levels,  charge 
density,  and  current  density.  Response  is  sought  for  excitation  pertinent 
to  the  proposed  system  application.  The  understanding  thus  gained  is  useful 
in  designing  the  devices  comprising  the  circuit.  In  the  third  modeling 
problem,  one  tries  to  link  the  physical  make-up  and  the  current  voltage 
behavior  of  each  device  at  its  terminals.  The  vehicle  for  describing  this 
behavior  is  usually  called  a  circuit  model  or  an  equivalent  circuit  The 
availability  of  adequate  circuit  models  is  essential  in  the  design  of  any 
circuit . 

Past  research  in  device  characterization  has  focussed  nearly  exclusively 
on  these  three  modeling  problems.  This  paper  deals  with  the  fourth  problem 
indicated  in  Fig.  1,  which  is  concerned  with  selection  of  models. 

The  maximum  size  of  semiconductor  circuit  amenable  to  computer  analysis 
and  the  efficiency  of  the  analysis  are  both  determined  substantially  by 


V 


51 


the  complexity  of  model  selected  to  present  each  device  in  the  circuit. 

The  choice  of  model  for  active  devices  such  as  the  bipolar  transistor  is 
particularly  crucial.  Since  the  invention  of  the  bipolar  transistor  in  1949, 
many  different  circuit  models  have  been  proposed  to  account  for  the  various 
phenomena  that  combine  to  determine  transistor  performance.  But  the  related 
problem  of  systematically  selecting  the  most  appropriate  model  to  represent 
each  transistor  in  a  given  circuit  environment  has  as  yet  received  scant 
attention.  The  present  paper  describes  a  method  for  systematic  model 
selection.  For  dc  circuits  operating  in  ordinary  environments,  it  has  proved 
capable  of  substantial  savings  in  cost  of  analysis,  and  it  potentially 
enables  simulation  of  larger  circuits  than  hitherto  possible.  *  Application 
this  method  to  the  simulation  of  irradiated  circuits  requires  extensions, 
some  of  which  are  suggested. 

II.  Reservoir  of  Available  Models 

To  provide  perspective,  we  review  certain  of  the  transistor  models 

available  for  selection.  First,  there  exist  the  classical  models:  the  Ebers- 

3 

Moll  model  and  its  equivalents.  In  this  paper  the  term,  Ebers-Moll  model, 
will  mean  the  model  in  the  original  sense , 

(1  +  b;>  rs  -  *s 

'  rS  (1  + 

■m 

h  +  XB  +  XC  ’  °’ 

in  which  the  slope  factors  associated  with  the  exponential  terms  are  unity 
and  in  which  the  parameters  of  the  square  matrix  are  constants,  independent 
of  the  current  or  voltages  at  the  device  terminals.  In  Eq.  (1),  kT/q  denotes 
the  thermal  voltage  and  8p  and  8^  denote  the  forward  and  reverse  current  gains, 
respectively.  This  form  of  the  Ebers-Moll  model  is  termed  its  transport 
version,  which  exhibits  certain  advantages  in  representing  modern  transistors 
relative  to  its  so-called  injection  version.^ 

The  Ebers-Moll  model  is  the  simplest  of  the  available  dc  transistor 
models,  and  this  is  its  main  attribute.  Risk  lies  in  its  use,  however,  for 
it  may  not  embody  enough  accuracy  to  yield  a  realistic  simulation.  To  avoid 
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such  risk,  workers  who  developed  package  computer  analysis  programs^  ^  derived 
more  accurate  models  by  retaining  the  form  of  the  Ebers-Moll  model  but  letting 
the  parameter  values  depend  on  current  and  voltage.  The  dependence  was 
determined  by  measurements  at  the  device  terminals.  Though  this  curve-fitting 
procedure  gave  greater  accuracy,  it  did  so  at  the  cost  of  increased  complexity 
and  its  consequent  disadvantages:  increased  memory  requirement  and  cpu  time, 
and  severe  limitations  of  the  size  of  circuit  amenable  to  simulation. 

Recently,  Gummel  and  Poon  proposed  a  model  that  retains  much  of  the  accuracy 
of  the  curve-fitting  approach,  yet  is  mathematically  more  compact.9,10 
Moreover,  the  Gummel-Poon  formulation  contains  a  set  of  models  that  facilitates 
trading  between  accuracy  and  simplicity.  This  is  a  highly  desirable  attribute.11 

III .  Model  Selection 

The  foregoing  models  constitute  a  basic  part  of  the  reservoir  of  available 
models  from  which  theanalyst  can  choose.  Our  objective  is  to  choose  from 
this  reservoir  for  each  transistor  in  a  circuit  under  study  the  simplest  model 
that  will  give  adequate  accuracy.  Thereby  we  avoid  the  severe  limitation  on 
circuit  size  and  the  excessive  cpu  time  and  memory  requirements  associated 
wit  i  representing  each  transistor  by  the  most  complex  model  available  as  well 
as  the  risk  and  the  prohibitive  cost  of  experimental  cut  and  try  associated 
with  representing  each  transistor  by  the  simplest  model  available.  The 
following  method  accomplishes  this  objective: 

1.  Read  in  the  network  configuration  using  the  Ebers-Moll  model  (with  para¬ 
sitic  resistances)  to  represent  each  transistor. 

2.  Analyze  the  network  using  an  analysis  program  of  choice. 

3.  Determine  for  each  transistor  whether  the  computer  currents  and  voltages 
at  its  terminals  are  consistent  with  the  approximations  underlying  the 
model.  This  self-consistency  check  is  accomplished  by  comparing  the 
computed  terminal  currents  and  voltages  with  onset  parameters,  defined 
below,  which  describe  the  onset  of  failure  of  the  approximations. 

4.  Increase  model  complexity  where  failure  occurs  in  accord  with  which  of 
the  onset  parameters  is  violated. 

5.  Analyze  the  network  again. 

6.  Continue  this  procedure  until  self-consistency  prevails  for  each  transistor 
in  the  network. 

IV .  Model  Failure  and  Onset  Parameters 

Implementation  of  this  method  requires  specification  of  the  onset  param¬ 
eters,  which  characterize  the  onset  of  model  failure.  This  we  now  do  for 
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the  Ebers-Moll  model.  For  simplicity  this  discussion  will  limit  consideration 
to  transistors  operating  in  the  static,  forward-active  model;  moreover,  we 
ignore  for  now  the  effects  of  a  radiation  environment.  For  these  conditions 
the  Ebers-Moll  model  fails  because  of  six  phenomena  neglected  in  its 
formulation: 

(1)  High  injection  in  the  base^'1"* 

(2)  Emitter  crowding1^ 

(3)  High  current  modes  in  the  collector1**-2^* 

(4)  Breakdown^ 2® 

(5)  Dominance  of  net  recombination  in  junction  transition 

29  10 

regions  and  at  surfaces  ' 

(6)  Early  effect  (base-width  modulation)’*1 

The  Early  effect  can  be  easily  incorporated  in  the  Ebers-Moll  formulation 
by  the  addition  of  a  single  parameter.  This  parameter,  termed  the  Early 
voltage,  is  defined  by  the  intersection  of  the  extrapolations  of  the  forward- 
active  Ic  -VCE  characteristics.  The  Early  effect  differs  from  the  other 
phenomena  listed  above  in  that  it  occurs  for  all  currents  and  voltages  in  the 
forward-active  region  of  transistor  operation.  In  contrast,  the  other  five 
phenomena  occur  only  over  certain  ranges  of  currents  and  voltages; 
consequently  the  onset  of  these  phenomena  can  be  described  by  specifying 
pertinent  critical  voltages  or  currents,  which  are  the  onset  parameters 
referred  to  earlier.  For  example,  the  usual  measurements^  employed  to 
characterizepunch  through,  avalanche  or  Zener  breakdown  will  determine  the 
critical  voltage  that  signals  the  onset  of  breakdown.  Further,  the  critical 
current  below  which  dominance  occurs  of  net  recombination  in  the  junction 
transition  region  or  at  surfaces  can  be  easily  found  from  the  slopes 
appearing  on  the  log  IB-VBE  characteristics . ^  The  onset  parameters  for 
the  other  three  phenomena  will  now  be  described. 

A.  High  Injection  in  the  Base 

High  injection  in  the  base  sets  in  when  the  charge  injected  in  the  base 

reaches  some  sizeable  fraction  F^g  of  the  net  ion  charge  present  there  in 

thermal  equilibrium.  Expressing  this  statement  by  charge-control  theory 

gives  the  onset  parameter  I  UTD  as 

C-HIB 

IC-HIB  "  FHIB  q  AE  NB/Tf 
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in  which  A_  denotes  the  eoitter  area,  N_  is  Gunnel's  number,  and  t.  Is 

®  C  f 

the  forward  charge-control  time  constant.  The  constant  FUI_  is  chosen  at 

HIo 

the  analyst's  discretion  in  accord  with  the  maximum  error  tolerable  in 

2 

assessing  model  adequacy;  typically  one  selects  values  of  1,  0.5,  or  0.25. 
If  the  collector  current  exceeds  the  phenomena  normally  associated 

with  high  Injection  set  in:  degradation  of  emitter  efficiency  and  current 
gain;  diminishing  of  base  resistance;  current  dependence  of  internal 
variables  such  as  lifetime  and  electric  field;  and  alteration  of  the 
Boltzmann  relation  between  barrier  voltage  and  injected  carrier  density. 

To  Illustrate  the  procedure  of  model  updating,  suppose  that  the 
Initial  computer  analysis  shows  for  a  particular  transistor  model  that 
the  collector  current  exceeds  but  that  no  other  oniet  parameter 

is  violated.  We  then  replace  the  errant  model  with  the  Gumrael-Poon  model, 
but  abbreviated  to  contain  only  those  parameters  necessary  to  account  for 
high  injection  in  the  base.  Other  models  are  similarly  updated  where 
necessary,  the  network  is  then  again  analyzed,  and  the  method  outlined 
earlier  is  followed  until  self-consistency  everywhere  prevails. 

B.  Emitter  Crowding 

14  J5 

According  to  present  theory,  '  emitter  crowding  begins  to  set  in 
when  the  collector  current  reaches  approximately 


1  .  y  . 

C-CR  CR  Rg/»8f 


(3) 


in  which  F^  again  denotes  a  constant  chosen  as  for  high  injection  at  the 
analyst's  discretion,  kT/q  is  the  thermal  voltage,  R_  is  the  transverse 

D 

base  resistance,  m  is  a  constant  depending  on  the  stripe  geometry,  and  Sj. 

Is  the  forward  current  gain.  If  the  collector  current  exceeds  I  ,  the 

C-CR 

effects  associated  with  emitter  crowding  start  to  become  severe:  the 
thresholds  of  high  injection  in  the  base  and  high  current  modes  in  the 
collector  arc  hastened,  with  accompanying  degradation  of  Beta  and  f^.,  and 
the  base  resistance  diminishes. 

To  update  a  model  to  represent  emitter  crowding,  one  can  use  a  version 
of  the  Gumtnel-Poon  model  abbreviated  to  contain  only  enough  complexity 
to  account  for  crowding.  The  resulting  model  is  moderately  complicated, 
however,  because  the  Gurarael-Poon  formulation  is  based  on  the  assumption  of 
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one-dimensional  flow.  It  accomplishes  a  first-order  modeling  of  the  effects 
of  multi-dimensional  flow,  such  as  those  accompanying  crowding,  by  curve¬ 
fitting  one  of  the  so-called  "push-out"  parameters.  Thus  computational 
advantage  and  improved  accuracy  both  may  derive  from  use  instead  of  a  model 
patterned  after  the  work  of  Chosh.  At  the  present,  insufficient  empirical 
and  computational  evidence  exists  to  suggest  which  method  of  updating  for 
crowding  is  superior. 

C.  HlRh  Current  Modes  in  the  Collector 

Two  basically  different  mechariisms  have  been  put  forward  to  explain 

the  degradation  of  Beta  and  f^  that  occurs  in  epitaxial  transistors  because 

of  high  currents  in  the  collector  region.  In  the  first  mode  of  operation, 

one  assumes  that  increasing  current  gives  rise  to  qussi-saturation  following 

16  20-23 

collapse  of  the  collector  space-charge  region.  *  In  the  second,  one 

assumes  Chat  increasing  current  leads  to  conditions  under  which  spacc- 

charge-limited  flow  is  avoided  by  the  occurrence  of  cither  base  widening  or 

17  18 

lateral  spreading  or  both.  * 

The  origin  of  these  two  different  modes  lies  in  the  two  opposing  tendencies 

indicated  in  Fig.  2  that  accompany  high  currents  in  the  collector  region. 

On  the  one  hand,  increasing  collector  current  reduces  the  effective  voltage 

applied  across  the  collector  space-charge  region;  thereby  this  region  tends 

to  collapse  as  current  rises,  leading  ultimately  to  quasi-saturation.  On 

the  other  hand,  increasing  collector  current  implies  the  existence  of  mobile 

electrons  whose  charge  tends  to  cancel  that  of  the  donor  ions;  thus  the 

collector  space-charge  region  tends  to  widen  with  increasing  current, 

producing  eventually  the  mode  associated  with  space-charge-limi ted  flow. 

24 

Which  mode  predominates  is  characterized  by  the  onset  parameter 


V 


crit 


<*) 


in  which  v 

sat 

the  collector 
the  thickness 


denotes  the  scatter-limited  velocity  of  mobile  electrons  in 

region,  u  denotes  the  low-field  mobility  there,  and  V„  denotes 
no  t 

of  the  low-conductivity  epitaxial  layer.  For 


Vcrit  4 
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Fig.  2  Depending  on  the  value  of  the  collector  voltage  relative 
to  the  value  of  the  onset  parameter  Vcrit,  the  collector 
space-charge  region  will  either  tend  to  collapse  of  to 
widen,  as  indicated  by  the  arrows.  Collapse  leads  toward 
the  node  of  quasi-saturation,  and  widening  leads  toward 
the  mode  associated. with  space-charge-lioited  flow. 
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quasi-saturation  occurs;  otherwise  space-charge-limited  flow  tries  to  start. 
When  the  collector  current  approaches  the  onset  parameters 

IC-HCC  "  FHCC^q  NDC  AE  Mno/WC^ 

either  quasi-saturation  or  the  mode  associated  with  space-charge-limited 
2  24 

flow  will  begin.  *  The  mode  that  prevails  is  found  by  the  test  involving 

V  .  Here  again  F..__  is  an  adjustable  constant  chosen  subject  to  the 
cn  t  hcc  2 

maximum  error  tolerable  in  assessing  model  adequacy.  The  first  term  in 
brackets  in  Eq.  (5)  denotes  the  conductance  of  the  high-resistivity  portion 
of  the  collector,  taken  apart  from  the  rest  of  the  transistor.  The  second 
term  in  brackets  accounts  for  the  field  dependence  of  mobility.  By  V_Dfr 
we  mean  the  total  voltage  (applied  plus  contact  potential)  across  the 
collector. 

In  updating  models  to  account  for  the  presence  of  high  current  modes 

in  the  collector,  one  should  carefully  note  which  mode  prevails.  If  the 

collector  current  exceeds  I_  and  the  test  with  V  . .  implies  the 

C-HCC  crit  v 

presence  of  the  mode  associated  with  space-charge-limited  flow,  then  base 

resistance,  beta  and  f  all  diminish  with  increasing  current.  Evidence 

17  1 

presently  available  suggests  that  the  functional  dependences  describing 
this  diminution  result  from  a  combination  both  of  base  widening  and  lateral 
spreading.  This  mode  therefore  requires  that  the  model  used  in  updating 
account  for  the  effects  of  multi-dimensional  flow.  As  was  discussed  earlier 
in  connection  with  emitter  crowding,  the  presence  of  multi-dimensional  flow 
raises  questions  concerning  the  adequacy  of  the  Gummel-Poon  formulation 
which  the  present  dearth  of  experimental  and  computational  experience  leaves 
unsettled.  If,  on  the  other  hand,  the  self-consistency  test  implies  the 
existence  of  the  quasi-saturation  mode,  the  associated  one-dimensional  flow 
suggests  that  one  can  update  with  confidence  using  an  appropriately 
abbreviated  version  of  the  Gummel-Poon  model. 

V.  Epilogue 

A  reservoir  of  transistor  models  exists  from  which  one  can  select. 

The  purpose  of  this  paper  has  been  to  describe  how  systematic  selection  can 
be  incorporated  into  the  computer  simulation  of  circuit  behavior  in  such 
a  way  that  the  models  selected  best  fit  the  circuit  environment  in  which 
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each  transistor  resides.  The  essence  of  the  method  involves  a  self-consistency 
test  for  the  assessment  of  model  adequacy.  In  this  test  an  initial  simulation 
yields  currents  and  voltages  which  are  compared  with  onset  parameters. 

Violation  of  any  onset  parameter  calls  for  an  updating  of  the  pertinent  model. 
The  model  used  as  the  replacement  in  the  next  simulation  is  not  the  most 
complex  model  available,  but  rather  a  model  only  complicated  enough  to 
represent  the  phenomena  implied  by  the  violated  onset  parameters.  Thus  Lhe 
method  puts  complexity  only  where  needed  and  only  as  much  as  needed.  In  this 
way  it  leads  to  the  selection,  for  each  transistor  in  a  circuit  being  simulated, 
of  the  model  of  least  complexity  that  will  give  adequate  accuracy.  Benefits 
follow:  basically  the  method  simplifies  computation,  saves  computational  cost, 
and  potentially  enables  the  simulation  of  larger  circuits  than  hitherto 
possible . 

It  is  sensible  that  any  modeling  techniques  be  conceived  and  tried  for 
conventional  circuits  before  subjecting  them  to  the  additional  difficulties 
imposed  by  a  radiation  environment.  The  applicability  of  the  method  described 
here  to  the  simulation  of  irradiated  circuit  behavior  has  yet  to  be  demonstrated. 
But  some  comment  about  potential  applicability  can  be  made. 

The  method  of  model  selection  applies  directly  to  the  dc  simulation  of 
transistor  circuits  suffering  permanent  radiation  damage.  The  likely  dominance 
of  net  recombination  current  from  the  transition  regions  can  be  included  by  use 
of  the  Gummel-Poon  formulation,  appropriately  abbreviated.  To  describe  the 
transient  behavior  of  a  circuit  exposed  to  time-varying  radiation,  however, 
one  must  generalize  both  the  reservoir  of  available  models  and  the  self- 
consistency  test  for  model  adequacy.  One  needs  a  hierarchy  of  models  of 
uifterent  accuracies  and  complexities  that  embodies  the  displacement  currents 
that  flow  and  the  photocurrents  generated  in  the  intrinsic  and  extrinsic 
portions  of  the  transistor.  The  test  for  model  adequacy  must  be  generalized 
to  assess  whether  the  transient  components  of  the  computed  currents  and 
voltages  are  consistent  with  the  approximations  underlying  each  model  used 
in  a  simulation. 
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FIGURE  CAPTIONS 


Fig.  1  The  four  modeling  problems  involved  in  modern  circuit  design. 

Fig.  2  Depending  on  the  value  of  the  collector  voltage  relative  to  the 
value  of  the  onset  parameter  ^cr^t»  the  collector  space-charge 
region  will  either  tend  to  collapse  or  to  widen,  as  indicated  by 
the  arrows.  Collapse  leads  toward  the  mode  of  quasi-saturation, 
and  widening  leads  toward  the  mode  associated  with  space-charge- 
limited  flow. 
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Radiation  Sensitivity  of  Amorphous  Semiconductors  (L.L.  Hench) 


STRUCTURAL  EFFECTS  ON  FAST  NEUTRON  RADIATION 
SENSITIVITY  OF  SEMICONDUCTING  GLASSES 
(H.  F.  Schaake,  A.  E.  Clark  and  L.  L.  Hench) 


Introduction 

Because  of  the  high  degree  of  structural  disorder  in  a 
glass,  amorphous  semiconductors  are  materials  whose  electrical 
properties  snould  be  relatively  insensitive  to  high  energy  radi¬ 
ation.  The  fact  that  a  certain  amount  of  structural  organiza¬ 
tion  is  present,  however,  suggests  that  a  suitably  high  radia¬ 
tion  dosage  may  affect  the  electrical  conductivity. 

The  structural  organization  in  a  glass  may  be  considered 
from  both  the  atomic  level  and  from  the  microscopic  level.  In 
the  first  class,  recent  pair  function  distribution  analyses  of 
amorphous  x-ray  diffraction  studies  have  shown  a  high  degree  of 
structural  correlations  as  far  as  the  seventh-neighbors^^  in 
glasses.^  ,2^  Geometrical  argumentsP’^  coupled  with  the  re¬ 
quirement  for  dense  packing, ha”e  suggested  that  there  is  even 
more  organization  in  the  glass  than  is  indicated  by  interpreting 
the  pair  function  distributions  as  random  arrays  of  bond  lengths 
and  angles.  That  is,  the  above  argu rents  suggest  that  the  exis¬ 
tence  of  a  p a r t i c u la r  -  a t omi c  coordination  at  one  point  in  the 
glass  results  in  a  very  high  probability  that  in  nearby  regions 
the  coordination  is  predictable,  i.e.,  often  the  same. 
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Conventional  heterogeneities  constitute  the  second  class 
of  organization  in  glasses.  Heterogeneities  may  result  from 

1)  microcrystallization  during  quenching  of  bulk  or  thin  films; 

2)  glass-glass  phase  separations;  3)  nucleation  and  growth  of 
crystals  from  a  glassy  phase;  and  4)  compositional  fluctuations. 
Such  features  are  extraordinarily  difficult  to  avoid  in-toto 
during  the  fabrication  annealing,  and  measuring  of  most  amor¬ 
phous  semiconductors. ^  ^ 

It  is  proposed  ♦'hat  when  a  high  energy  radiation  dosage 
is  sufficiently  large  to  alter  either  the  atomic  structural 
organization  in  the  glass  or  the  mi c ros t r uc t ur al  features, 
changes  in  the  electrical  properties  will  be  observed.  Such 
changes  have  been  observed  in  heterogeneous  semiconducting  glass¬ 
es  previously  rep o r ted .  ^  ^ ^ ^ 

In  the  present  study,  the  effect  of  fast  neutron  irradia¬ 
tion  on  the  electrical  properties  of  33  moleZ  K^PO^-b?  moleX 
V2°5  seraiconducting  glasses  has  been  investigated.  The  degree 
of  structural  organization  has  been  varied  by  subjecting  the 
glasses  to  different  thermal  treatments  before  irradiation.  A 
correlation  between  various  structural  features  and  fast  neutron 
damage  thresholds  is  considered  and  irradiation  dependent  chan¬ 
ges  in  electrical  conduction  mechanisms  are  proposed. 

Measurements  and  Analysis 

The  conductivity  as  a  function  of  frequency  of  a  series 
of  samples  was  measured  before  and  after  various  fast  neutron 
(70.1 MeV)  radiation  doses.  Details  of  sample  preparation,  con¬ 
ductivity  bridge  measurements  and  irradiation  techniques  are 
given  elsewhere. ^  ^  The  67%  K^PO^  samples  were 

melted  in  Pt  crucibles  in  air  and  rapidly  quenched  in  steel 
molds.  The  samples  analyzed  in  this  study  were  in  the  as- 
quenched  condition,  heat  treated  one-half  hour  at  288°C,  and 
heat  treated  30.8  hours  at  288°C.  As  previously  reported,  these 
thermal  treatments  produce  a  glass-in-glass  phase  separation  and 
microcrystallization,  respectively. ^  ^ 
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There  is  evidence  of  the  formation  of  some  sort  of  bar¬ 


rier  at  the  electrode-sample  interface  giving  rise  to  a  low 

/  o ) 

frequency  relaxation  mechanism.  It  is  impossible  to  distin¬ 

guish  between  interfacial  and  heterogeneous  relaxation  mecha- 

(9) 

nisms  on  the  basis  of  overall  conductivity  measurements  alone. 
Furthermore,  if  both  interfacial  and  heterogeneous  relaxation 
mechanisms  are  present,  then  the  interpretation  of  the  low  fre¬ 
quency  conductivity  will  be  complicated  by  the  intermixing  of 
these  effects.  Complicating  the  situation  further  is  the  fact 
that  the  conductivities  and  dielectric  constants  of  the  con¬ 
stituent  phases  are  frequency  dependent,  leading  to  an  even  more 
complex  overall  behavior. 

A  degree  of  simplification  can  be  achieved  if  certain  con¬ 
ditions  are  met.  The  conductiv  .  ty  for  both  heterogeneous  and 

(9) 

electrode  barrier  mechanisms  will  be  given  by  an  equation: 

<cs-cJw2rR 

o(u)  -  o'  +  - — r —  (1) 

1+“  TR 

If  to  >>  1/t  ,  then  Eq  .  (1)  becomes 


o  (u>) 


o' 


o 

o 


(2) 


i.e.,  a  plateau  is  reached  in  the  frequency  vs.  conductivity 
curves  (assuming  that  o',  cg,  e#  and  tr  are  independent  of  fre¬ 
quency  in  this  range).  If  only  electrode  barrier  effects  are 
present,  oq  is  the  D.C.  conductivity  of  the  sample  in  the  ab¬ 
sence  of  the  barrier.  If  there  are  heterogeneous  effects,  then 

o  is  as  given  in  Eq .  (2)  with  the  parameters  as  defined  in 

o 

Reference  9.  A  plateau  was  always  observed  in  the  present  mea¬ 
surements  (see  Fig.  1). 

Changes  in  the  conductivity  at  still  higher  frequencies 
will  be  due  to  changes  in  the  conductivities  and  dielectric 
constants  of  the  phases  present.  The  A.C.  conductivity  defined 
by 


°AC<“) 


ou  (u>) 
Meas 


(3) 
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therefore  is  due  to  frequency  dependent  mechanisms  operating 
within  the  phases. 

The  analysis  of  the  radiation  dependent  conductivity 
changes  will  center  largely  around  0^(01)  and  oq.  In  all  cases 
the  conductivity  could  be  represented  by  an  equation  of  the 
form 


o(u)  -  A  exp (-E^/kT)  +  Bu"  expC-Ej/kT) 

+  Cu2  exp(-E2/kT)  (4) 

The  first  term  is  o  ;  the  second  two  terms  comprise  o,_(u). 

o  AC 

In  the  present  as  well  as  in  previous  investigations  n  had  a 
value  from  0.8  to  1.3. 

Experimental  Results 

In  the  quenched  glass  (Fig.  2)  a  dose  of  1  x  1017  nvt  pro- 

0  8  2 

duces  no  change  in  the  u>  *  conductivity,  while  the  u  contri¬ 
bution  decreases.  Further  Irradiation  to  a  dose  of  2.1  x  1017 

2 

nvt  results  in  the  complete  destruction  of  the  u>  conductivity 

0 . 8 

in  the  observed  frequency  range  and  an  Increase  in  the  u 

conductivity  by  a  factor  of  -3. 

In  the  glass  heat  treated  for  one-half  hour  at  288°C,  a 

dose  of  5  x  10*^  nvt  leaves  the  A.C.  conductivity  unchanged 

0 . 8 

(Figs.  3  and  4).  The  activation  energy  for  the  u  ’  conduc- 

2 

tivity  remains  at  0.1  cv,  while  that  for  the  u  conductivity 

remains  at  0.0S  ev.  o  ,on  the  other  hand,  Increases  by  5Z, 

o 

while  the  activation  energy  decreases  slightly  from  -0.21  ev 
to  -0.19  ev  (Fig  5).  The  results  of  higher  dosages  on  a  simi¬ 
lar  glass  is  shown  in  Fig.  6,  After  a  dosage  of  1  x  10*7  nvt, 

0  8  2 
the  u  *  conductivity  Increases  by  a  factor  of  3,  while  the  u 

conductivity  increases  by  a  factor  of  almost  13.  Further  ir¬ 
radiation  to  a  dose  of  2.7  x  10  nvt  results  in  the  con- 

2 

ductlvity  decreasing  to  1.5  times  the  original  and  the  u  con¬ 
ductivity  decreasing  to  10  rimes  the  original.  At  the  same 
time,  o  behaves  as  follows: 
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Pre-irradiation  1.43  x  10~4 

17  , 

1  >  10  nvt  1.65  x  10 

2.7  x  1017  nvt  0.29  x  10~4 

Th“  <fata  for  the  glass  heat  treated  30.8  hours  Is  shown 

in  Fig.  7.  Tlieie  is  no  detectable  w°  * 8  conductivity  in  any 

samples.  A  dose  of  1  x  10  nvt  causes  the  w2  conductivity  to 

increase  by  a  factor  of  2.  Further  radiation  to  a  level  of 

2.1  x  10  results  in  no  change  in  the  u2  conductivity,  but 

docs  introduce  a  new  low  frequency  dependence  proportional  to 
1.3.  .  ( 7\ 

"  *  As  recently  discussed.'  '  the  behavior  of  o  1*  ns  fol¬ 

lows  : 

Pre-irradiation  5.2  x  10-3 

1  X  10  nvt  6.8  x  10“3 

2.1  x  1017  nvt  5.85  x  10"3 

Discussion 

The  data  presented  suggests  that  the  detectable  threshold 
for  neutron  irradiation  alteration  of  electrical  properties  in 
the  33  KPO^-67  V205  glass  is  5  x  103f  nvt.  In  the  glass  heat 
treated  for  one-half  hour,  this  dosage  resulted  in  no  change 
in  the  A.C.  conductivity,  and  a  «mall  increase  in  o  .  In  the 
quenched  glass  and  in  the  glass  heat  treated  foi  30.8  hours, 
this  threshold  is  under  1  x  1037  nvt;  examination  of  the  data 
suggests  that  it,  too,  is  about  5  x  1 0  3  **  nvt.  Since  the  30.8 
hour  glass  contains  -0.4  volume  fraction  crystals,  this  means 
that  over  a  range  of  from  '0  to  ‘0.4  volume  fraction  crystals, 
the  radiation  damage  threshold  is  5  x  10*^  nvt. 

Considering  the  controvery  which  apparently  exists  over 
the  source  of  the  u2  behavior  of  the  A.C.  conductivity  at  higher 
frequencies,'  the  present  data  is  very  interesting.  Of 
particular  interest  is  the  fact  that  the  w0,8  ',.id  u>2  conductiv¬ 
ity  possess  different  activation  energies,  both  of  which  arc 
smaller  than  the  oq  activation  energy.  It  is  also  seen  that 
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the  u0*6  and  u2  conductivities  change  by  different  amounts  with 
sufficiently  high  radiation,  indicating  that  they  are  in  fact 
due  to  different  mechanisms. 

In  this  research,  u  A.C.  conductivity  has  been  observed 
in  every  vanadate  glass  and  glass-ceramic  measured,  while  it 
has  never  been  observed  in  a  variety  of  ionically  conducting 
glasses  and  glass-ceramics  investigated  with  the  same  dielectric 
facility. 

An  «  A.C.  conductivity  dependence  is  not  predicted  by  a 

random  hopping  center  model. However,  as  discussed  in  the 

Introduction,  a  glass  is  not  a  random  array  of  atoms,  but  con- 

2 

tains  regions  ot  equivalent  ionic  coordination.  The  u  conduc¬ 
tivity  can  be  attributed  to  the  regions  of  structural  cquiva- 
lence.  An  u  dependence  will  be  observed  if  there  exists  in 
the  glass  a  considerably  larger  number  of  short  relaxation 
times  than  would  be  present  in  a  random  array  of  hopping  cen¬ 
ters  (the  hops  between  centers  randomly  distributed  gives  rise 
to  the  <j0'8  A.C.  conductivity).  Such  a  situation  oppears  to  be 
possible  in  a  real  glass  due  to  the  high  degree  of  short  r»ngc 
order. 

To  demonstrate  this,  we  shall  digress  briefly  to  consider 
the  case  of  amorphous  silicon.  This  digression  is  necessary 
because  of  the  lack  of  structural  information  on  the  vanadate 
glasses.  Pair  function  analyses  of  any  three  or  four  element 
glass  hn'j  not  been  reported.  However,  considerable  structural 
Information  has  been  amassed  for  amorphous  silicon.  Grigorovicl 
and  others  ^  have  argued  that  the  "mistake"  in  silicon  which 
leads  to  an  amorphous  rather  than  a  periodic  structure  is  the 
occurrence  of  adjacent  Sl^  tetrahedra  in  the  eclipsed  (cis-) 
rather  than  the  correct,  staggered  (trans-)  configuration,  i.e., 
a  rotation  of  60°.  This  is  an  error  in  third  nearest  neigh¬ 
bors  in  the  radial  distribution  function;  first  and  second 
neighbor  distances  remain  unchanged.  Looking  at  these  molecular 
units,  the  ground  state  energy  in  the  cis-Si  will  be  higher  than 
in  trans-Si.  The  c i s-con f i gurat ion  forces  the  wave- f unc t ion  to 
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have  sore  chan  che  trans-configuration  urvature;  hence  the 
kinetic  energy,  and  by  the  virlal  theorem,  che  total  energy  is 
higher. 

Grigorovici  has  also  argued  that  these  mistakes  do  not 
occur  randomly;  in  order  to  achieve  a  high  densitv  of  packing 
the  mistakes  occur  in  groups  --  that  is,  the  existence  of  a  cis- 
grouplng  about  one  Si  atom  occasions  the  existence  of  one  or 
more  c  i  s-groupi ngs  about  nearby  Si  atoms.  In  the  extreme  case 
of  grouping  of  these  structural  mistakes,  we  are  led  to  the 
Vornoi  dodecahedron.  Amorphous  silicon  may  therefore  be  de¬ 
scribed  as  a  mixture  of  clusters  of  mistakes  and  regions  of 
crystalline-coordinated  Si.  We  shall  now  make  the  following 
assumption:  That  th®  energy  of  the  first  excited  state  of  the 
cis-conf iguration  is  lower  than  the  first  excited  state  of  the 
t rans -c on f i gu ra t 1 on .  (The  validity  of  this  assumption  is  cur¬ 

rently  being  tested  by  direct  calculation;  it  certalnlv  seems 
reasonable.)  Thus,  when  we  investigate  the  wave funct ions  of 
che  conduction  states  in  the  glass,  those  in  the  conduction 
band  or  localized  Just  below  it,  we  will  find  that  those  having 
the  lowest  energy  have  large  contributions  from  the  cis-group- 
lng  "molecular"  wave f unc t i ons ,  and  small  contributions  from  the 
trans-grouping  "molecular"  wave f uncc ions .  Thus,  if  the  above 
assumption  is  correct,  the  localized  states  at  the  bottom  of 
the  conduction  band  arc  localized  on  clusters  of  mistakes. 

Two  types  of  hops  may  occur  for  a  corrier  located  in  the 
localized  states  at  the  bottom  of  the  conduction  band*:  it  may 

*NOTE:  When  the  overlap  potential  energy  (Jj.  ■  < i  !  r  ^  |  J  >  for 

a  simple  onc-electro.i  model)  between  two  states  in  the 
same  cluster  sufficiently  exceeds  the  e 1 ec t ron- 1  a 1 1 ice 
b  Indin;  energy  Ef  for  a  single  state  E*  »  l/2Kr^, 
where  n  is  the  dielectric  constant  and  r  £  £  is  the 
radius  of  static  (rn  -  <i|r|i>),  then  the  electron 
will  be  localized  on  both  sites;  i.c.,  the  minimum 
energy  wavefunction  is  |i+J>.  However,  a  transition 
to  i-j>  «ay  still  take  place,  since  the  dipole  moment 
<i-j|x|i+j>  is  not  zero,  and  the  following  argument  is 
valid. 
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hop  to  another  cluster,  or  it  may  hop  to  another  site  in  the 
same  cluster. 

For  a  hop  through  a  given  distance,  the  former  process 
will  result  in  a  longer  relaxation  time  than  the  latter.  This 
is  because  the  overlap  between  localized  states  on  different 
clusters  is  considerably  smaller  than  the  overlap  between  two 
states  on  the  same  cluster.  Therefore,  as  the  hopping  dis¬ 
tance,  d,  decreases,  the  average  r  laxation  time,  r(d),  de¬ 
creases,  first  in  agreement  with  the  random-center  distribution 
model,  then  more  rapidly,  since  the  probability  is  that  shorter 
hops  occur  within  a  mistake  cluster.  This  results  in  a  larger 
number  of  hops  with  short  relaxation  times,  and  a  smaller  num¬ 
ber  of  hops  with  moderate  relaxation  times,  and  hence  to  u2 
ant^  w  conductivity.  The  u  ®  conductivity  may  be  attributed 
to  hops  between  different  clusters  of  mistakes,  while  u2  con¬ 
ductivity  results  from  hops  within  a  cluster. 


A  similar  situation  to  that  discussed  above  for  amorphous 
silicon  should  exist  in  the  vanadate  glasses.  The  basic  struc¬ 
tural  unit  in  V-0  glasses  is  known  to  be  a  vanadium  ion  surrounded 
by  six  oxygen  ions  in  octahedral  coordination.  From  the  macro¬ 
scopic  standpoint,  the  oxygen  ions  form  a  densely  packed  struc¬ 
ture  (in  the  crystalline  vanadium  oxides,  the  oxygen  ions  form 
a  close-packed  structure,  distorted  slightly  by  the  vanadium 
ions),  with  vanadium  ions  distributed  in  octahedral  interstices 
to  satisfy  charge  requirements.  The  difference  between  a  V4+ 
and  a  V  ion  is  therefore  best  distinguished  by  the  vanadium 
ion  radial  distribution  function.  In  the  glasses  under  study, 
the  existence  of  other  ions  (K,  P)  with  different  charges  and 
coordination  modifies  this  simple  interpretation. 

As  in  silicon,  it  seems  probable  that  if  a  mistake  occurs 
which  causei!  a  V  site  to  have  the  vanadium  ion  coordination 
of  a  V  +  site  [which  we  shall  designate  (V4)5+),  then  a  number 
of  surrounding  sites  will  have  the  same  error  leading  to  clus¬ 
ters.  Note  that  long-range  charge  neutrality  can  be  met  by 
having  the  surrounding  region  vanadium  deficient.  Because  of 
the  excess  positive  charge  in  the  vicinity  of  (V4)5+  sites,  the 
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localized  states  at  the  bottom  of  the  conduction  band  (assuming 

one  exists)  will  be  peaked  on  these  sites. 

As  in  silicon,  the  clustering  of  mistakes  would  lead  to 

0.8  ,  2 

u  and  w  A.C.  conductivity.  The  observed  difference  in  acti¬ 
vation  energy  between  the  two  regions  of  u  dependence  (Figs. 

3  and  * )  suggests  that  there  is  a  variation  in  the  average  ener¬ 
gy  levels  of  the  clusters.  In  other  words,  there  appears  to  be 
a  small  variation  in  the  energy  levels  of  a  single  cluster. 
However,  if  the  energy  levels  within  each  cluster  are  averaged, 
then  there  will  be  a  larger  variation  in  these  averages.  Thus, 
as  the  temperature  is  increased,  hops  become  possible  between 
morn  clusters,  yielding  a  higher  activation  energy  than  for 
hops  within  a  cluster.  This  variation  in  levels  between  clus¬ 
ters  could  conceivably  result  from  the  P^,  K20,  or  V*  +  im“ 
purity  ions  in  the  glass. 

The  foregoing  is  summarized  in  the  qualitative  electronic 
structure  illustrated  in  Fig.  9  for  the  glass  heat  treated  one- 
half  hour  at  288°C. 

A  radiation  dosage  of  5  x  10^  nvt  in  the  one-half  hour 
heat  treated  glass  results  in  only  a  small  change  in  the  o 

o 

activation  energy  (Figs.  3,  4  and  5).  This  corresponds  to  a 
decrease  in  the  energy  difference  between  the  Fermi  level  and 
the  band  edge.  The  effect  of  this  level  of  radiation  is  to  con¬ 
tinue  the  process  begun  by  the  heat  treatment.  There  is  also 

a  small  Increase  in  the  magnitude  of  o  .  This  has  been  pre- 

60  ® 

vioualy  shown  through  Co  y-r experiments  to  be  due  to  yray 

ionization  of  carriers  through  Compton  scattering. 

Further  irradiation  to  a  dosage  of  1  x  1017  nvt  (Fig.  6) 

does  change  the  A.C.  conductivity,  increasing  both  the  u7  and 
0  B 

the  u  conductivity.  This  would  occur  if  the  glass  structure 
became  more  disordered,  leading  to  the  creation  of  more  clus¬ 
ters.  An  increase  in  the  number  of  donor  ions  participating 

also  evidently  occurs,  resulting  from  the  continuing  y-ray  ex¬ 
posure,  leading  to  the  higher  oq  and  part  of  the  increased  A.C. 
conductivity.  The  oq  activation  energy  may  also  change  when 
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the  dosage  Is  increased  to  2.1  x  IQ17  nvt,  the  clusters  of  mis¬ 
takes  evidently  become  more  disordered,  leading  to  a  further 
decrease  in  the  u  conductivity,  and  a  decrease  in  the  u®'® 

conductivity.  The  band  edge  energy  E  also  increases  due  to 

c 

tie  increased  disorder,  causing  more  localized  states  to  exist, 

resulting  in  a  decrease  in  o  . 

o 

In  the  quenched  glass  (Fig.  2),  the  disorder  is  evidently 
sufficiently  great  that  radiation  causes  the  glass  to  approach 
the  random  hopping  center  model;  that  is,  the  clusters  become 
sufficiently  disordered  that  the  number  of  hops  contributing 
to  the  u  '  conductivity  greatly  increases. 

The  30.8  hour  heat  treated  glass  presents  a  different  pic¬ 
ture  (Fig.  7).  If  the  clusters  of  mistakes  are  large  and  well 

separated,  then  the  A.C.  conductivity  will  be  proportional  to 
2 

u  as  observed.  These  clusters  may  be  large  enough  to  be 
recognized  os  a  separate  phase  as  indicated  by  x-ray  diffrac¬ 
tion  results.  Radiation  apparently  first  Increases  the  number 
of  clustered  regions.  Since  they  remain  well  separated,  this 
results  mainly  in  an  Increase  in  the  u  conductivity.  Further 
irradiation  then  produces  a  splitting  up  of  the  clusters,  re- 
suiting  in  the  appearance  of  u  ’  conductivity.  That  the  lower 
frequency  conductivity  is  proportional  to  u1 ’ 3  may  have  one  of 
two  explanations:  (1)  Since  the  glass  is  electrically  hetero¬ 
geneous,  if  one  phase  has  a  conductivity  proportional  to  u®  ‘  ® 

2 

and  the  other  to  u  ,  then  the  composite  conductivity  when  the 
conductivities  of  the  two  phases  are  about  equal  will  be  given 
by : 

°A  .  C  .  ^  “  Cu" 

where 

n  •  v  (0.8)  +  ( 1  — v  )2 
c  c 

with  vc  the  volume  fraction  of  the  precipitate  phase,  and  ns- 
suraing  the  precipitate  phase  has  u  conductivity.  With  v  • 

0.4, 
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n  "  1.5 

which  1 6  about  what  is  observed.  (2)  An  alternative  explanation 
Is  that  as  new  clusters  are  created,  they  a  re  not  randomly  dis¬ 
tributed  but  ore  clustered  together.  This  means  tl  at  shorter 
relaxation  times  are  more  numerous  than  in  the  random  cluster 
distribution  model,  resulting  in  a  frequency  coefficient  greater 
than  0.8,  but  less  than  2. 

Summary 

1.  Due  to  the  existence  of  structural  organization  in 
amorphous  semiconductors,  there  exists  a  threshold  radiation 
dosage  above  which  changes  in  electrical  properties  are  ob¬ 
served.  For  the  vanadate  glasses  in  the  system  33  KPO.j-67 

this  level  is  5  x  10^  nvt. 

2 

2.  The  u  conductivity  is  attributed  to  electron  hops 
within  clusters  of  structural  mistakes,  that  is,  within  regions 
of  the  glass  where  the  atomic  coordination  deviates  markedly 
from  the  crystalline  coordination. 

3.  The  observed  changes  in  the  A.C.  and  D.C.  conductivity 
can  be  attributed  to  changes  in  donor  concentration  as  a  result 
of  Compton  scattering  from  y-rays  and  from  changes  in  atomic 
structure  resulting  from  neut ron- induce d  damage. 
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Figure  Captions 


Fig.  1 

Fig  2 
Fig.  3 
Fig.  4 
Fig.  5 
Fig.  6 
Fig.  7 
Fig.  8 


-  Total  conductivity  versus  frequency  for  a  33  moleZ 
KPOj-67  raoleZ  V205  glass  heat  treated  30.8  hours. 
oq  Is  as  described  in  the  text. 

-  A.C.  conductivity  (o,\.C.)  of  quenched  glass  for  vari¬ 
ous  radiation  doses. 

-  A.C.  conductivity  (  0  a  .  C  .  ^  °  f  glass  heat  treated  0.5 
hour.  No  Irradiation. 

-  A.C.  conductivity  (oa.c.)  of  glass  heat  treated  0.5 
hour.  Dosage  5  x  1016  nvt. 

-  0o  versus  reciprocal  temperature  for  glass  heat 
treated  0.5  hour. 

-  A.C.  conductivity  ( o a . C . 1  of  glass  heat  treated  0.5 
hour  for  various  doses. 

-  A.C.  conductivity  (oa.C.)  of  glass  heat  treated  30.8 
hours  for  various  doses. 

-  Proposed  electronic  structure  for  glass  heat  treated 
0.5  hour.  Energies  Indicated  are  approximate,  nnd  In 
electron-volts.  Ep  is  the  Fermi  level:  Ec  the  criti¬ 
cal  energy  separating  band  and  localized  states. 
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B.  EFFECTS  OF  MICROSTRUCTURE  ON  THE 

RADIATION  STABILITY  OF  AMORPHOUS  SEMICONDUCTORS 
(L.  L,  Hench,  A.  E.  Clark  and  H.  F.  Schaake) 


Introduction 

Because  of  the  structural  randomness  of  amorphous  semi¬ 
conductors  they  offer  the  promise  of  radiation  insensitive 
electronic  properties.  Studies  have  demonstrated  that  there 
is  only  minor  alteration  of  the  electrical  properties  of  non¬ 
crystalline  V2°5*P2°5  and  K2°“V2°5~P2°5  amorPhous  semiconduc¬ 
tors  at  fast  neutron  fluences  of  4  x  10*2  nvt  *»nd  y-ray  dosages 
of  1.2  x  108  rads. (1,2) 

(2  3) 

However,  it  has  been  shown  for  both  oxide  *  and  chal- 
(4  s ) 

cogenide  '  amorphous  semiconductors  that  structural  hetero¬ 
geneities  are  present  in  many  compositional  systems.  The  data 
indicate  that  heterogeneities’  result  from:  1)  mf.c  rocrys  ta  1 1  i  - 
zation  during  quenching  of  bulk  or  thin  films;  2)  glass-glass 
phase  separation;  3)  nucleation  and  growth  o>f  crystals  from  a 
glassy  phase;  and  4)  compositional  fluctuations.  Can  such 
structural  heterogeneities  influence  the  radiation  sensitivity 
of  these  materials?  Several  recent  reports''  show  chalcog- 

enide  thin  film  switching  devices,  which  contain  filamentary 

16  2 

heterogeneities,  are  stable  up  to  fluences  of  10  n/cm  and 
3.3  x  102  R  of  fission  y-rays.  However,  heterogeneous  vanadia 
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based  amorphous  scniconduc tors  in  bulk  forn  exhibit  marked 
changes  in  electronic  behavior  at  fluences  of  1  x  10^  nvt  and 
above. <2-10> 

Consequently,  one  objective  of  the  present  paper  is  to 
extend  the  range  of  the  previous  work  by  reporting  the  extent 
of  radiation  sensitivity  of  a  33  mole  X  KPOj,  67  Bole  X 
glass  when  crystallized  to  a  large  volume  fraction,  -0.4. 
Secondly,  a  basis  for  selecting  compositions  for  Improved  ir¬ 
radiation  stability  is  suggested. 

Experimental  Procedure 

The  glasses  were  melted  in  air  in  Pt  and  melt  cast  at 

900°C  into  graphite  molds.  The  1.8  cm  diameter  discs  were 

heat  treated  at  288°C  for  30.8  hours  in  air  after  polishing. 

(2  11  12) 

Previous  studies  ’  ’  of  thermal  treatment  of  glasses  in 

the  K20_P205"V2°5  systcm  have  shown  that  a  90  X  phase  separa¬ 
tion  is  enhanced  by  a  thermal  treatment  for  1/2  hour  at  288°C. 
The  matrix  appears  to  order  during  heating  with  the  D.C.  con- 
ductivity  increasing  by  10  .  Longer  treatments  at  288°C  -esult 
in  the  precipitation  of  P205-2(V205>  crystallites  until  at 
30.8  hours  -0.4  volume  fraction  of  crystals  are  present.  D.C. 

conductivity  increases  during  this  process  to  a  2*»°C  value  of 
, -3  -1  -1 

5  x  10  ohm  cm  .  Evaporated  double  guard  ring  gold  elec¬ 
trodes  were  annealed  at  150°C  for  1/2  hour.  Measurements  were 
made  as  described  elsewhere . 

Irradiation  was  done  in  the  Wright-Patterson  Air  Force 

Base  test  reactor  having  a  fast  neutron  flux  capability  of 
13  2 

1.5  x  10  n/cm  (>0.1  :neV)  .  Cadmium  wrapping  protected  the 
gold  electrodes.  The  sample  temperature  was  maintained  near 
50°C  during  the  run.  Samples  held  at  50°C  without  radiation 
showed  no  property  changes. 

Results 

The  conductivity  frequency  dependence  at  25°C  of  the 
288°C/30.8  hour  sample  is  shown  in  Fig.  1  before  and  after 
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cumulative  exposures  of  1.0  x  1017  nvt  and  2.1  x  1017  nvt. 

The  total  conductivity  first  increases  with  nvt  and  then  de¬ 
creases.  The  conductivity  increase  is  similar  to  that  seen  in 
the  quenched  and  heat  treated  glasses  previously  studied . ^3 ,2^ 
The  decrease  in  ctot  at  higher  nvt  was  also  observed  in  the 
288°C/0.5  hour  glasses  at  nvt  >2  x  1037.^33^ 

The  A.C.  conductivity,  oA_(w)  •  a  (to)  -  a  ,  before  and 

al  meas  a 

After  irradiation  is  shown  to  be  a  function  of  u  in  Fig.  2. 

of  the  multiphase  material  is  considered  to  be  the  average 
of  the  A.C.  conductivities  of  the  phases  and  a  is  defined  as 

Aft® 

shown  in  Fig.  1.  There  is  no  detectable  u  *  conductivity  in 
any  samples.  A  flucnce  of  1  x  10  nvt  causes  the  u>  conduc¬ 
tivity  to  increase  by  a  factor  of  2.  Further  radiation  to 

17  2 

2.1  x  10  results  in  no  change  in  the  u  conductivity,  but 

does  introduce  a  new  low  frequency  dependence  of  w3'3.  The 

behavior  of  o  is  as  follows: 
o 


Pre-i rr  diat ion 
1  x  1017  nvt 
2.1  x  103  7  nvt 


oq  *  5.2  x  1033  ohm  3cm  3 
oq  “  6.8  x  10  3  ohm  3cm  3 
oq  ■  5.9  x  10  3  ohm  3cra  3 


The  loss  angle  of  the  material  is  greatly  increased  with 
1  x  1C  nvt,  as  shown  in  Fig.  3.  After  2.1  x  1037  a  large 
loss  peak  appears  in  the  t an  6- f req uency  plot. 


Discussion 

The  results  show  that  fast  neutron  irradiation  affects: 

3>  oDC  magnitude,  2)  oAC-u  dependence,  and  3)  tan<S  peaks.  The 
increases  in  may  be  attributed  to  an  increase  in  mobile 

carrier  concentration  due  to  y-ray  ionization  processes.  As 
shown  y-ray  dosages  of  1.2  x  10  rads  produce  a  15%  increase 
3n  °DC  3n  ^2^5~^2^5  Brasses  •  A  disordering  of  the  glass  ma¬ 
trix,  reversing  the  thermal  treatment  effects,  is  proposed  for 
the  conductivity  decrease  after  2.1  x  1037  nvt.  The  matrix 
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disordering  would  Isolate  high  conductivity  crystallites  pro¬ 
viding  barriers  suitable  for  Maxwe 1 1-Wagne r-Wi 1  la rs  (Ml.’S) 
peaks  which  appear  at  2.1  x  10*^  nvt. 

Explanations  for  the  variations  in  u  dependence  are  more 
speculative.  Localized  state  pairs  nay  be  present  in  the 

r 

heterogeneous  material  which  would  result  in  an  u  dependence 

(14) 

via  adiabatic  hopping.  The  irradiation  may  produce  a  suf¬ 

ficient  redistribution  of  nearest  neighbors,  as  reported  for 
fused  silica^*’^  that  a  complex  ^  behavior  produced.  It  is 
also  posfible  that  excess  carrier  generation  or  other  effects 
could  produce  an  electrode  barrier  behavior  altering  the 
dependence. 

The  nvt  threshold  for  degradation  of  electrical  properties 
related  to  heterogeneities  is  plotted  as  a  function  of  soften¬ 
ing  point*  for  heterogeneous  K.O-V  0,-P_0  and  Li_0-2Si0 
(1916)  l  l  J  i  1  l  i. 

glasses  '  in  Fig.  4.  In  both  systems  the  onset  of  changes 

in  the  MW 3  loss  peaks  is  taken  as  nvt  threshold.  The  threshold 
for  significant  structural  change  of  fused  silica  and  crystal¬ 
line  quartz  is  used  for  the  S  if^  data  point^*^  since  partially 
devitrified  SiO^  has  not  been  studied. 

Figure  4  shows  that  as  the  softening  point  increases  the 
threshold  for  damage  to  heterogeneities  increases.  This  is 
reasonable  since  the  softening  point  is  related  to  the  magni¬ 
tude  of  bond  energies  in  the  material.  Tiius  a  larger  cumula¬ 
tive  number  of  events  is  requires  to  damage  the  boundaries 
and  crystallites  of  higher  bond  energy  materials. 

Extrapolation  of  the  plot  in  Fig.  4  indicates  that  damage 
to  ordered  heterogeneities  in  chalcogenidc  systems,  e.g.,  on- 
state  filaments,  may  occur  at  >5  x  10^  nvt.  The  higher 
threshold  for  glasses  with  larger  crystallites  indicates 

that  the  damage  threshold  of  devices  may  be  increased  as  the 
on-state  filament  diameter  is  increased.  The  correlation 


*Although  other  criteria  might  be  used  for  comparison,  e.g., 
Tx ,  Tg»  ,  an  advantage  of  the  softening  point  is  that  it  is 
sensitive  to  previous  thermal  history  and  heterogeneities  in 
the  material. 
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shown  in  Fig.  4  suggests  that  another  way  to  improve  radiation 
tolerance  is  to  increase  the  softening  points  of  amorphous 
semiconductors  by  coraoosi t ional  modification. 
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Flp.u vc  Captions 


1  Effect  of  fast  neutron  irradiation  on  the  total  conduc¬ 
tivity  of  a  33  mole  Z  ^  PO^  -  67  mole  Z  ?2°5  glass  pre¬ 
viously  heat  created  for  30.8  hours  at  288°C. 

2  Change  in  A.C.  conductivity  of  a  heat  treated  vanadia- 
phosohate  glass  with  fast  neutron  irradiation. 

3  Influence  of  fast  neutron  irradiation  on  the  loss  tangent 
of  a  partially  crystallized  vanadla-phosphatc  glass. 

4  A  correlation  between  the  dllatometrlc  softening  point 
of  glasses  with  the  threshold  flucncc  for  fast  neutron 
structural  changes. 
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FREQUENCY  HZ 


33  Mol  %  (K2P03)-V205  G 
Treated  30.8  hrs  at  288*C 


FREQUENCY  (HZ) 


iOO 


Glass  Softening  Pomf  vs  Fast  Neutron  Damage 
Threshold  for  Structural  Change  (nvt) 


C.  ELECTRONIC  MATERIALS  INTERFACIAL 
CHARACTERIZATION  PROGRAM 
(G.  J.  Scott  and  L.  L.  Hench) 


Introduction 

Previous  research  reported  in  this  contract  program  demon¬ 
strated  that  the  dielectric  properties  of  electronic  substrates, 
radomes,  and  window  materials  are  strongly  influenced  by  rela¬ 
tive  humidity  and  temperature.  This  behavior  was  attributed 

to  penetration  of  water  vapor  along  grain  boundaries  and  micro¬ 
pores  forming  a  polar  condensate  which  contributes  to  inter¬ 
facial  and  dipolar  losses  in  the  dielectric  properties  of  the 
material . 

A  review  study  presented  in  the  last  progress  report  of 
this  program  discussed  a  variety  of  means  suitable  for  charac¬ 
terizing  the  surface  of  ceramic  powders.  Understanding  powder 
surfaces  is  an  especially  important  concern  of  electronic  cer¬ 
amics  processors  and  users  today.  This  is  because  it  lias  been 
shown  that  agglomeration  of  powders  controls  the  microstructural 
heterogeneity  of  fine-grained  sintered  ceramics  such  as  insula¬ 
tor  substrates,  ferrites,  ferroelectrics ,  etc.  Agglomeration 
of  powders  is  a  surface  phenomena  which  is  very  poorly  under¬ 
stood.  Consequently,  only  by  obtaining  a  more  fundamental 
understanding  of  powder  surfaces  will  it  be  possible  to  control 


a gg lome r a t i on  and  eliminate  heterogeneous  sintered  microstruc¬ 
tures  ar.d  thereby  improve  the  reliability  and  performance  of 
a  variety  of  thin  film  and  bulk  devices  requiring  ceramic  com¬ 
ponents  . 

It  was  shown  in  the  review  study  that  indirect  means  of 
powder  surface  characterization  are  the  only  practical  tools 
available  for  bulk  studies  at  the  present  time.  Two  of  the 
most  powerful  indirect  methods  for  surface  characterization 
are  differential  thermal  analysis  adsorbtion  isotherms  and  di¬ 
electric  adsorbtion  isotherms.  We  have  recently  completed  an 
experimental  facility  which  couples  these  two  techniques  to¬ 
gether  for  surface  characterization  studies.  One  of  the  con¬ 
sequences  of  the  combined  surface  analytical  approach  is  hav¬ 
ing  to  analyze  the  formidable  amount  of  data  collected  in  such 
a  study.  Data  are  accumulated  for  eleven  parameters  and  vari¬ 
ables.  Calculations  must  be  made  on  the  data  and  a  large 
number  of  cross-plots  constructed  in  order  to  interpret  the 
significance  of  the  results. 

In  order  to  make  such  an  undertaking  tractable,  a  com¬ 
puter  program  has  been  written,  debugged,  and  operated.  The 
following  is  a  description  of  that  program. 

Overview 

1.  The  "DIEL"  set  of  programs  gives  graphical  results  of 
complete  studies,  including  temperature  dependence,  of  dielec¬ 
tric  property  response,  on  a  single  sample.  The  five  programs 
are  designed  to  accommodate  data  collected 

a.  at  up  to  five  temperatures, 

b.  during  adsorptive  or  desorptive  modes  (and  here¬ 
after  mode  refers  to  the  ad-  or  desorptive  run  of  an  ex¬ 
periment  at  a  given  temperature), 

c.  up  to  20  relative  humidities  per  temperature  and 
mode  combination  (hereafter  T-MC) , 

d.  and  up  to  25  frequencies  per  relative  humidity 
and  T-MC. 
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2.  These  programs  can  plot  up  to  16,200  points  on  nine 
types  of  cuts  through  an  eleven-dimensional  space,  with  up  to 
24,200  specified  data  values,  on  up  to  270  sets  of  axes. 

3.  The  plots  are  of  nine  types: 

a.  l->3:  loss  tangent,  charging  permittivity,  and 

loss  permittivity  as  ordinates  versus  log  frequency  as 
the  abscissa,  parametrically  at  up  to  three  relative 
humidities,  for  a  given  T-MC.  The  usual  symbolic  names 
for  the  ordinate  variables  are  tan  del,  e',  and  e"  re¬ 
spectively  . 

b.  4->9  :  charging  and  loss  permittivity  and  capaci¬ 

tance  change  are  plotted  against  the  quantity  of  water 
adsorbed.  Symbolic  names  are  respectively  e',  e",  Ac, 

and  milligrams  per  gram.  The  parameter  for  these  graphs 
is  the  set  of  six  special  frequencies  --  data  for  type  4, 

6  and  8  graphs  is  taken  at  high  frequency,  viz.,  1  MHz, 

500  KHz,  100  KHz;  for  types  5,  7  and  9  at  lower  frequency, 
20  KHz,  5  KHz  and  500  Hz. 

4.  There  are  five  programs  in  the  set.  Their  names  and 
f unc  tions  are : 

a.  DIELPNCH  -  data  reduction;  punch  a  deck  for  read¬ 
ing  data  and  calculated  values  into  a  direct  access  stor¬ 
age  unit  . 

b.  DIELCELL  -  arrange  data  according  to  the  tempera¬ 
ture,  mode,  relative  humidity  and  frequency  breakdown, 
explained  above,  into  the  data  cell  (a  magnetic  chip-on- 
disc  memory  unit;  or  into  a  magnetic  tape. 

c.  REPAIR  -  change  several  isolated  data  values  in 
the  data  cell,  as  necessary. 

d.  D1  -  draw  the  graph  of  types  1  through  3. 

e.  D2  -  same  for  types  3  to  9. 

5.  Detailed  explanations  of  these  programs,  construction 
rules  for  their  data  decks,  and  listings  of  the  programs  follow. 
An  explanatipn  of  the  job  control  language  ("JCL")  cards  re¬ 
quired  for  access  to  the  programs'  software  is  an  appendix. 
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6.  Five  subprograms,  which  do  the  actual  work  of  putting 
the  graphs  on  paper,*  are  shared  by  D1  and  D2.  We  will  refer 
to  them  as  they  are  referred  to  during  operation  of  D1 ,  so  that 
they  will  physically  separate  Dl’s  description  from  D2's. 


DIELPNCH 

1.  Purpose  -  to  produce  a  printout  oc  the  raw  data  and 
resulting  calculated  information,  and  punch  a  data  deck  suit¬ 
able  for  reading  graph-drawing  data  into  the  data  cell. 

2.  Input  -  a  manually  key  punched  deck  containing  the 
samples,  dimensions,  partial  pressures  of  water,  weight  in 
vacuum  and  a  humid  atmosphere,  and  frequency-conductivity- 
capacitance-vacuum  capacitance  data. 

3.  Output  -  the  printout  and  deck  referred  to  above. 

4.  Operation 

a.  Allocate  storage  for  us«.  i,y  the  subroutine  which 
adds  vacuum  capacitance  data  for  entries  lacking  it.  Read 
the  "vacuum  data"  frequencies  and  corresponding  capaci¬ 
tances  from  the  appropriate  number  of  cards.  Then  read 
a  "case  heading"  card.  Calculate  the  values  of  those  vari 
ables  particular  to  each  case.  Then  print  out,  and  have 
the  machine  punch  cards  containing  "heading"  values,  punch 
ing  only  those  essential  to  graph  drawing.  Read  the  first 
of  the  several  cards  containing  the  frequency-dependent 
data  of  each  case.  Test  to  determine  the  signaling  ef¬ 
fect  of  the  frequency,  capacitance  and  vacuum  capacitance 
values  just  read,  and  call  the  routine  which  calculates 
or  recalls  the  vacuum  capacitance,  if  none  was  read.  Then 
for  F>5,  i.e.,  a  data  value,  the  frequency-dependent  vari¬ 
ables  are  calculated.  If  the  vacuum- to-humid  capacitance 
change  turns  out  to  be  less  than  zero  (a  physically  un¬ 
realistic  event),  then  the  capacitance  change  is  assigned 


*From  the  practical  standpoint. 


a  zero  value.  The  output  is  completed  by  printing  and 
punching  the  frequency-dependent  data. 

b.  If  the  frequency  value  is  >5,  then  the  program's 
control  is  transferred  to  reproducing  the  frequency  card 
just  read,  propagating  the  appropriate  signal  to  "DIELCELL." 
If  the  frequency  value  equals  five,  then  execution  is 
ended,  after  printing  the  list  of  frequencies  and  vacuum 
capacitances  read  and/or  calculated. 

5.  Variables  names,  in  order  of  appearance: 

FD  -  The  set  of  data  frequencies  at  which  a  vacuum 
capacitance  has  been  experimentally  determined. 

CVD  -  The  set  of  vacuum  capacitances  experimentally 
determined . 

WR  -  The  number  of  milligrams  of  water  adsorbed  from 
the  atmosphere  per  gram  of  sample  weight,  the  Weight  Ratio. 

L0GPP0  -  The  logarithm  to  the  base  10  of  the  ambient 
partial  pressure  of  water  divided  by  the  saturated  partial 
pressure  of  water,  at  the  sample  temperature. 

IVS  -  The  reciprocal  of  the  capillary  volume  squared . 

NFR  -  The  number  of  frequencies  for  which  vacuum 
capacitance  data  has  been  experimentally  determined. 

T  -  The  sample  thickness. 

D  -  The  sample  diameter. 

P  -  Ambient  partial  pressure  of  water. 

P0  -  Saturated  partial  pressure  of  water  at  the  sample 
temperature . 

W  -  Weight  of  the  sample  in  the  ambient  atmosphere. 

W0  -  Weight  of  the  sample  after  baking  out  in  vacuum. 

TEMP  -  The  sample  temperature. 

SAMPLE  -  The  sample  number. 

C0  -  The  theoretical  or  "geometric"  capacitance  of 
the  sample. 

RH  -  The  relative  humidity  of  the  atmosphere. 

V  -  The  capillary  volume  of  the  sample,  equal  to  the 
number  of  cc  occupied  by  the  adsorbed  water. 
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F  -  The  frequency  of  alternating  current  applied  to 
the  sample,  generally  100  to  1  million  hertz,  inclusive. 

G  -  The  conductance  (micro-ohms)  of  the  sample  at  the 
above  frequency. 

C  -  The  capacitance  (picofarads)  at  the  given  fre¬ 
quency. 

CV  -  The  vacuum  capacitance  appropriate  to  the  par¬ 
ticular  frequency. 

TD  -  The  loss  tangent,  commonly  called  "tan  del." 

EP  -  e 1 ,  the  charging  permittivity. 

EDP  -  e",  the  loss  permittivity  of  the  sample. 

DC  -  The  change  of  capacitance  on  adsorption  of  water 
from  the  atmosphere. 


6 .  Data  deck . 


Card 

Columns 

Variable 

type 

Format 

(inclus ive  ) 

code  name 

Remarks 

1 

12 

1  &  2 

NFR 

The  number  of  frequencies 
at  which  vacuum  capaci- 

tance  data  has  been  taken. 

2 

3F1 0 . 0 

1-10 

FD  ( 1 ) 

The  highest  of  the  above 
frequencies . 

11-20 

FD  (  2  ) 

The  next  highest  fre¬ 
quency. 

21-30 

FD  (  3  ) 

Continue  punching  fre- 

• 

• 

quencies  in  decreasing 

• 

• 

order  until  all  "NFR" 

• 

frequencies  are  punched. 
Use  as  many  cards  as 
necessary . 

3 

8F10.0 

1-10 

CVD(l) 

The  vacuum  capacitance  for 
the  highest  freouencv. 

11-20 

CVD ( 2 ) 

Carry  on  punching  vacuum 

• 

• 

capacitances  until  each 

• 

• 

vacuum  capacitance  corres- 

• 

• 

ponding  to  the  respective 
frequency  is  recorded  on 
a  type  3  card  or  cards. 
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Card 

Columns 

Variable 

type 

Forma  t 

(inclusive ) 

code  name 

Remarks 

4 

8  F 1 0 . 5 

1-10 

T 

"Case  heading"  card. 

11-20 

I) 

21-30 

T 

31-40 

T0 

41-50 

W 

51-60 

W0 

61-70 

TEMP 

71-80 

SAMPLE 

5 

F20. 10 

1-20 

F 

This  is  the  frequency- 

21-40 

G 

dependent  card,  normally 

41-60 

CV 

having  the  frequencies 

61-80 

CV 

arranged  in  decreasing 

order.  There  are  usually 
10-15  of  these  cards  per 
case  . 


DIELCELL 

1.  Purpose  -  to  arrange  the  calculated  data  punched  by 
the  DIELPNCH  program,  In  the  data  cell,  such  that  all  tempera¬ 
ture,  mode,  relative  humidity  and  frequency  sequential  corres¬ 
pondences  are  maintained;  to  establish  the  values  of  other 
quantities  needed  as  data  di s t in gu i she r s  or  by  the  CALC0MP 
plotting  routines. 

2.  Input  -  the  deck  from  DIELPNCH  preceded  by  cards  con¬ 
taining  the  axis  names,  the  number  of  temperatures,  the  humidi¬ 
ties  at  each  temperature  and  mode,  and  interleaved  with  the 
sequence  numbers  of  the  special  frequencies  (section  3-b  above) 
specific  to  each  case  (temperature,  mode,  and  relative  humidity 
comb i na  t ion ) . 

3 .  Output 

a.  A  printout  of  the  number  of  frequencies  for  each 
case  and  an  optional  (5300  line)  namelist  printout  of  the 
contents  of  the  data  cell,  as  arranged  therein. 

b.  Properly  arraiged  data  sets  designed  for  conveni¬ 
ent  use  by  the  drawing,  programs. 
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4.  Operation 

a.  Allocate  storage  in  the  configuration  which  will 
be  used  for  the  drawing  programs.  Establish  namelists 
containing  the  variables  used  by  each  of  the  drawing  pro¬ 
grams.  Set  all  data  spaces  equal  tc  zero,  to  eliminate 
uncontrolled  and  undesired  information  from  the  data  set. 
Read  18  cards  with  the  axis  titles  centered  in  the  first 
20  spaces,  then  another  card  with  the  number  cf  tempera¬ 
tures  at  which  data  has  been  taken.  Read  a  card  on  which 
is  punched  the  number  of  relative  humidities  appropriate 
to  each  temperature-mode  combination  ("T-MC").  This  set 
of  numbers  must  be  determined  by  inspection  of  the  print¬ 
out  of  DIELPNCH . 

b.  The  program  reads  two  "case  heading"  cards.  The 
first  contains  the  sequence  numbers  of  the  "special"  fre¬ 
quencies  which  are  the  parameters  for  drawing  the  type 
4-9  graphs.  The  value  of  the  case's  relative  humidity 
and  the  number  of  milligrams  per  gram  are  then  read  from 
the  following  card,  which  was  the  first  card  per  case 
punched  by  DIELPNCH.  We  next  have  the  machine  establish 
the  count  of  frequency  cards  for  each  case.  After  count¬ 
ing  the  current  frequency's  sequence  number,  the  values 
of  frequency,  capacitance  change,  e',  e"  and  tan  del  are 
read  in  sequence.  As  before,  the  value  of  F  in  comparison 
with  5  is  used  to  control  the  next  operation.  If  F  is 
greater  than  5,  hence  a  data  value,  the  particular  space 
in  storage  involved  is  filled  with  the  logarithm  to  the 
base  10  of  the  frequency,  other  data  values  are  stored, 
the  frequency  counter  is  stepped  by  one,  and  another  fre¬ 
quency  card  is  read.  Eventually,  the  program  encounters 

a  card  whose  frequency  is  less  than  or  equal  to  5.  At  the 
end  of  each  case  except  the  last,  an  F  value  of  zero  is 
encountered.  This  causes  the  program  to  store  the  number 
of  frequencies  that  are  in  the  case  and  to  print  out  the 
temperature,  mode  and  relative  humidity  indices,  the  number 
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of  frequencies  and  the  relative  humidity  value  Just  cal¬ 
culated.  Then  the  program  adjusts  the  index  for  relative 
humidity,  mode  or  temperature,  as  appropriate,  and  returns 
to  read  the  two  case  heading  cards  and  more  frequency 
cards.  If,  at  the  end  of  the  data  deck  punched  by  DIEL- 
PNCH,  an  F  value  of  5  is  encountered,  then  the  number  of 
frequencies  for  the  last  case  is  entered  in  storage  and 
the  last  printing  of  the  above  indices,  number  of  frequen¬ 
cies  and  relative  humidity  is  made. 

c.  The  next  job  that  the  program  undertakes  is  to 
calculate  the  set  of  origin  location  values,  giving  the 
x  and  y  coordinates  of  each  origin  in  sequence.  (See  page 
7,  the  TD  program  write-up.  Report  No.  7.)  Then  the  en¬ 
tire  set  of  data  calculated  and  recorded  in  this  program 
is  transferred,  array  by  array,  from  the  machine's  core 
into  the  data  cell  for  use  by  the  drawing  programs.  Then 
a  card  bearing  a  signal  as  to  whether  to  print  the  name- 
list  containing  the  data  (as  it  appears  in  the  data  cell) 
is  read.  If  this  signal  equals  zero,  then  the  namelist, 
consisting  of  some  5300  lines  of  printout,  is  not  printed 
out,  and  execution  is  ended.  However,  if  the  signal  is 
not  equal  to  zero,  then  the  namelists  are  printed  before 
execution  is  ended. 

5.  Variable  names,  in  order  of  appearance,  not  already 
named  above. 

TITLE  -  the  set  of  axis  titles  arranged  according  to 
graph  type  and  axis. 

LSP  -  the  set  of  "special  frequency"  sequence  numbers 
specific  to  the  graph  A+9  parameter  frequencies  and  speci¬ 
fic  to  a  case. 

KM  -  the  set  of  maximum  number  of  cases  per  T-MC. 

NF  -  the  number  of  frequencies  per  each  case. 

0  -  the  set  of  origin  locations  with  respect  to  the 
one  preceding. 
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Card 

ly-pe 

1 


I  -  the  index  which  marks  a  temperature’s  sequence 
number.  Maximum  value  is  5. 

J  -  the  index  which  marks  the  mode's  sequence  number. 
Maximum  value  is  2.  A  set  of  values  of  1  and  J  identifies 
a  T-MC. 

K  -  the  index  which  narks  the  relative  humidity  se¬ 
quence  number,  per  T-MC.  Maximum  value  is  20.  A  set  of 
values  of  I,  J  and  K  identifies  a  case. 

L  -  the  index  which  marks  the  frequency  sequence  num¬ 
ber  within  each  case.  Maximum  value  is  25. 

NSP  -  the  index  identifying  one  of  the  six  special 
frequencies,  the  parameters  for  graph  types  4  through  9. 

NT  -  the  number  of  temperatures  at  which  data  has  been 
taken  for  a  particular  specimen. 

NRH  -  the  number  of  relative  humidities,  specific  to 
a  given  T-MC. 

0X  -  the  number  of  inches  in  the  x  direction  to  the 
next  origin's  location. 

0Y  -  the  same  for  the  y  direction. 

Nl.Nz  -  integers  assisting  set-up  of  the  origin  loca¬ 
tion  array. 

QPRINT  -  the  indicator  of  whether  to  print  namelists. 

6.  Data  deck 


Columns  Variable 

Forma_t - (inclusive)  code  name _ Remarks 


5A4 


5A4 


1-20 


1-20 


The  first  "1"  means  that 
(l,l,l-*5)  this  title  refers  to  an 
x  axis,  the  second  means 
that  it  refers  to  a  type 
1  graph,  and  1  -»  5  refers 
to  each  block  of  four 
alphabetic  characters 
contained  therein.  The 
actual  content  of  this 
title  is  "Log  Frequency , " 
centered  on  the  first  20 
space  block. 

Title  "hog  frequency"  for  plot- 

(1,2,  l-*5)  ting  against  c',  whereas 
the  above  title  was  for 
plotting  against  tan  del. 
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Card 

tjj>c 

Foraat 

Columns 

(inclusive) 

Variable 
code  name 

P.ema  rks 

SA4 

1-20 

Title 
(2, 9, 1-5) 

"Delta  capacitance,"  for 
use  on  the  vertical  axis 
of  the  last  graph  type. 

2 

11 

1 

NT 

The  number  of  tempera¬ 
tures  for  which  data  has 
been  taken. 

3 

2413 

1-3 

4-6 

7-9 

• 

• 

• 

KM ( J , F) 

J-l  ,2 

1-1  , NT 

This  is  the  maximum  num¬ 
ber  of  relative  humidi¬ 
ties  appearing  per  T-MC . 
Note  that  the  entries 
must  be  r i gh t - j us t i f  led 
in  the  3-colunn  block 
that  each  of  these  vari¬ 
ables  uses. 

4 

613 

• 

• 

1-3 

LSP 

(1 ,K,J,I) 

The  sequence  number  of 
the  first  special  fre¬ 
quency  ,  1 .  c . ,  1MHz . 

• 

I  <*P 

(2  *K , J  ,  I) 

Sequence  number  at  which 
500  KHz  is  found. 

16-18 

LSP 

(6  ,  K  ,  J  ,  I ) 

Sequence  number  of  the 
last  special  frequency, 
500  Hz,  for  the  Kth  rela¬ 
tive  hunldity  of  the 

I-jth  x-MC.  The  indices 
for  the  case  nay  be 
punched  as  a  guide  to 
easy  handling  of  the  deck 
are  any  columns  beyond  18 

5 

19X.5.1 

6X.F6.2 

20-25 

RH 

This  case's  relative 
humidity,  for  printout 
only . 

31-37 

WR 

(K.J.I) 

The  ng/g  value  which  is 
to  be  plotted  for  this 
case . 

This  card  is  the  case 
heading  card  punched  by 

DI  ELPN’CII . 

no 


Card  Coluens  Variable 

type  Format  (Inclusive)  code  nnae 


Rena  rks 


6  SFI2.3  1-12 


13-24 

25-36 

37-48 


49-60 


F  The  particular  frequency 

(L,K,«,1)  value  identified  with 

the  "I.”  petition  in  the 
sequence  f  o .  the  case 
identified  with  the  R 1 ^ 
relative  Humidity  and/or 
ni  ingraft,  p  •  r  gran 
value,  the  .1 1  '■  node,  and 
the  1  1  L  t  enpc  ra  t  u  re . 

DC 

(L.K.J.l) 

rr 

(L.K.J.l) 

edp 

(L.K.J.l) 

TO 

(L.K.J.l) 

This  is  the  frequency- 
dependent  card  punched 
by  DI  r.LPNCI! . 


Several  groups  of  types  4,  5  and  6  cardo  forr  the  bulk  of  the 
data  deck  for  D1ELCELL. 


7  79X.F1.0  80  OPR  I  XT  If  the  value  of  QIM’IST 

is  zero,  then  no  nanc- 
litt  will  be  printed 
out  . 


REPAIR 

1.  Purpose  -  to  do  the  software  manipulating  required 
to  change  Isolated  values  or  grouns  of  values  of  the  data  in 
the  data  cell. 

2.  Input  -  the  changed  values,  both  bv  none  end  proper 
magnitude.  See  this  paragraph  for  subroutine  CHARCR.  immedi¬ 
ately  f  e  1  loving. 

3 .  Out pu  t 

a.  The  revised  values  to  the  proper  location  in  the 

data  cell. 

b.  Optional  printout  of  the  old  and  now  values. 

Details  arc  discussed  In  the  section  for  CHA.'h'.R. 


Ill 


4.  Operation  -  establish  blocks  of  main  core  storage  for 
use  by  both  REPAIR  and  CHANCE.  Read  the  entire  contents  of 
the  data  sets  designed  for  the  D1  and  D2  programs'  use.  Call 
the  subroutine  which  actually  changes  the  values,  CHANCE.  Upon 
return  from  CHANCE,  rewind  the  data  sets,  write  the  proper 
values  therein,  and  finish  execution. 

5.  Variable  names  are  the  same  as  in  the  DTELCELL  data 

deck . 

CHANCE 

1.  Purpose  -  assign  proper  values  to  space*?  in  the  data 
set  which  were  given  Incorrect  contents  when  established  during 
the  operation  of  DIELCELL.  Manipulation  of  the  data  set  (in 
the  magnetic  disc  "data  cell"  or  a  magnetic  tape)  is  done  by 
the  REPAIR  main  routine  above. 

2.  Input  -  this  may  be  done: 

a.  by  direct  arithmetic  assignment  statements, 

b.  calculated  arithmetic  assignment  statements,  or 

c.  rending  exact  or  functionally  related  values  from 
cards . 

3.  Output  -  the  revised  values  automatically  placed  in 
the  main  core  storage  shared  with  REPAIR. 

4.  Method  of  operation  -  the  possibilities  are  best  de¬ 
scribed  through  an  example.  Suppose  errors  in  transcribing 
data  from  notebooks  to  data  coding  sheets,  misreading  of  in¬ 
strument  ranges,  keypunch  errors,  or  svstem  malfunctions  have 
resulted  in  the  following  errors: 

a.  The  7th  through  12th  m i 1 1 i gr ans / g r am  values  for 
the  second  temperature,  adsorption  mode  run  are  0.312 
mg/g  too  high. 

b.  The  values  for  the  c'  and  r"  for  the  third  tem¬ 
perature,  the  desorption  mode  and  the  ninth  relative 
humidity  and  sixth  frequence  in  sequence  arc  missing. 

That  t '  value  is  ha,  -calculated  to  be  27.372,  and  the 
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value  of  TD  (6,9,2, 3)  appears,  is  checked  and  found  to  be 
correct.  Note  we  also  know  that  tan  del  *  e"/e'. 

c.  All  24  log  frequencies  for  the  seventh  relative 
humidity  and  the  first  temperature  adsorption  run  are 
nonsystematically  incorrect. 

d.  To  correct  these  errors  a  satisfactory  procedure 
would  be: 

(1)  use  iterative  arithmetic  assignment  statements 
to  subtract  0.312  from  WR  (7  12,1,2). 

(2)  assign-by-equation,  e"  (6, 9, 2, 3)  =  27.372 
TD  (6, 9, 2, 3). 

(3)  read  the  correct  frequencies  from  a  card  deck 
and  store  their  (common)  logarithms.  The  part  of  the 
original  deck,  extracted  for  this  job,  containing  the 
proper  frequencies  would  be  suitable,  and  eliminate 
any  need  for  additional  key  printing. 

e.  Since  it  is  advisable  to  have  visible  assurance 
that  the  changes  actually  took  place,  the  old  and  new 
values  of  interest  should  be  printed  out. 

5 '  Variable  meanings  -  the  same  as  in  DIELCELL. 

6.  Data  deck  -  it  is  understood  that  the  data  deck 
structure  listed  here  refers  to  the  above  example  only.  The 
needs  of  other  jobs  may  result  in  no  data  deck  or  data  decks 
of  several  hundred  cards. 


Card  Columns  Variable 

type  Format  (inclusive)  code  name _ Remarks 


F20. 10 


1-20 


(L  ,  7  , 1  ,  1  ) 


Listed  from  the  proper 
place  in  the  DIELCELL 
data  deck,  24  cards  rc- 
qui red . 


D 1 

1.  Purpose  -  to  draw  report- ready  graphs  of  tan  del,  e  1  , 
c  ,  at  various  r-XIC's,  drawing  three  parametric  curves  at  three 
relative  humidities. 
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2.  Input  -  the  above  data  set;  a  single  card  to  control 
an  optional  printout  of  the  data. 

3 .  Output 

a.  Up  to  210  report-ready  graphs. 

b.  A  printout  to  show  the  progress  (and  propriety) 
of  execut ion . 

4.  Operation  -  there  are  seven  subroutines,  which  do  the 
following  jobs : 

a.  MAIN  (the  name  will  be  used  interchangeably  with 

D 1  )  controls  the  allocation  of  storage,  the  temperature - 
mode  flow  through  the  data,  and  the  flow  through  the 
three  graph  types  for  each  T-MC.  MAIN  also  initializes 
the  graphing  sheet,  the  three  sets  of  points  (nascent 
curves)  for  each  graph,  and  sets  the  position  of  the  ori¬ 
gin  for  the  next  graph. 

b.  SETF3  converts  the  dielectric  quantities  into 

abscissa  and  ordinate  variables  for  graphing  All  further 
manipulations  of  the  data  will  be  made  in  terms  of  curve 
number,  first  through  third,  and  point  number  along  the 
curve.  i 

c.  JUSTFY  determines  whether  there  is  data  adequate 
to  justify  continuing  on  to  draw  a  graph.  It  also  deter- 

the  range  of  the  data  values  for  the  present  curves. 

d.  SCALE  adjusts  the  range  of  the  variables  to  make 
drawn  axes  both  efficient  and  easy  to  read  displavers  of 
the  data.  That  this  is  done  for  a  set  of  curves  is  the 
main  difference  between  the  TD  programs  and  the  DIEL  (see 
Report  No.  7  for  ID  program).  This  routine  is  a  much  am¬ 
plified  version  of  the  TD  "SCALIN’"  routine. 

e.  PREP,  as  in  TDDRAW,  draws  titled  and  scaled  axes, 
and  calculates  and  prints  enougli  tic  values  to  implement 
easy  readability. 

f.  POINT  sets  the  curves  on  the  axes,  by  drawing  a 
triangle  centered  on  the  curve  l's  points,  a  square  on 

i 

curve  2  '  s  ,  anti  a  diamond  on  curve  3's. 
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g.  SHFL  arranges  the  sequence  of  point  drawings  to 
reduce  graph  drawing  time. 

h.  Details  of  the  operation  of  D1  are: 

(1)  allocate  storage,  then  read  a  card  which  con¬ 
tains  a  signal  whether  to  print  out  the  contents  of 
the  data  set  used  by  the  program.  If  the  last  column 
of  the  data  card  is  blank  or  contains  a  zero,  the  data 
set's  contents  will  not  be  printed  out.  Initialize 
the  CALC0MP  package,  the  graphing  sheet,  the  graph 
position  number  (equal  to  zero)  aid  the  graph  position 
advancing  criterion  number  (equal  to  one). 

(2)  establish  and  iterate  the  temperature  and  mode 
indices  so  that  all  temperature-mode  combinations  are 
used  in  order  (i.e.,  adsorption,  desorption, 

adsorption,  desorption,  etc.).  Initialize  the 
graph  type  number  (equal  to  zero),  and  assign  a  scalar 
value  to  the  maximum  number  of  relative  humidities 
proper  to  this  T-MC .  If  there  are  no  relative  humidi¬ 
ties,  iterate  the  temperature  and/or  mode  index  and 
start  again  at  the  second  sentence  of  this  paragraph. 
For  a  finite  number  of  relative  humidities,  add  one  to 
the  graph  type  number,  and  set  the  per-graph  relative 
humidity  control  index  to  zero. 

(3)  if  the  graph  position  advancing  criterion 
number  is  not  zero,  a  new  set  of  axes  needs  to  be 
drawn,  and  the  graph  position  number  is  stepped  by  one. 
In  case  the  advancing  criterion  number  had  previously 
been  set  to  zero,  it  is  now  reset  to  one,  since  we 
assume  that  a  new  graph  position  will  be  used  the  next 
time  through.  Then,  the  data  from  previous  sets  of 
points  is  wiped  out  by  setting  the  number  of  points 
index- fo r-each-o f- the - th ree-suc ce e ding-c urve s ,  and  all 
75  available  abscissa  and  ordinate  values  equal  to 
zero.  Then,  the  origin  is  moved  to  its  next  location. 
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(4)  if  the  graph-advancing-criterion  is  zero  (and 
here  we  glimpse  the  functioning  of  the  later  routines), 
the  preceding  steps  have  been  taken,  but  a  graph  has 
not  been  drawn,  due  to  there  being  no  data  to  plot. 
Therefore,  the  old  graph  position  and  the  zeroing  of 
the  curve-point-locacion  variables  suffices  for  the 
present  case. 

(5)  write  a  message  in  the  printout  showing  the 
values  of  the  graph  type  index,  the  graph  position 
number,  and  the  temperature-mode  indices.  Then  call 
SETF3,  sending  the  appropriate  array  (TD,  EP  or  EDP) 
for  use  as  ordinate  values  for  the  graph  about  to  be 
drawn,  according  to  the  graph  type  number. 

(6)  as  seen  in  (4),  there  can  be  a  "good"  return 
from  SETF3  or  a  "bad"  one,  depending  upon  whether  the 
graph  has  been  drawn  or  not.  A  "bad"  return  causes 
the  T-M  indices  to  be  iterated.  "Good"  returns  from 
type  1  and  2  graphs  are  handled  by  testing  whether  the 
curves  have  been  drawn  for  all  relative  humidities. 

If  yes,  then  control  is  sent  back  to  where  the  graph 
type  index  is  stepped  (last  sentence  of  paragraph  (2) 
above),  so  that  it  eventually  becomes  three.  A  "good" 
return  from  SETF3  after  drawing  the  type  3  graph  again 
elicits  the  test  for  having  drawn  curves  at  all  rela¬ 
tive  humidities.  A  "no,"  meaning  that  there  are  more 
graphs  to  be  drawn,  sends  control  back  to  the  top  of 
(4). 

(7)  when  graphs  for  all  T-MC  have  been  drawn,  the 
CALC0MP  package  is  emptied  and  closed,  and  the  program 
stops  execution,  writing  a  message  in  the  Operating 
System  Monitor's  log  file,  showing  that  the  program 
has  enled  according  to  design. 

List  of  symbols  not  already  introduced. 

S  -  the  set  of  tic  spacings. 

ZS  -  the  set  of  scale  factors. 
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SD  -  the  set  of  units  of  axis  variable  per  tic  divi¬ 
sion,  henceforth  called  "tic  division." 

NPR  -  the  ordinal  number  of  tics  between  those  which 
are  marked,  i.e.,  print  a  label  every  "NPRt*1"  tic. 

X  -  the  array  containing  the  values  of  the  abscissas 
for  each  point  on  each  curve. 

Y  -  the  similar  array  for  the  ordinates. 

ZMAX  -  starts  out  being  the  largest  data  value  to  the 
graph,  per  axis. 

ZMIN  -  similarly,  except  the  smallest. 

NPTS  -  the  array  containing  the  number  of  points  on 
each  of  the  three  curves  for  each  graph. 

NG  -  the  graph  position  sequence  number,  or  merely 
"graph  number."  Maximum  value  =  210. 

NGADVQ  -  the  index  which  controls  the  advancing  of  the 
graph  position  sequence  number.  Mnemonic  with  the  ques¬ 
tion  "Should  NG  advance?" 

NTP  -  the  index  for  the  type  of  graph. 

KST  -  the  relative  humidity  sequence  number  which 
starts  the  parametric  curves  on  a  given  graph. 

KND  -  likewise,  except  the  end. 

NS  -  the  curve  number  on  a  particular  graph. 

6.  Data  deck 

Card  Columns  Variable 

type  Format  (inclusive)  code  name _ Remarks _ 

1  79X.F1.0  80  QPRINT  If  the  value  of  QPRINT 

is  zero,  no  namelist 
will  be  printed  out. 


Subroutine  SETF3 

Purpose  —  to  fill  generalized  abscissa  and  ordinate 
variables  with  those  values  of  tan  del,  e',  e",  and  log  F, 
which  will  draw  a  maximum  of  three  distinguishable  curves  on 
the  set  of  prepared  axes. 

I 

1 

1 

1 
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2.  Input  -  arguments  in  the  subroutine  list. 

a.  The  ordinate  array. 

b.  The  sequence  number  of  the  starting  relative 
humidity . 

c.  The  current  graph  number. 

d.  The  type  number  of  the  graph. 

e.  The  temperature  sequence  number. 

f.  The  mode  sequence  number. 

g.  A  dummy,  not  used  on  input,  variable. 

h.  The  number  of  relative  humidities  pertaining  to 
this  T-MC. 

i.  A  dummy  variable. 

j.  Also  implied  as  input  are  the  various  scaling  and 
graph  drawing  variables  which  are  in  the  common  storage 
shared  by  all  subroutines  of  Dl.  The  SETF3  common  block 
of  storage  implicitly  makes  input  of  the  log  frequency 
set  of  values,  and  the  number  of  f requencies-per-case 
values . 

3.  Output  -  the  normal  method  of  operation  makes  the  out¬ 
put  of  SETF3  go  to  subroutine  JUSTFY.  Output  to  this  routine 
includes  the  graph  sequence  number,  the  graph  number,  the  total 
number  of  curves  which  are  to  be  drawn.  On  return  to  MAIN, 
however,  if  the  number  of  frequencies  should  be  zero,  then  the 
last  numerical  argument  in  the  argument  list  is  given  a  value 
of  zero,  while  the  first  dummy  argument  (the  seventh  argument 
in  the  argument  list)  is  given  the  sequence  number  of  the  last 
relative  humidity  which  was  Used  in  curve  drawing. 

4.  Operation  -  the  first  time  through  per  T-MC;  the 
starting  relative  humidity,  sequence  number  is  given  a  value  of 
1  and  the  finishing  relative  humidity  sequence  number  is  given 
a  value  which  is  the  least  of  three  or  the  number  of  relative 
humidities.  However,  on  the  second  and  subsequent  times 
through,  the  starting  relative  humidity  sequence  number  is  in¬ 
creased  by  3.  The  ending  relative  humidity  sequence  number  is 
given  a  value  vhich  is  the  least  of  three  greater  than  the 
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value  it  had  the  last  time  or  the  total  number  of  relative 
humidities  for  the  T-MC.  Then  the  curve  number  index  is  set 
at  zero,  the  number  of  curves  is  calculated  and  for  each  rela¬ 
tive  humidity,  the  curve  sequence  number  is  increased  by  1. 

The  number  of  points  are  in  the  NSth  curve  and  the  abscissa, 
and  ordinate  point  locations  are  filled,  the  abscissas  with 
the  proper  value  of  log  frequency  and  the  ordinate  with  the 
proper  value  of  the  input  o r d ina t e - t ype  array.  After  the  re¬ 
quired  number  of  curves  has  had  its  abscissa-ordinate  loca¬ 
tions  determined,  the  program  writes  a  message  giving  the  range 
of  the  relative  humidity  sequence  numbers  and  the  number  of 
curves,  then  the  program  is  sent  to  subroutine  JUSTFY. 


5.  List  of  symbols  not  already  named. 

0RD  -  a  dummy  array  name  set  up  to  receive  the  values 
of  one  of  the  ordinate  arrays  as  named  in  MAIN. 

M  -  a  dummy  variable  whose,  value  is  given  to  NGADVQ 
on  return  to  MAIN. 

NCVS  -  the  number  of  curves  to  be  drawn  on  a  specific 
graph . 


NPNTS  - 
lar  curve . 


the  number  of  points  to  be  drawn  on  a  pai ticu- 


Subroutine  JUSTFY 

1.  Purpose  -  JUSTFY  is  designed  to  justify  carrying  on 
to  do  the  work  which  is  required  to  scale  and  draw  ticked  and 
labeled  axes.  The  specific  justification  criterion  is  that  at 
least  one  of  the  points  have  a  nonorigin  location,  inferring 
that  there  is  indeed  data  to  be  plotted.  In  the  process, 
JUSTFY  determines  the  largest  and  smallest  values  of  both  the 
abscissa  and  ordinate  variables. 

2.  Input 

a.  The  graph  position  sequence  number,  the  graph 
type  nu*  >er,  and  the  number  of  curves  which  are  to  be 
drawn . 
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b.  The  ten  arrays  in  that  block  of  storage  common  to 
all  other  routines. 

3.  Output 

a.  Normal  operation  -  JUSTFY  sends  to  subroutine 
SCALE  the  values  of  NG,  NTP  and  NCVS,  as  well  as  the 
largest  and  smallest  data  values  per  axis  via  the  block 
of  blank  common  storage. 

b.  "Bad"  operation  -  if  no  data  is  found  by  JUSTFY, 
control  returns  to  the  statement  in  subroutine  SETF3, 
which  sets  the  graph-numb er-advancing-criterion  integer 
equal  to  zero. 

4.  Operation  -  allocate  storage.  For  the  x  axis,  ini¬ 
tialize  the  range-end  variables  so  that  they  are  certain  to  be 
changed  by  a  data  value.  Determine  the  highest  and  lowest 
values  of  x  for  the  three  curves.  Repeat  this  process  for  the 
variables  on  the  y  axis.  The  justification  for  proceeding  is 
that  both  axis  maximum  values  be  different  from  zero  (see  para 
graph  b,(4)  of  section  Dl) .  If  justification  succeeds,  write 
the  largest  and  smallest  data  values  per  axis,  and  proceed  to 
subroutine  SCALE.  If  justification  fails,  write  a  failure 
message  and  return  to  MAIN  by  the  route  which  precludes  advanc 
ing  the  graph  position  number. 

5.  Flowchart 


Subroutine  SCALE 

1.  Purpose  -  SCALE  selects  interval  bounds  which  minimize 
both  dead  space  on  the  axis  and  distraction  by  unusual,  overly 
significant  tic  markings. 

2.  Input  -  the  numbers  of  the  current  graph,  its  type  and 
the  number  of  curves. 

3.  Output  -  proper  values  of  the  scaling  variables  (the 
first  four  of  the  blank  common  blocks),  as  well  as  revised 
values  of  the  high-  and  low-end  variables  for  each  axis. 


120 


4.  Operation 

a.  To  the  storage  allocations  which  have  been  used  so 
far,  add  thosje  for  the  array  of  axis  lengths.  Then,  for 
each  axis  in  turn,  convert  the  maximum  and  minimum  array 
values  to  scalars  for  speedier  referencing.  Determine 

the  order  of  magnitude  of  the  bounding  data  values  and 
their  ratio.  If  the  minimum  data  value  should  be  zero, 
assign  the  maximum  data  value  to  be  the  ratio  (as  though 
the  minimum  data  value  was  one).  Next,  determine  what 
fraction  of  the  order-of-magnitude  is  used  by  the  maximum 
axis  value  and  assign  a  proper  upper  bound,  from  a  selec¬ 
tion  of  six.  Similarly,  select  an  appropriate  value  for 
the  lower  bound  of  the  axis. 

b.  The  work  of  scaling  the  current  axis  is  done  in 
the  following  order: 

(1)  calculate  the  number  of  axis  variable  units 
per  tic  division,  then  the  number  of  tics  between 
labels.  Assure  that  the  axis  bounding  values  are  inte¬ 
gral  multiples  of  the  tic  division. 

(2)  test  to  see  that  more  than  six  tics  will  be 
drawn.  If  not,  halve  the  tic  division,  and  check  to 
be  sure  that  the  axis  range  is  nonzero.  If  it  should 
be,  return  to  SETF3  on  a  second  "bad"  route.  When  both 
these  tests  are  passed,  the  Scale  factor  and  the  number 
of  inches  of  actual  tic  spacing  may  be  calculated. 

When  both  axes  have  been  scaled,  the  success  message, 
showing  the  axis  bounds  and  the  scaling  variables  val¬ 
ues,  is  printed  out.  Control  then  goes  to  PREP. 

5.  Symbols  not  already  introduced. 

R  -  The  set  of  axis  lengths. 

NAX  -  The  axis  number;  "1"  for  abscissas,  "2"  for  ordi¬ 
nates. 

ZX  -  A  specific  value  for  the  maximum  value  on  an  axis. 
ZN  -  The  similar  minimum. 
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A  -  The  order  of  magnitude  one  greater  than  the  maxi¬ 
mum  data  value . 

B  -  The  similar  for  the  minimum  data  value. 

RATIO  -  The  maximum  data  value  divided  by  the  minimum 
data  value. 

N  -  An  integer  which  assists  in  determining  the  frac¬ 
tion  of  its  order  of  magnitude  that  the  data  variable 
actually  occupies. 

NTHI  -  The  number  of  tics  required  to  exceed  by  one 
the  maximum  data  value. 

NTL0  -  The  number  of  tics  required  to  bound  the  mini¬ 
mum  data  value  on  the  low  side. 

RANGE  -  The  actual  magnitude  of  the  axis  variable  which 
is  to  be  on  the  axis. 

Subroutine  PREP 

1.  Purpose  -  to  cause  CALC0M?  drawing  of  a  set  of  ticked 
and  labeled  axes,  printing  enough  tic  values  for  easy  readability. 

2.  Input 

a.  The  graph  number,  type  and  the  number  of  curves. 

b.  The  variables  in  the  blank  common  block. 

3.  Output 

a.  The  ticked  and  labeled  axes  on  the  graph  sheet. 

b.  PassLng  subroutine  P0INT  the  number  of  curves,  as 
well  as  the  abscissa  and  ordinate  variable  values  in  the 
blank  common  block. 

A.  Operation 

a.  Allocate  storage  as  before,  with  additional  stor¬ 
age  for  a  one-dimensional  array  containing  the  contents 
of  a  specific  title.  For  each  axis  in  turn,  assign  scalar 
values  to  array  variables.  Assign  values  to  variables 
which  locate  the  titles,  and  control  axis  and  title  draw¬ 
ing  angles.  Next  draw  the  axis  (subroutine  AXIS2 )  ,  then 
write  out  the  axis  title. 
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b.  Then  the  tic  labeling  parameters  are  determined. 

The  origin  tic  value  and  the  number  of  tics  arc  calculated 
and  the  tic  at  which  labeling  is  to  start  is  assumed  to  be 
the  first  to  be  labeled.  The  origin  tic  value  is  tested 
as  to  whether  it  is  an  integer  multiple  of  both  the  tic 
division  and  the  labeling  ordinal.  If  not,  the  second  or 
subsequent  tic  value  is  determined  and  the  same  test  is 
given  to  it.  Each  time  this  process  is  repeated  the  tic 
number  at  which  ticking  is  to  be  started  is  increased  by 
one.  Eventually  the  program  finds  a  value  which  is  an 
integral  multiple  of  both  tic  division  and  labeling  ordinal 
and  so  the  tic  number  for  starting  labeling  takes  its  value 
for  that  graph. 

c.  If,  however,  the  first  tic  will  be  more  than  one- 
third  the  way  down  the  axis,  the  tic  division  par,a;:is® ter  is 
halved,  and  the  other  scaling  parameters  are  adjusted. 

Then  the  scaling  process  is  repeated.  A  message  is  writ¬ 
ten,  telling  of  the  change.  Next,  depending  on  the  value 
of  the  tic  division  parameter,  the  labeling  ordinal  is 
fixed,  i.c.,  NPR'*"* (NAX)  or  N'MARK  ■  5  means  "print/mark” 
every  fifth  tic's  label. 

d.  The  tic  labeling  iteration  parameters  having  been 
determined,  at  each  appropriate  tic 

(1)  determine  the  tic  values. 

(2)  convert  the  numeric  tic  value  to  an  alphabetic 
label,  for  use  by  the  label  drawing  routine. 

(3)  determine  the  location  of  the  tic  value  ac¬ 
cording  to  its  value,  axis  and  number  of  characters, 
so  that  the  tic  label  will  be  centered  upon  its  proper 
tic. 

(4)  the  program  branches  here  to  accommodate  the 
printing  of  short  tic  labels  (four  characters  or  less) 
or  long  ones  (five  to  eight  characters).  When  all 
tics  on  both  axes  arc  drawn,  then  the  program  proceeds 
to  subroutine  P0INT.  The  PREP  success  message  gives 
the  value  of  the  last  tic  label  on  the  v  axis. 
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5.  Symbols  not  already  niaed: 

TL  -  A  one-diacnslonal  array  rontaining  a  specific 
title. 

SDE  -  Scalar  value  for  a  particular  axis'  tic  division 
value . 

NP  “  Scalar  value  of  the  nuabe r-be tween-pr int ed- t 1 c- 
values  array  ("labeling  ordinals"  array)  for  the  particu¬ 
lar  axis. 

SE  -  Tic  spacing  for  the  particular  axis. 

R  ~  Lengths  of  the  particular  axis. 

A  -  Angle  at  which  the  axis  is  to  be  drawn. 

XST  -  Abscissa  value  of  the  lower  left  hand  corner  of 
the  block  on  the  paper  which  contains  the  title. 

YST  -  The  similar  ordinate  value. 

AL  -  The  angle  at  which  the  lettering  in  the  title  is 
placed . 

NHALVE  -  The  nuaber  of  times  the  program  has  had  to 
halve  the  tic  division,  so  as  to  allow  for  an  adequate 
nuaber  of  marked  tics. 

TV0  -  The  tic  value  at  the  origin  or  other  first- 
labeled-tic. 

NTICS  -  The  number  of  tics  on  the  axis,  disregarding 
the  tic  which  is  placed  at  the  origin. 

NSTCNC  -  The  Number  of  tic  at  which  to  j>tarT  tiC 
label INC. 

Q  -  The  test  variable  which  determines  whether  a  tic 
value  is  an  integc,  multiplier  of  the  tic  division  value 
and  the  printing  ordinal. 

NQ  -  The  value  to  which  noninteger  Q's  arc  rounded  down. 

NMARK  -  The  scalar  tic  labeling  ordinal. 

N  -  The  sequence  number  of  tics  to  be  labeled. 

H  -  The  number  away  from  the  origin  of  a  tic  which  is 
to  be  labeled. 

TV  -  The  value  represented  by  a  particular  tic,  stored 
in  a  floating  point  form. 
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ITV  -  The  tic  value  raundod  down  to  the  next  lowest 
Integer. 

NC  -  The  number  of  characters  in  a  tic  label. 

T.TVIII  -  That  portion  of  an  integer  tic  value  between 
10,000  and  99,999,999. 

CTR  -  The  distance  a  tic  label  has  to  be  displaced  in 
the  -x  direction,  so  that  the  value  is  center'd  on  the  tic. 
This  is  proportional  to  the  number  of  characters. 

XI, Y1  -  The  abscissa  and  ordinate  values  of  the  lower 
left  hand  corner  of  the  block  of  paper  which  is  occupied 
by  the  tie  label. 

LAB  -  The  alphabetic  (binary  coded  digital)  content 
of  a  tie  label. 

LABII I  ,  LABL0  -  The  respective  portions  of  a  tic  label 
greater  and  less  than  10,000. 

Subroutine  P0IKT 

1.  Purpose  *  P0INT  marks  the  position  of  the  data  points 
on  the  prepared  set  of  axes. 

2.  Input 

a.  From  DREP  -  the  number  of  curves. 

b.  From  the  block  of  blank  common  storage,  the  values 
of  the  abscissa  and  ordinate  variables,  and  their  scale 
factors,  and  origin  values. 

3.  Output  is  the  set  of  points,  arranged  in  up  to  three 
curves,  distinguished  by  using  different  symbols  for  the  points 
of  each  curve. 

A.  Operation  -  Storage  allocation  consists  of  the  block 
of  blank  common.  For  each  curve  in  sequence,  check  to  deter¬ 
mine  if  its  points  should  be  reversed  in  order.  (See  subrou¬ 
tine  SHFL.)  Determine  the  CALC0MP  symbol  number  and  height 
appropriate  to  the  curve  number,  then  calculate  the  actual  x 
and  y  coordinates  (in  inches)  of  the  center  of  a  point  about 
which  the  symbol  will  be  drawn.  The  symbol  for  the  first  curve 
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is  a  triangle;  for  the  second,  a  square;  and  for  the  third, 
diamonds.  Actual  execution  of  the  drawing  is  done  by  the 
CALC0MP  routine  ’'TPRINT."  A  running  total  of  the  number  of 
points  is  kept,  and  after  all  three  curves  are  drawn,  the  suc¬ 
cess  aessage  gives  the  nunbe/  of  points  actually  drawn. 

After  the  third  graph  is  completed,  a  CALC0MP  message  is 
drawn  on  the  graph  sheet,  showing  the  order  of  drawing  the 
curve  s . 

5.  Symbols  not  already  named: 

NSY  -  The  symbol  number,  as  given  in  the  CALC0MP 

repertoire . 

HT  -  The  height  of  the  symbol. 

A,B  -  The  abscissa  and  ordinate  of  the  center  of  the 

point . 

NPT  -  The  total  number  of  points  drawn. 

Subrout ine  SHFL 

1.  Purpose  -  SHFL  changes  the  order  in  which  points  are 
drawn,  greatly  reducing  the  total  distance  that  the  CALC0MP 
pen  must  travel. 

2.  Input  -  SHFL  uses  the  curve  number  as  given  to  it  by 
subroutine  P0INT,  and  the  values  in  the  NPTS  array  (of  blank 
common)  so  as  to  properly  interchange  the  x  and  y  array  values. 

3.  Output  -  SHFL  returns  a  revised  sequence  of  x  and  y 
array  values  to  P0INT. 

4.  Operation  -  This  is  a  classic  algorithm  in  which  one 
of  the  values  to  be  interchanged  is  assigned  to  a  holding  vari¬ 
able.  Then  the  value  of  the  other  variable  is  given  to  the 
first,  after  which  the  value  of  the  holding  variable  is  as¬ 
signed  to  the  second.  This  process  is  cut  off  before  the  mid¬ 
dle  or  farther  point  in  the  original  sequence  is  changed,  and 
control  is  returned  to  P0INT. 
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5.  Symbols  not  already  named: 

M  -  The  sequence  number  of  the  "low-end"  point. 

NHI  -  Same  for  the  "high-end"  point. 

T1.T2  -  The  ordinate  and  abscissa,  respectively,  hold¬ 
ing  variable. 

D2 ,  MAIN 

1.  Purpose  -  to  draw  report  ready  graphs  of  e',  e",  and 
A  capacitance  vs.  mg/g  (the  WR  (K,J,I))  array  at  all  possible 
T-MC's.  Graphs  are  drawn  parametrically  at  1  MHz,  500  kHz  and 
100  kHz;  then  20  kHz,  5  kHz  and  500  Hz. 

2.  Input  -  The  data  set  established  by  DIELCELL,  plus  a 
single  card  to  control  an  optional  printout  of  the  data. 

3.  Output 

a.  Up  to  60  report  ready  graphs. 

b.  A  printout  to  show  the  progress  (and  propriety) 

of  execution. 

4.  Operation  -  There  are  seven  subroutines,  two  of  which 
(PREP  and  SHFL)  are  identical  with  their  counterparts  in  Dl. 
Three  others  are  virtually  the  same  as  in  Dl.  The  main  and 
setup  routines  are,  of  course,  quite  different. 

The  details  of  D2  MAIN's  operation  are:  allocate  storage, 
and  establish  a  namelist.  Read  the  data  from  the  data  set  and, 
if  desired,  write  its  contents,  depending  upon  the  just-read 
value  of  QPRINT .  Initialize  the  CALC0MP  package,  the  graph 
sheets  and  the  control  indices.  Determine  the  appropriate 
t e mp e r a t u re -mo d e  combination,  iterating  the  temperature  and 
mode  indices  so  as  to  cover  all  T-MC's. 

Set  and/or  step  the  graph-type  index,  and  set  a  new  origin 
when  indicated  by  a  nonzero  value  of  NGADVQ .  Print  a  message 
identifying  each  graph.  Eliminate  unessential  and  uncontrolled 
numbers  from  the  graph  drawing  block  of  storage.  A  computed 
"GO  TO"  statement  directs  the  program  to  set  up  one  of  the 
last  six  graph  types,  according  to  the  value  of  the  graph-type 
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index.  Note  that  provision  for  two  types  of  returns  from  sub¬ 
routine  SETL6  is  made.  After  all  T-MC ' s  have  been  graphed,  the 
CALC0MP  package  is  emptied  and  closed,  and  the  program  stops 
execution,  writing  a  message  in  the  Operating  System  Monitor's 
Log  File,  showing  that  the  program  has  ended  according  to  de¬ 
sign. 

Subroutine  SETL6 

1.  Purpose  -  to  fill  generalized  abscissa  and  ordinate 
variables  with  those  values  of  e' ,  e",  Ac  and  milligrams  per 
gram  needed  to  draw  curves  at  three  high  or  three  low  frequen¬ 
cies  on  a  set  of  prepared  axes. 

2.  Input  -  arguments  in  the  subroutine  statement's  list: 

a.  1)  a  generalized  ordinate  array,  which  can  hold 
the  contents  of  any  of  the  above  four-dimensional  arrays; 

2)  the  starting-special-frequency-number;  3)  the  current 
graph  number;  4)  the  type  number  of  the  graph;  5)  the 
temperature;  6)  mode  sequence  numbers;  and  7)  dummy  vari¬ 
ables  not  used  on  input. 

b.  Other  implicit  inputs  are  in  the  SET49  and  blank 
common  blocks. 

3.  Output  -  the  normal  method  of  operation  makes  the  out¬ 
put  of  SETL6  go  to  Subroutine  JUSTFY.  Output  to  this  routine 
includes  the  graph  sequence  number,  the  graph  number,  the  total 
number  of  curves  which  are  to  be  drawn.  On  return  to  MAIN, 
however,  if  the  case's  number  of  frequencies  was  zero,  then 

the  last  numerical  argument  in  the  argument  list  is  assigned 
a  value  of  zero. 

4.  Operation  -  allocate  storage  and  initialize.  Test  to 
determine  that  there  is  a  nonzero  number  of  relative  humidities, 
i.e.,  cases,  at  this  T-MC.  If  so,  for  each  of  the  special  fre¬ 
quencies  entered,  whether  of  the  high  frequency  or  the  low 
frequency  group,  and  for  all  the  relative  humidities,  test  to 
assure  that  data  have  been  taken  for  at  least  one  of  those 
frequencies.  If  not,  write  a  message  (only  once  per  T-MC) 
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that  there  is  no  data.  If  no  relative  humidity  for  the  T-MC 
has  data  for  the  frequency  group,  return  so  as  to  change  only 
the  graph  type. 

If  there  is  at  least  one  case,  and  at  least  one  point  of 
data,  assign  the  proper  values  for  the  particular  frequency 
applicable  to  the  case.  Form  the  sum  of  abscissa  and  ordinate 
as  a  check.  If  there  is  at  least  one  curve  with  at  least  one 
nonorigin  point  to  be  drawn,  proceed  to  Subroutine  JUSTFY. 

5.  Variables  not  already  named: 

Z  -  The  sum  of  data  point  coordinate  values. 

KBAD  -  The  count  of  relative  humidities  for  which  no 
high  (usually)  frequency  data  has  been  taken. 

NSP1,NSP2  -  These  values  are  the  remaining  values  in 
the  set  of  special  frequency  indicating  numbers. 
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//DI ELPNCH  JOB  (XX XX*XX, 010*06*  3000)  f  •  i  t  SBUOR  AYGER  H  t  *  ,  CL  A  SS  =  S 

/•PASSWORD 

//  EXcC  F4GCXS 

//FORT. SYS  IN  CD  * 

ccccccccccccccccccccctccccccccccccccccccccccccccccccccccccccccccc 

DI ELPNCH  CALCULATES  ALL  QUANTITIES  FOR  COMPLETE 
DIELECTRIC  ANALYSIS  GF  A  SAMPLE,  THEN  PRINTS  AID 
PUNCHES  THEM  IN  TEE  FORM  REQUIRED  FDR  DIELCELL. 

THERE  IS  NO  LIMIT  ON  THE  AMOUNT  OF  DATA  WHICH  CAN 
BE  HANDLED. 

TU  USE  DIELPNCH — 

KEYPUNCH  A  CARD  CONTAINING  THE  COUNT  OF  FRE¬ 
QUENCY  VALUES  FOR  WHICH  VACUUM  CAPACITANCE  VALUES 
HAVE  BEE  i  DETERMINED.  KEYPUNCH  (A)  CARDS  CONTAI¬ 
NING  ALL  THE  FREQUENCY  VALUES  APPEARING  IN  THAT 
LIST,  BEING  SURE  THEY  ARE  IN  DECREASING  ORDER. 

THESE  FREQUENCIES  SHOULD  BE  PUNCHED  EIGHT  TU  A  CARD 
IN  TEN  SPACE  BLOCKS,  TILL  ALL  ARE  RECORDED.  KEYPUNCH 
THE  CORRESPONDING  VALUES  JF  VACUUM  CAPACITANCE 
IN  ANOTHER  SET  OF  CARCS.  PLACE  THESE  CARL'S 
IN  FRCNT  OF  THE  CASE-HEADING  AND  FREQUENCY-DEPENDENT 
CARDS,  PUNCHED  FRCM  THE  CATA  CODING  SHEETS. 


ALLOCATIONS 

COMMON  FC( 130) ,CVO( ICO) , NFR 
REAL  LQGPPO , IVS 

READ  THE  VACUUM  CAPACITANCE  DATA. 

READ! 5 , I ) NFR 

1  FORMAT! 12) 

RE AD ( 5 , 2 ) ( F  D ( N ) , N=1 , NFR ) 

2  FORMAT ( 8 F 10 • 0 ) 

READ ( 5 , 2 ) ( C VU ( N ) , N= 1 , NFR  ) 

READ  A  CASE  HEADING  CARD,  AND  CALCULATE  THE  QUAN¬ 
TITIES  PARTICULAR  TU  EACH  CASE. 

3  READ (5, AIT , D,P,PO,W,WO, TEMP, SAMPLE 

4  FORMAT (oF 10. 5) 

Ca=8.654F-14*D*D»3. 1415V/4./(2.54*T) 

RH=10?.*P/P0 
LCGPPU=ALCGlv'(  P/PO) 

WR= 1000 . * ( w-WO )/ WO 

V=100.*„R 

IVS=1./V«V 
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noo  r>oo  n  o  rs  o  onoo  o  o  o  non  onoo 


PRINT  OUT  THE  CASE-PARTICULAR  QUANTITIES*  AND 

WRITE  COLUMN  HEADINGS  FOR  AC  FREQUENCY-DEPENDENT  DATA. 

WRI  TE(6,5)TEMP,P,.RH,WR.  SAMPLE*  L0GPP0,V,  I  VS 

5  FORMAT!  1H1///'  T  EMPERATURE=  •  F5 . 1 •  PARTIAL  PRESSURE  H20=' 
$F6. 2  *  '  RELATIVE  HUMIDITY  =  'F5.1,.'  MG  PER  G  =  '  F6.2/  '  AND' , 

1  *  SAMPLE  IS  NUKBER'Fb.l,  'LOG  P/PO  =  'F 6. 3 , 'CAPI LLARY  VOLUME' 

2  ,  *  =  ' F6 • 4  * '  IVS  =?  <  F8.1////3X,4HFRE0t5X,  5HOMEGA , 7X,4HC0ND , 
36X,3HCAP,6X,7HCAP  V AC , 3 X » 9H0ELT ACAP » 2X * 10HE PS  PRIME  • 
411HEPS  DEL  PRM , 2X » 7HT  AN  OEL ) 

47X ,  7HT  AN  OEL) 

PUNCH  THUSE  ESSENTIAL  TO  GRAPH  DRAWING. 
WRITE!7,6)TEMP,P,RH,WR 

6  FORMAT!'  T=  '  F4 ,  l  *•  *  +=*F5.2»'  RH='F5.I»'  WR  =  'Ffa.2) 

READ  A  FREQUENCY-DEPENDENT  CARD. 

7  READ!5,3)F,G,C,CV 

8  FORMAT ! 4F20 • 10 ) 

DETERMINE  THE  SIGNAL  APPROPRIATE  TO  THE  FREQUENCY 
AND  CAPACITANCE  VALUES. 

IF ! F-5. ) 15,15,9 

9  I  F ! C ) 7 , 7 , 10 

RETRIEVE  OR  INTERPOLATE  UNPUNCHED  VACUUM  CAPACITANCES 
BY  USING  ACOCV. 

10  I F ( CV ) 1 2 , 1 1 ,  12 

11  CALL  ADDCV ! F , CV ) 

CALCULATE,  ... 

12  0=2. *3. 14159*F 
TD=1.E6*G/U/C 
EP=1.E-12*C/C0 
EDP=EP*TD 
DC=C-CV 

IFIUC.LT. 0. )CC=0. 

...PRINT  OUT,  AND... 

WRI TE ! 6 , 1 3 ) F , 0, G, C, CV, DC, EP , EDP , TD 

13  FORMAT(/2F10.1,7F11. 3) 
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/ 


.  •  • PUNCH  THE  FREQUENCY- DEPEND ENT  VARIABLES,  AND  RE¬ 
TURN  TO  READ  ANOTHER  FREQUENCY-DEPENDENT  CARD. 

WRITE (7, 14  IF, DC, EP, EDP, TD 

14  FORMAT (  5F 1 2 • 3  I 
GO  TO  7 

FOR  F  LESS  THAN  OR  EQUAL  TO  5,  THE  CASE  HAS  ENDED. 
THE  SIGNAL  IS  PROPAGATED  TO  DIELCE^L  BY  THIS 
PUNCH  COMMAND. 

15  WRITE(7,0)F,G,C,CV 

THE  PROGRAM  IS  SENT  FOR  ANOTHER  CASE  CARD  WITH  F 
LESS  THAN  5,  OR  TO  THE  END  OF  EXECUTION  FOR  F  =  5, 
AFTER  PRINTING  THE  VACUUM  DATA. 

I F ( F-l • I  3,3,16 

16  PRINT, 17  (  FD(  N)  ,  CVDI  N  I  ,  N=.  1 ,  NFR  )  ,-NFR 

17  FORMAT (5(F10.0,F8.3)  ) 

STOP  1 

END 
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/ /OIELCELL  JOB  ( XXXX , XX, 0 10, 10, 1COO ) , » SiSLARRAYGER S$$ » , CLASS=M 

/•PASSWORC 

//  EXEC  F4GCXM 

//FORT.SYSIN  DC  * 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 


OIELCELL  SETS  UP  THE  DATA  PUNCHED  BY  DIELPNCH  INTO 
THE  FORM  USED  BY  THE  DRAWING  PROGRAMS.  TWO  MAG¬ 
NETIC-MEMORY  DATA  SETS  HOLD  THE  DRAWING  PROGRAMS* 
CATA. 


C 
C 
C 
C 
C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  ALLOCATIONS. 

c 


INTEGER  TITLE 

COMMON  LSP(6,20»2»5),KM(2»5)».TD(25,20»2»5),EP(25,20»2,5), 
*EDP(25,20,2,5),  DC( 25, 20,2,  5), F(25,  20,2,  5),  WR(20,2,5)  , 

INF (20, 2,5), 0(2, 210), TITLE! 2, .9, 5) 

NAMELIST/ONE/NT ,KM,NF, F,TD»  EP» EDP, 0, TI TLE 
NAMEI.I  ST/TWO/NT  ,KM,LSP,WR,  EP  ,.E  DP  ,.DC ,  0,  T I TL  E 
C 

C  CLEAR  THE  DATA  STORAGE  OF  UNCONTROLLED,  UNDESIRED 

C  INFORMATION. 

C 

DO  2  1=1,5 
DO  2  J=  1 , 2 
KM ( J , I ) =0 
DO  2  K=1 , 20 
DO  1  NSP=1 , 6 

1  LSP(NSP,K,J, I ) =0 

NF ( K, J ,  1 1  =  0 
WR(K, J , I ) =0 . 

DO  2  L=  1 , 2  5 
TD(L,K,J, I ) =0 . 

E  P  (  L  ,  K  ,  J  ,  I  )  =0  . 

EDP ( L , K , J , I ) =0 . 

DC  (  L,  K ,  J,  I  )=f). 

2  F(L,K, J, I )=0. 


FOR  THE  NEXT  FOUR  COMMENTS  ,  THE  PROGRAM  READS  THE  ... 


...€IGHTE5N  AXIS  TITLES,  A  PERMANENT  PART  OF  THIS 
PROGRAM'S  DECK. 


READ( 5 , 3 ) ( ( ( TITLE( NAX, NTP,N),N=l,5),NTP=l,9) ,NAX=1 ,2) 
3  FORMAT ( 5A4 ) 


...NUMBER  OF  TEMPERATURES,  AND  THE  NUMBER  OF  RELA¬ 
TIVE  HUMIDITIES  FOR  EACH  TEMPERATURE-MODE  COMBINATION. 
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READ! 5 , 4 )  NT 

4  FORMAT  (ID 

RE ADI  5, 5) (IKM( J, I) , J=l,2), 1=1, NT) 

5  FORMAT (1013) 

00  15  1=1,  .mT 
00  15  J  =  1 »  2 
NRH=KH( J,I) 

I F ( NRH ) 14, 14,6 

6  00  14  K=1,NRH 

...THE  TWO  CASE-HEADING  CARDS.  THE  PRECEEDING  THREE 
CARDS  CONTENTS  RESULT  FROM  INSPECTING  THE  DIELPNCH 
PRINTOUT. 

READ! 5, 7)  (LSP(  NSP,K, J,  I  ),NSP=U,6) 

7  FORMAT (613) 

READ!  5«8)  RH»  WK(K,JV.I  ) 

8  FORMAT ( 19X,  F5 . 1 , 6X , F6 . 2 ) 

t.=0 

COUNT  THE  NUMBER  OF  FREQUENCIES  PER  CASE. 

9  L  =  L  +  1 

REAC  A  FREQUENCY-DEPENDENT  CARD*  AND  EVALUATE 
ITS  SIGNAL. 

READ(5»10)F(L»K«.Jt  I  )  t  DC  ( L,  K,  J ,  I ) ,  EP  (  L .  K,  J ,  I  )  , 
$EOP(L*K,JtI )«TD(L*K,J, I ) 

10  FORMAT ( 5F12 .3) 

IF(5.-F(L,K,J,I  Ull, 16, 12 

FUR  F  GREATER  THAN  5,  I.  E.,  A  DATA  VALUE,  STORE  ITS 
COMMON  LOGARITHM,  AND  READ  ANOTHER  CARD. 

11  F(L,K,Jf I )= AL0G101  F ( L, K, J,  I  )  ) 

GO  TO  9 

IF  F  IS  LESS  THAN  5,  THIS  CASE  ENDS,  BUT  OTHERS 
FCLLCW . 

12  NF(K,J, I)=L-1 

WRI TE ( 6 , 13 )  I , J, K,NF( K, J,  I  )  ,RH 

13  FORMAT  (  10Xf  *  FOR  I,Jf  K  =  »I2,2(*,  M2),',  NF  =  »I2, 

$•  AND  «H  =  *  F5 . 1 , »  » ) 

14  CONTINUE 

15  CONTINUE 

WHEN  F  EQUALS  FIVE,  NO  MORE  DIELCTRICS  INFORMATION 
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16  NF ( K, J  ,  I ) =L- 1 

PRINT, 13  I,  JyKyNFlKyJ,  I  ) , RH 

THIS  LOOP  CALCULATES  THE  LOCATIONS  OF  THE  ORIGINS. 


0X  =  0. 

OY  =  8  •  S 

DO  17  N=  1 » 2  10,  3 
N1=N* 1 
N2=N+2 
0! 1 ,N)-OX 
0( i,Ni)=0. 

01  1,N2)=0. 

0(2 ,N) =0. 

0(2, N1 ) =0Y 
0(2,N2)=0Y 
0X=11. 

17  0Y=-0Y 

WRITE  THE  DATA  INTO  THE  DATA  CELL,  IN  SETS  DESIGNED 
FOR  USE  BY  D1  AND  02. 

WRITE!  1 )  NT  ,KM,0,.TITLE,NF,F;  TD,  EP,  EDP 
WRITE!  2)  NT,KM,0,.TITLE,LSP,WR,EP,  EDP,  DC 

READ  THE  CARD  WHICH  DETERMINES  WHETHER  THE  NAME- 
LISTS  ARE  TO  BE  PRINTED. 

KEAD(5,18)QPRINT 

18  FORMAT  ( 79X , FI .0 ) 

I F ( CPRiNT.EQ.O. )  STOP  1 
WRITE! 6 ,ONE ) 

WR I TE ( 6 , 19 ) 

19  FORMAT ( 1H1 ) 

WR I T  E ! 6 , TWO ) 

STOP  2 


THE  TITLE  CARDS  AND  JCL  BELOW  ARE  AN  UNVARYING  PART 
OF  THIS  PROGRAM  * S  DECK. 


//GO.FTOlFOOl  DO  OSN=THRE£ . D AY. FOR. Di , UN  I T*2 321, 

//  SPACE=(CYL, (10,5) ),VOL=f  LR=OCELL 1, DI SP=( NEW, KEEP ) 

/ /GO.  FT02F00 1  00  DSN  =  THR  EE .DAY  .FOR  *D2,UNITs232l, 

//  SPACE2 (CYL, (10,5) ), VOL =S ER=DCELL 1 , DI SP= ( NEW, KEEP ) 

OR  IF  TAPES  SHOULD  BE  REQUIRED,  ESTABLISH 
THEM  WITH  THE  FIRST  FOUR  CARDS  BELOW,  BUT 
FOR  SUBSEQUENT  REFERENCES,  BY  REPAIR,  DI  OR 
D2,  REPLACE  THE  2ND  AND  4TH  CARDS  BY,  RES¬ 
PECTIVELY,  THE  5TH  AND  6TH,  USING  THE 
VOLUME  SERIAL  NUMBER  SUPPLIED  BY  THE 
COMPUTING  CENTER. 


//GO.FTOlFOOl  OD  UNI T=TAPE9 , OCB= ( DEN=2 ) , LABEL=( l , SL ) , 
//  OSN«DRAW.ONE.TAPE,DISP=( NEW, KEEP) 

//GO.  FT02F00 1  DO  UNIT*.TAPE9,  DCB=  (  DEN=2  ) ,  LABEL»(  1 ,  SL)  , 
//  DSN2ORAW.TWO.TAPE,DISP  =  (  NEW,  KEEP ) 

//  DSN* DRAW .ONE. TAPE,DISP=(OLD,KEEP), V0L=SER*02XXXX 
//  DSN=DRAW.TWO.TAPE,DISP=( OLD, KEEP), VOL  *SER*02XXXX 
//GO. SYS  I N  DD  * 

LOG  FREQUENCY  NAX  AND  NTP  ARE 

LOG  FREQUENCY 


LOG  FREQUENCY 
MG  H2C  PER  G 
MG  H2C  PER  G 
MG  H2C  PER  G 
MG  H20  PER  G 
MG  H20  PER  G 
MG  H2C  PER  G 
TAN  DELTA 
EPSILON  PRIME 
EPSILON  DOUBLE  PRIME 
EPSILCN  PRIME 
EPSILON  PRIME 
EPSILON  DOUBLE  PRIME 
EPSILON  DOUBLE  PRIME 
DELTA  CAPACITANCE 
DELTA  CAPACITANCE 


1,1, 

1,2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 
2,1 
2,2 

2.3 

2.4 

2.5 

2.6 

2.7 

2.8 

2.9 


RESPECTIVELY. 


*  THE  REST  OF  THE  DATA  DECK. 

ft 


THE  l  IN  COLUMN  60  CAUSES  DATA  PRINTUUT.  1 

/*EOJ 
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SUBROUTINE  AC0CV(A,0) 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccct 

ACCCV  PROVIDES  A  VALLE  FUR  VACUUM  CAPACITANCE  WHERE 
NONE  HAVE  BEEN  PUNCHEO.  THE  VALUE  MAY  EITHER  OE  RE¬ 
CALL  EC  FROM  THE  ORIGINAL  SET,  OR  CALCULATED  BY 
LINEAR  INTERPOLATION. 


ALLOCATIONS 

COMMON  FC!I00),CV0(100),NFR 
CO  3  N=  1 » NFR 

THE  INPUT  FREQUENCY,  CALLED  '  'A*  *  BY  THIS  ROUTINE,  IS 
REPEATEDLY  COMPARED  WITH  THE  LIST  OF  DECREASING-OR- 
CEREC  FREQUENCIES,  UNTIL  ONE  SMALLER  OR  EGUAL  IS  FOUND. 

IFIFC(NJ-A)  1,2,3 

FOR  THE  ♦♦SMALLER'*  CASE,  THE  APPROPRIATE  VALUE  OF 
VACUUM  CAPACITANCE  IS  CALClLATED  BY  LINEAR  INTERPO¬ 
LATE. 

1  8*CVC(N-l)-(CV0(N-l  )-CVC(Nn»IFD(N-l)-A)/lFO(N-l)-FD(N)) 
B*FLCAT( IFIX(B*iOOC. J I/1COC. 

GC  TO  * 

FOR  THE  EQUAL  I T  Y  CASE,  THE  VACUUM  CAPACITANCE  IS 
THE  CORRESPONDING  ORIGINAL  ONE. 

2  B=CVCIN) 

RETURN 

3  CCNTINUE 

PUT  INTERPOLATES.  VALUES  INTO  THEIR  PLACE  IN  THE  LIST. 

4  P®  FC ( N ) 

R»CVC(M 
F  C  (  N  )  *  A 
CVC(M=C 

CC  5  M=N,NFR 
O®  FC ( M ♦ 1 ) 

S®CVC( M*  1 ) 

FC(MM)=P 

CVU(M*1)=R 

P®C 

5  R=S 
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noon 


SINCE  THE  LIST  HAS  BEEN  INCREASED  IN  LENGTH  BY  ONE, 
THE  INCREASE  IS  MAGE  EXPLICIT  HERE. 

NFR*NFR«1 

RETURN 

END 
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^/REPAIR  JOB  (XXXX»XXf010,01»0000 ) , » •  UBUDRAYGERt*** ,CLASS=M 

/•PASSWORD 

//  EXEC  F4GCXM 

//FORT.SYSIN  OD  • 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 

c 

c 

c 


REPAIR  OOES  THE  MANIPULATION  OF  THE  MAGNETIC 
OATA  SETS  WHILE  SMALL-  OR  LARGE-SCALE  ALTERATIONS 
OF  THEIR  CONTENTS  ARE  MADE. 


C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C 

C  ALLOCATIONS. 

C 

INTEGER  TITLE 

COMMON  LSP(6,20,2,5),KM(2, 5 ) , TD( 25 , 20, 2 , 5 ) ,  EP ( 25, 20 , 2 , 5 ) , 

$E0P(25,20,2,5),0CC25,20,2,5),FC25,20,2,5),WR(20,2,5) , 

INF (20, 2, 5), 0(2, 210), TITLE! 2, 9, 5), NT 
C 

C  READ  THE  CONTENTS  OF  THE  DATA  SETS. 

C 

READ  (L)  NT, KM, 0, TITLE, NF,F, TO, EP.EDP 
READ  (2)  NT, KM, 0, TITLE, LSP.WR, EP.EDP, DC 
C 

C  *  *  CHANGE  *  *  IS  UNIQUE  TO  EVERY  RUN,  AND  CORRECTS  ONLY 

C  THOSE  PORTIONS  OF  THE  OATA  SET  REQUIRED. 

C 

CALL  CHANGE 
C 

C  REWIND  AND  REFILL  THE  OATA  SETS. 

C 

REWIND  2 
REWIND  3 

WRITE! I )  NT, KM, 0,TITLE,NF,F,TD, EP.EDP 
WRITE12)  NT,KM,0,TlTLE,LSP,WR,EP,EUP,DC 
STOP 
END 


o  o  o  o  n 


I 


SUBROUTINE  CHANGE 

cccccr.ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


C  CHANGE  PERFORMS  THE  OPERATIONS  REQUIRED  TO  COR- 

C  RECT  BAO  PORTIONS  OF  THE  DATA  IN  THE  DATA  SET. 

C 

c  THE  OPERATIONS  OF  THE  PROGRAM  LISTED  BELOW  ARE  AN 

C  EXAMPLE  ONLY,  MEANT  TO  SHOW  THE  OPTIONS  AVAILABLE. 

C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  ALLOCATIONS 

C 


cc  ccccccccccccccccccccccccccccc 


COMMON  LSPI6,20,2,5),KM(2, 5 > , TD( 25 , 20, 2 , 5) . CPI  25 , 20 ,2 , 5 ) , 
SEOPI25, 20, 2, 5), OC (25, 20, 2, 5), FI25, 20,2,5), MR (20,2, 5) , 
INF(20,2»5),0(2,2I0)»TITLE(2,9,5),NT 
G 

C  PRINT  OUT  THE  INCOMING  VALUES. 

C 


WRI TE ( 6 , 1 ) 

1  FORMAT  (•  ON  ENTRY—*) 

WRITE(6,2)(WR(K,i,2),K=i,20),(F(L,7,i,  1),L=*1,25) , 
*EDP(6,9,2,3) 

2  FORMAT!*  WR(K,1,2)  =  ' /4( 5G24. 5 ) // *F I L , 7, 1 , I )  *  •/ 
J5I5G24.5)//* EDPI6, 9,2, 3)  =  *F20. 10//// / ) 

C 

C  CHANGE  THE  VALUES. 

C 

00  3  K*7, 12 

3  WR(K,I,2)*WR(K,1,2)-0.312 
EDP(6,9,2,3)*27.372*T0(6,9,2,3) 

DO  5  L*l,24 

READ! 5 • 4 )  A 

4  FORMAT ( F20. 10) 

5  F(L,7,1,1 )3 ALOGLOI A ) 

C 

C  PRINT  OUT  THE  CORRECTED  VALUES. 

C 

WRI TE ( 6, 6 ) 

6  FORMAT  I •  AFTER  CORRECTION*) 

WRITE(6,2)( WR(K,1,2),K*I,20), (F(L,7,I, I),L*1,25), 
$E0P(6,9,2,3) 

RETURN 


1A0 


o  n  o  o  n  o  n 


I 


c 

C  JCL  FOR  REPAIR/CHANGE. 

END 

/* 

//GO.FTOlFOOl  DO  DSN  =  THREE .DAY. FOR .D L , UN  I T  =  232 1 » 
//  VOL*SER=CCELLl »  DISP=OLD 

//GO.FT02F001  OD  OSN  =  THREE. DAY . FOR . D2 . UN  I T  =  2 32 1 , 
//  VUL*SER=OCELLl»OISP=OLD 

//GO. SYS  I N  UD  • 

••••♦the  data  deck*«*»* 

/*eoj 
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// D1  JOB  (XXXXtXX, 035,06#  1000), ' $SBUDRRAYGER$SS » , CLASS=M 

/•PASSWORD 

//  EXEC  F4GCXM 

//FORT •  SYS  I N  DD  * 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

C  D1  DRAWS  REPORT-READY  GRAPHS  OF  THE  PERMITTIVITIES  C 

C  AND  LOSS  TANGENT  VS.  LOG  FREQUENCY,  SCALING  THE  C 

C  AXES  FOR  EACH  GRAPH.  C 

C  C 

C  c 

C  THE  MAIN  PROGRAM  READS  THE  DATA,  AND  CONTROLS  THE  C 

C  FLOW  THROUGH  ALL  THE  RELATIVE  HUMIDITIES  UF  EACH  C 

C  TEMPERATURE-MODE  COMBINATION,  DRAWING  ALL  GRAPHS  C 

C  OF  A  GIVEN  TYPE  (VIZ.,  WITH  A  GIVEN  ORDINATE),  C 

C  BEFORE  PROCEEDING  TO  THE  NEXT  TYPE.  C 

C  TO  USE  Dl—  t 

C  DETERMINE  WHETHER  THE  DATA  PRINTOUT  IS  DESIRED.  C 

C  PUNCH  THE  DATA  CARD  ACCORDINGLY.  USUALLY,  NO  PRINT-  C 

C  OUT  IS  NEEDED,  CALLING  FOR  A  BLANK  OR  A  ZERO  IN  THE  C 

C  DATA  CARD'S  LAST  COLUMN.  C 

C  c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  ALLOCATIONS 

c 

INTEGER  TITLE, KM(2»5) 

REAL  ORDI25,20,2,5  ) , TD( 25, 20, 2 , 5) , EDP ( 25, 20, 2, 5 ) , EP 1 25 ,20, 2 
S , 5 ) ,Q( 2,210 ) »R( 2 ) 

COMMON  S(2),ZS(2),SD(2),NPR(2),X(3,25) ,Y(3,25) ,ZMAX(2) ,ZMIN 
$(2)»TITLE(2»9, 5)»NPTS(3) 

COMMON  /SETI3/F(25,20,2,5) ,NF{ 20,2,5) 

NAMELIST  /UNE/NT,KM,NF,F,TD,EP,EDP,Q, TITLE 
C 

C  READ  THE  DATA  FROM  THE  DATA  CELL,  AND  PRINT  THE 

C  CONTENTS  IF  DESIRED. 

C 

READ  (I)  NT,KM,0,T  ITLE,NF, F,TD,EP,EDP 
READ  (5,1)  QPRINT 
I  FORMAT  { 79X , FI .0 ) 

IF ( QPRINT .NE.O. )  WRITE(6,0NE) 

C 

C  INITIALIZE  THE  CALCOMP  PACKAGE,  THE  GRAPH  SHEETS,  AND 

C  THE  CONTROL  INDICES. 

C 

CALL  PLSTRT ( 1111*11) 

CALL  PLOT ( 3 . , I • 50 , -3 ) 

NG-0 

NSP=0 

NGADVQ-I 
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noon  o  ooo  ooooo  non  ooooo 


FOR  EACH  TEMPERATURE-MODE  COMB  I NA TI ON,  TEST  TO  SEE  IF 
THE  NUMBER  OF  CASES  IS  GREATER  THAN  ZERO,  AND  STEP 
THE  GRAPH  TYPE  INDEX. 


DO  14  1=1, NT 
DO  14  J=  1 , 2 
NT  P=0 

NRH=KM ( J , I ) 
IF(NRH)  14,14,2 
2  NTP=NTP+1 
KST =0 


SET  A  NEW  ORIGIN  WHEN  NECESSARY. 

3  I F ( NGADVQ )4, 6,4 

4  NG=NG*1 
NGADVQ=1 

ELIMINATE  ALL  UNESSENTIAL  AND/OR  UNCONTROLLED  NUMBERS 
FROM  THE  BLOCK  OF  STORAGE  WHICH  WILL  BE  FILLED  WITH 
GRAPH  DRAWING  INFORMATION. 


DO  5  NS=  1 , 3 
NPTS ( NS  1  =  0 
DO  5  N= 1 , 25 
X ( NS » N )=0 • 

5  Y ( NS ,N ) =0. 

CALL  PLOT (0( l,NG),0(2,NG),-3) 

PRINT  A  MESSAGE  IDENTIFYING  EACH  GRAPH. 

6  WRITE(6,7)NTP,NG, I  ,  J 

7  FORMAT  (  '  0  * ,  '  NTP  ='I2,',  NG  =  M3,',  I  =  M2,»,  J  =  M2) 

B  I  F  (  NTP-2  )  9,10,12 

PROCEED  WITH  SETTING  UP  THE  GRAPH  OF  THE  PRESENT  TYPE, 
WITH  THE  PROPER  ORDINATE  ARRAY  (FIRST  ARGUMENT)  AND 
FREQUENCY  RANGE  (SECOND  ARGUMENT). 

9  CALL  SETF3(  TD» KST, NG, NTP, I ,J,KND,NRH, NGADVQ,  14,  11) 

10  CALL  SETF3( EP, KST, NG, NTP, I ,J,KND,NRH, NGADVQ,  14,  11) 

11  I F ( KND-NRH)  3,2,2 

12  CALLSETF3(EDP,KST,NG, NTP, I, J,KND,NRH, NGADVQ,  14,  13) 

13  I F( KND-NRH)  3,14,14 

14  CONTINUE 
CALL  PLTEND 
STOP  111 
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CHOOSE  ONE  OF  THE  SAMPLE  DATA  CARDS  BELOW 


/ 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


**THE  JOB  CONTROL  LANGUAGE  CARDS  BELOW**** 
**  BETWEEN  ••/*»'  AND  ••/•EOJ'*  **** 
**ARE  PLACED  BEHIND  THE  DECK  FOR  SUB-  **** 
••ROUTINE  SHFL  WHEN  ACTUALLY  RUNNING  **** 
**  THE  PROGRAM,  **** 


END 

/* 

//GO.FTOIFOOI  DD  DSN = THREE. DAY, FOR .D1»UNIT=2321# 

//  VOL=SER=DCELLI,DISP=OLD 

/ /GO. FT18F00I  DD  UNI T=TAPE7 , DCB= ( RECFM=U , BLKS IZE= 1440 , DEN=1 , TRTCH=E T ) , 
H  DISP=OLD,LABEL={ * BLP), VOL=SER=OUTPT, DSN=CALCOMP 

//GO. SYS  I N  DD  * 

A  BLANK  IN  COLUMN  80  DELETES  DATA  PRINTOUT, 

THE  I  IN  COLUMN  80  CAUSES  DATA  PRINTOUT.  1 

/*EOJ 


/ 
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SUBROUTINE  SETF3(ORO,KST,NG,NTP, I , J,KND,NRH,M,*t* ) 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

C  SETF3  DETERMINES  WHICH  TRIO  OF  RELATIVE  HUMIDITIES* 

C  DATA  IS  TO  BE  GRAPHEOt  AND  FILLS  THE  THREE  CURVES' 

C  ORDINATE  AND  ABSCISSA  ARRAYS. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  ALLOCATIONS. 

C 

COMMON  S(2) .ZSI2).SDI2),NPR(2),X(3,25),Y(3,25) ,ZMAX(2) ,ZMIN 
$<2),TITLE(2,9f5),NPTS<  3) 

COMMON  /SETI3/Ft25,20,2#5),NF(20,2»5) 

REAL  0RD(25,20.2.5) 

C 

C  DETERMINE  THE  SEQUENCE  NUMBERS  OF  THE  STARTING  AND 

C  FINISHING  RELATIVE  HUMIDITIESt  ACCORDING  TO  HOW 

C  MANY  HAVE  BEEN  GRAPHED  BEFORE. 

C 

IF(KST-l) It 2»2 

1  KST*1 

KND=MI NO ( 3»  NRH ) 

GO  TO  3 

2  KST=KST+3 
KND*MIN0(KND+3»NRH) 

3  NS*0 

NCVS=KND-KST*l 

C 

C  FILL  THE  ORDINATE  AND  ABSC I SSA- VALUE  ARRAYS  FOR  EACH 

C  RELATIVE  HUMIDITY'S  CURVE  WITH  DATA  VALUES. 

C 

DO  A  K*KST • KND 
IF(NF(K,J,I J.EQ.O)  GO  TO  6 
NS-NS+I 

NPNTS=NF(K, JtI ) 

NPTSI NS )*NPNTS 
DO  4  L  =  1 » NPNTS 
X(NS.L)*F(L,K, J,IJ 

4  Y(NS»L)*ORD(LiK,J, I) 


145 


u  o  o  o  o 


o  o  o  noon 


THIS  MESSAGE  IDENTIFIES  THE  ORDINALS  OF  THE  RELATIVE 
HUMIDITIES  TO  BE  GRAPHEOi  AND  THE  NUMBER  OF  CURVES. 

WRI TE ( 6 i 5 )  KST.KND.NCVS 

5  FORMAT  {  *  K  RANGES  FROM*  I3f*  TO  *  la.*."!**1  CURVES  WILL* 
St* BE  DRAWN.' ) 

GO  ON  TO  JUSTIFY  GRAPH  DRAWING. 

CALL  JUSTFY  (NG.NTP.NCVS,  6.  7,  2) 

C  *****  IN  D2,  USE  THIS  FORM—  ***** 

CALL  JUSTFY ( NG»  NTP* NCVS »  6,  7) 

6  M=0 
RETURN  l 

7  RETURN  2 
END 
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/ 


ooo  ooon  r*>  o  o  o  o  ooooooooooo 


SUBROUT INEJUSTFY(NG»NTP,NCVS»*»*»* ) 

C  *****  im  02 1  USE  THIS  FORM —  ***** 

SUBROUTINE  JUSTFYI NG,NTP,NCVS, *,* ) 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

THIS  ROUTINE  DETERMINES  THAT  AT  LEAST  ONE  OF  THC 
POINTS  SENT  TO  IT  HAS  A  NON-ORIGIN  LOCATION,  INFER¬ 
RING  THAT  THERE  IS  INDEED  DATA  TO  BE  PLOTTED,  AND 
THUS  THAT  DRAWING  TICKED  AND  LABELLED  AXES  IS  JUS¬ 
TIFIED. 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


ALLOCATIONS. 

INTEGER  TITLE 

COMMON  S ( 2 ) , ZS ( 2), SD( 2 ) , NPR ( 2), XI 3,25) ,Y(3,25) , ZMA X ( 2 )  , ZMI  N 
$(2),TITLE(2»9,5)»NPTSI3) 

FOR  THE  X-AXIS,  IN  I T I At  IZE  THE  RANGE-END  VARIABLES  SO 
THAT  THEY  ARE  CERTAIN  TO  BE  CHANGED  BY  A  DATA  OR 
ASSIGNED  VALUE. 


ZX=-9.9E70 
ZN=+9. 9E70 

CHECK  ALL  THREE  CURVES  FOR  THE  HIGHEST  AND  LOWEST  X- 
VALUES. 

DO  1  NS=1, NCVS 
NPNTS=NPTS(NS) 

DO  1  N= 1 » NPNTS 
ZX=AMAXL(ZX,X(NS,N ) ) 

1  Z  N=  AM  I N 1  (  ZN » X  ( N  S » N  )  ) 

ZMAXI 1 )=ZX 
ZMIN(1)=ZN 

REPEAT  THESE  TWU  PROCEDURES  FOR  THE  Y-AXIS. 

ZX=-9. 9E70 
ZN-+9.9E70 
DO  2  NS=1,NCVS 
NPNTS=NPTS( NS) 

DO  2  N- I, NPNTS 
ZX=AMAX1(ZX,Y(NS,N  )  ) 

2  Z  N= AM  I N1 ( ZN , Y ( NS , N  )  ) 

ZMAX ( 2  )  =  ZX 

ZMI N( 2 ) -ZN 
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l_)  s_>  O  O  O  O  •_) 


o  n  o  o 


C 

C  A  VALUE  OF  FALSE  FOR  THIS  TEST  IS  THE  JUSTIFICATION 

C  FOR  PROCEEDING,  AND  PRINTING  OUT  THE  RANGE  OF  DATA 

C  VALUES. 

C 

IF  ( ZMAX  ( 1) .EQ.O..OR.ZMAX12J.EQ.O. )  GO  TO  4 

IF  JUSTIFICATION  FAILS,  THE  "BAD"  RETURN  ROUTE 
IS  ENTERED  AFTER  THE  FAILURE  MESSAGE  PRINTOUT. 


WRITE  (6,3)  ZMAX.ZMIN 

3  FORMAT ( •  THE  LARGEST  DATA  VALUES  PER  AXIS  ARE'  2G12.3, 

WHILE  THE  SMALLEST  ARE»2G12.3) 

CALL  SCALE(NG,NTP,NCVS,  8,  7) 

C  ###**.  in  D2,  USE  THIS  FORM— 

CALL  SCALE(NG, NTP.NCVS,  7) 

4  PRINT, 5  NG 

5  FORMAT (•  THERE  IS  NO  DATA  FOR  THE  GRAPH  WHICH  WOULD  HAVE' 

BEEN  THE  *I2,«TH.  CONTROL  IS  RETURNED  TO  MAIN  .') 

6  RETURN  1 

7  RETURN  2 

8  RETURN  3 
END 
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noooo  noon  o  o  o  o c o o o o o o o 


C 


SUBROUTINE  SCALE(NG,NTP,NCVS, * ,•  ) 

*#•**  IN  02*  USE  THIS  FORM —  •**•* 

SUBROUTINE  SCAL E ( NG, NT P , NCVS , *  I 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SCALE  SELECTS  INTERVAL  BOUNDS  WHICH  MINIMIZE  BOTH 
DEAD  SPACE  ON  THE  AXES,  AND  DISTRACTION  BY  OVERLY 
SIGNIFICANT,  MENTALLY  UNCOMFORTABLE,  AXIS  MARKINGS. 


ALLOCATIONS. 

INTEGER  TITLE 

COMMON  S(2),ZS(2),S0(2),NPR(2),X(3,25),Y(3,25),ZMAX(2) ,ZMIN 
S ( 2 ) ,TITLE(2,9,5),NPTS( 3) 

DIMENSION  R( 2) 

DATA  R/8.,5./ 

FOR  EACH  AXIS  IN  TURN  ... 

DO  20  NAX=1 *  2 
ZX=  ZMAX ( NAX ) 

ZN=  ZMI N ( NAX ) 

DETERMINE  THE  ORDER  OF  MAGNITUDE  OF  THE  BOUNDING  DATA 
VALUES,  AND  THEIR  RATIO. 

A=10»*( l+IFIXI ALOGIO(ZX) ) ) 

IF(ZN)  2,2,1 

1  B*IO»*( 1+IFIXIAL0GI01ZN) ) ) 

RAT IO*ZX/ZN 

GO  TO  3 

2  B*0. 

RAT IO*ZX 

DETERMINE  WHAT  FRACTION  OF  ITS  ORDER  IS  USED  BY  THE 
MAXIMUM  AXIS  VALUE,  AND  SET  THE  UPPER  BOUND  ACCOR¬ 
DINGLY. 


3  DO  6  N=2 , LO , 2 
IF(N-IO)  5, A, 4 

4  ZX* A 

IFIZX.LT. 0.78»A)  ZX=0.80*A 
IFIZX.LT. 0.6 2* A)  ZX*0.60*A 
GO  TU  8 
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o  o  o  n  o 


5  IF(ZX.LT.A/FLOAT(lO-( ) )  GO  TO  7 

6  CONTINUE 

7  I F ( N. EQ.8 )  ZX* A/2. 

IFIN.EQ.6)  ZX»A/4. 

IF(N.EQ.4)  ZX= A/5. 

IFIN.EQ.2)  ZX=A*0.15 

C 

C  ASSIGN  AN  APPROPRIATE  VALUE  TO  THE  LOWER  BOUND. 

C 

8  IFIRATIO. GT. 1000. .AND.ZN.LT. 100. )  ZN=0. 

I F ( RAT IO.GT • 1000. )  GO  TO  13 
IFIA-10.)  9,9,10 

9  ZN*0. 

GO  TO  13 

10  00  11  N32,8 , 2 

I F ( ZN. GT .6/ FLOAT ( N ) )  GO  TO  12 

11  CONTINUE 

12  IF (N.LE.4)  ZN*B/FLOATI N) 

IFIN.EQ.6)  ZN=0 . 15*B 

I F ( N. GE.8 )  ZN=B/10. 

C 

C  CALCULATE  THE  NUMBER  OF  VARIABLES  (••TIC  DIVISION1*), 

t  ANO  THE  NUMBER  OF  TICS,  BETWEEN  LABELS. 

C 

13  SD ( NAX ) =A/50. 

NPR I NAX ) *5 

IF (RAT IO.GT. 8. )SD(  NAX )*A/20. 

IF (RAT IO.GT. 8. ) NPR (NAX) =4 
C 

C  ASSURE  THAT  THE  AXIS  BOUNDS  BECOME  AN  INTEGER  HUL- 

C  TIPLE  OF  THE  TIC  DIVISION. 

C 

14  NTHI*i+IFIX(ZX/SD(NAX)) 

NTLO=IFIX(ZN/SD(NAX) ) 

I F ( (NTHI-NTLO) .LE. 10 )NPR(NAX 1*2 
C 

C  IF  LESS  THAN  SIX  TICS  WOULD  RESULT,  HALVE  THE  TIC  DI- 

C  VISION. 

C 

IFINTHI-NTLO-6) 15,  15,16 

15  SD  ( NAX  )  =SD( NAX ) /2. 

GO  TO  14 

16  ZMAX ( NAX )3NTHI »SD(  NAX )  • 

ZMIN(NAX)=NTL0*SD(NAX) 

C 

C  ASSURE  THAT  THERE  IS  A  NONZERO  AXIS  RANGE.  IF  THERE 

C  SHOULD  BE,  RETURN  01  TO  SETF3,  TO  TRY  THE  NEXT  TRIO 

C  OF  RH*  S.  IN  02,  RETURN  TO  MAIN  AND  CHANGE  GRAPH 

C  TYPES. 

C 
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o  o  o  o 


RANGE=ZMAX(NAX )-ZM I f J t  MAX ) 

IF  !RA<IGE>  17,17,19 

17  WRITE  (6,1H)ZMAX(NAX),ZMIN(NAX  ),  (  T  I  TLc  ( NAX  ,N  TP  ,  N »  •  N*1 , 5 
IB  FORMAT  I  '  ZMAX  =  aG12.3. •«  WHILE  ZMIN  =  *G12.3,*  SO  THAT* 

S,*  THERE  IS  APPARENT!  Y  A  NEGATIVE  OR  ZERO  RANGE  FOR  *, 

1* PLOTTING  */T20,5A4) 

RETURN  1 

CALCULATE  THE  SCALE  FACTOR,  AND  THE  NUMBER  OF  INCHES 
AC7UA'  TIC  SPACING. 

19  ZS ( NAX )=Rl NAX )/ RANGE 
StNAX)*ZS(NAX)*SDC NAX) 

20  CONTINUE 
C 

C  PRINT  OUT  THE  SUCCESS  MESSAGE,  SHOWING  THE  BOUNDS  AND 

C  THE  SCALING  VARIABLES*  VALUES. 

C 

WRITE  16,21 )NG, ( T I TLEl 2 , NTP ,N ) , N* 1, 5) , ( T ITLE ( 1, NTP, N > ,N*1 ,5 ) 
S, ( (S(NAX),SD(NAX), ZMAX ( NAX ),ZMININAX) ) ,NAX»1,2) 

21  FORMAT! *  FOR  GRAPH  *13,*,  OF  *5A4,*  VS  »5A4/T10, 

i*  ON  THE  X-AXIS  S  *  *G12.3,*,  SD  *  *G12.3,*,  ZMAX  *  » 
1G12.3,*,  ZMIN  =  *  G12 .3/T7,  *  AND  ON  THE  Y-AXIS  *6X,G12.3, 
27X,G12.3,2I9X,G12. 3) ) 

C 

C  GO  ON  TO  DRAW  THE  TICKED  AND  LABELLED  AXES. 

C 

CALL  PREPING, NTP, NCVS) 

RETURN  2 
END 
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o  o  o  n  o  o  o  o  o  non  noon  oooooooo 


/ 


SUBROUT  INE  PREPING, NTP.NCVS) 

ccccccccccccccccccrccccccccccccccccccccccccccccccccccccccccccccct 

PREP  CRAMS  THE  TICKED  AND  LABELLED  AXES  IN  PREPA¬ 
RATION  FUR  CURVE  CRAW  INC* 


ALLOCATIONS 

integer  title, tl 

CUPPCN  SI2I,ZSI2»,SCI2),NPR(2).XI3,25).VI3,25I,2MAXI2) ,ZHIN 
S(2),TITLE(2,9,5),NPTS(3) 

01  PENS  I CN  TL 129 ) , R ( 2 ) 

DATA  R/3.,5./ 

FOR  EACH  AXIS,  ASSIGN  SCALAR  VALUES  TO  ARRAY  VARI¬ 
ABLES,  FOR  FASTER  REFERENCE. 

CO  33  NAXa  1,2 
SCE«SCINAX) 

NP»NPR ( NAX ) 

SE*S ( NAX ) 

P*R ( NAX  I 

ASSIGN  VALUES  TO  TITLE  AND  AXIS-DRAWING  VARIABLES. 

IF! NAX- 1 1  1,1,2 

1  A«0  • 

XST*2 • 36 

YST«-0.656 

AL*0. 

AXIS2  CRAWS  THE  AXIS,  HERE  THE  ABSCISSA  ... 

CALL  AXIS2(0.,3.,P,A,SEJ 
GO  TO  3 

2  A*270. 

XST=-1.C0 
YS  T*0 • 86 
AL=  *90 . 

...  AND  HERE  THE  ORDINATE • 

CALL  AX1S2(0.«5.,P,A,SE) 

WRITE  CUT  THE  AXIS*  TITLE. 

3  CU  A  N«l,5 

4  TL(N)*TITLE(NAX,NTP,NI 

CALL  CPRINTI TL, 20.XST.YST, .16A.AL  ) 
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no  on 


n  o  o  o  onorio 


/ 


C 

C  THE  NEXT  FIVE  COMMENTS  LEAL  WITH  DETERMINING  THE 

C  TIC  LA3ELLINC  PARAMETERS • 

C 

NHALVE=0 

CE TERM  I  ME  THE  ORIGIN  TIC  VALUE.  ANO  WHETHER  IT  IS 
AN  INTEGER  MULTIPLE  OF  BOTH  THE  TIC  DIVISION  AND 
THE  LABELLING  CRCINAL. 

5  TVC=ZMIN(NAX) 

NT ICS= I F I Xt  P/SE  I  ♦  l 
NSTCNG=1 

6  C*TVC/SCE/FLOAT(NP) 

NC*IFlX(C».COOl  ) 

IF(C-FLCAT(NC)-1.E-A)8,8,7 

IF  NCT,  SET  UP  TO  DETERMINE  WHETHER  THE  SECOND  OR 
SUDSECUENT  T  IC  MIGHT  DE . 

7  NSTCNG=NSTCNP+1 
IFINSTCNG.EC.20)  GC  TO  II 
TVC=TVOfSCE 
GO  TO  6 

C 

C  IF  THE  FIRST  TIC  WILL  BE  MORE  THAN  ONE  THIRD  THE  WAY 

C  COWN  THC  AXIS,  HALVE  THE  TIC  DIVISION,  AND  ADJUST 

C  THE  CTHER  SCALING  PARAMETERS. 

8  IF(FLCAT(NSTCNG).LE.0.33A»FLOA TIN TICS). AND. NTICS.GT.5) GO  TO  10 
SCE*SLE/2. 

SE=SE/2. 

NTH  I  *1  ♦  IF  IX  I  ZM.  AX  (NAXJ/SDE) 

NTLC«IFIX(2P!N(NAX  )/SCE) 

ZM AX  I NAX ) »SC E*FLOAT ( NTH D 
ZMIN(NAX)=SCE*FLUATINT10) 

NTICS*NTHI-NTLO 

RANGE*  ZM  AX  (N  AX  )-ZM  IMfJAX) 

ZS(NAX)=R(NAX)/RANCf 
IF(NTICS.GT.10.ANC.NP.GT.2)NP«NPZ2 
IFINP.GT.2)  NP*NP/ 2 
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m  r  non  noon 


...  WRITE  A  MESSAGE  CkA  h I NO  AY  TENT  I  UN  TO  THl  CHANGE, 

AWO  REOE  TERM  IflL  THl  TIC  LAB  ELL  I  MG  PARAMETERS. 

WRITEI6,  9)  NAX,NTP,Se,$OE,MP 

9  FORMAT ( 30X, • S, SC,  NPR  $03  (  *  I  2 ,  ' , • l 2 , ' )  HAVE  BEEN  HALVED, '/»  T 

IHEIR  VALUES  ARE  MOW  *  1PE10.2i  *  t  '  1PE10.  2.  * .  M2,'  RESPECTIVELY.'} 

NHALVt  =  NH AL VE*  1 
IF(KHALVE-S»  5,5,10 

10  IFISLE.LT.  1  .  )  (MARK*  I  F  I X  (  1 0.  •  •  (I  F  I  X  ( ALUG 1  0  (  SDE  )))  /  SDE  » 

IFlSCE.GE.k . )  NMARK  *MP 
IFINKAHK.CT  ,5J  N'MARK * S 
IF(NMARK)  33,33,11 

U  WRITE  (6,12)  MAX,NSTCMG,WTICS,NHARK 

12  FORMAT  (  T 10,  'FOR  MAX  *M1,».  NSTCNG,  NT  ICS,  AND  NMARK  ARE  •  I  2  ,  '  ,  •  I  2  , 
1*.  AND  M2,',  RESPECT1VELYM 
I F ( NS  f  CMC .EC .20 )  CO  TO  33 

AT  EACH  APPROPRIATE  TIC  ... 

CO  32  N« NSTCNG, NT  ICS, MM ARK 
M*N-1 


DETERMINE  THE  TIC  VALUE. 

TV»SCfc  «  FLOAT  (Ml  MM  INI  MAX  ) 

1  T  V*  I  F  I  X  (  TV*O.A) 

I F (  I T  V )  13,13,16 
13  I  F  (  TV)  l*« ,  1  A  ,  1  5 
i<*  ITVsfl 

nc  *  l 

C  CONVERT  THE  NUMERIC  INFORMATION  (I  TV)  TO  ALPHABETIC 

L  (LAJ),  AS  R ECU  I  RED  BY  THE  TIC-LABEL  DRAWING  RUU  TINE. 

C 
C 

C  FOR  SMALL  LABELS  ... 

t 

CALL  ZUCC  1*4 C LASJ ,  I  T V  ) 

GO  TO  2 1 

15  NCM«  mxUHSI  ALOGIO(TV)  )  ) 

ITV*1F I X ( Tv • 10  *  «  * ( NC  )  ) 

CALI  BCCINT (LAB, l TV  ) 

GO  TO  21 

(6  I f  t I  TV- 10000)  17, 1H,  13 
U  CALI  ZCCCINILAB, I  TV) 

GO  TO  19 

19  ITVFI* 1TV/10C30 
1  TVLO  I  TV 


15A 


<**  o  o  noo  o  o  o  n  ooo 


/ 


•••  CR  LARGE* 

CALL  ZBCCIN(LABFI, ITVHI J 
CALL  BCDINT ( LA3LD,  ITVLO) 

19  NC=IFIX( ALOCiO(TV) )  *1 

20  IF  ( I T V-1C000 )  21,27,27 

21  IF(NAX-l)  22,22,26 

22  IFINC-3)  23*23,24 

CETERMINE  THE  LCCATION  OF  THE  LABEL  WITH  RESPECT  TO 
THE  TIC,  DEPENDING  ON  TIC  VALUE  AND  AXIS. 

23  CTR=.362-.05»FL0AT  (NO 
GO  TO  25 

24  CTR=0. 

25  Xl*SE»FLOAT (H)-CTR 
Yl*-0.20 

GO  TO  31 

26  X 1  =  - . 4 8 
Yi=SE«FLOAT (MI-.05 
GC  TO  31 

27  IF  (NAX-l)  28,28,29 

28  XHI^SE»FL0AT(M)-.05«FL0AT(NC) 

XLOXhl*.  l»4. 

Yl*-.20 
CO  TO  30 

29  XH I 88 
XLO  *  -.4C 
Yl*SE«FLCAT (HI-.05 

CRAW  THE  TIC  LABEL • 

30  CALL  CPRINT(LABFI,4,XHI,Y1,0.1C,0. ) 

CALL  CPRINT(LABL0,4,XL0,Y1,0.1C,0. ) 

GC  TO  32 

31  CALL  CPRINTILAB,  4,X1, Yl*.  10,0.  ) 

32  CONTINUE 

33  CONTINUE 

WRITE  (6,34)  I  TV 

34  FORPATC  THE  LAST  TIC  LABEL  ON  THE  Y-AXIS  !  S  *  1 8 ) 

PROCEED  TO  CRAW  THE  POINTS  UN  THE  PREPARED  AXES* 

CALL  PCINT(NCVS) 

RETURN 

ENC 
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SUBROUTINE  POINT(NCVS) 

ccccccccccccccccc  cccccccccccccccccccccccccccccccccccccccccccccccc 
c 

c  POINT  PLACES  THE  PRUPER  SYMBOL  AT  THE  PROPER  PLACE 

C  WITH  RESPECT  TO  THE  PREPARED  AXES,  AND  COUNTS  THE 

C  NUMBER  OF  POINTS  DRAWN  ON  THE  GRAPH.  IT  ALSO  DRAWS 

C  A  MESSAGE  ON  THE  GRAPH  SHEET,  SHOWING  THAT  THE  TRIA- 

C  GLES  MARK  THE  FIRST  CURVE,  SQUARES  THL  SECOND,  AND 

C  DIAMONDS  THE  THIRD. 

C 


C 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

ALLOCATIONS. 


COMMON  S (2) , ZS  t 2), SD<  2),NPRt  2),X( 3, 25), Y( 3,25) ,ZMAX(2)  ,ZMIN 
i(2),TITLE(2,9,5),NPTS( 3) 

NPT  =.0 

DO  A  NS= 1 , NCVS 

REVERSE  THE  ORDER  OF  DRAWING  POINTS  SO  AS  TO 
ELIMINATE  WASTE  PEN  MOTION.  FOR  THE  D2  PROGRAM, 

THE  RELATIONAL  OPERATOR  .NE.  APPEARING  BELOW  IS 
REPLACED  BY  .EQ.. 

IFCNS..NE.2I  CALL  SFFL(NS) 

IFINS-2)  1,1,2 

1  NSY^iNS 
HT=0 . 1 
GO  TO  3 

2  NSY=6 
HT*0. 2 

3  NPNTS*NPTS(NS) 

DO  A  N= 1 , NPNTS 

I  F  (  X  (  NS  ,-N)  .EQ.O..AND.Y(NS,N).EQ.O.  )  GO  TO  4 

CALCULATE  THE  X  AND  Y  COORDINATES,-  IN  INCHES. 

A>(X(NSf N)-ZMIfl(l) )  »ZS(  1  ) 

B*(Y(NS,N)-ZMIN(2) )*ZS( 2) 

TPRINT  DRAWS  THE  PROPER  SYMBOL,  CENTERED  ON  (A,B), 

CF  HEIGHT  HT. 

CALL  TPRINT (NSY, A, B,HTJ 
NPT*NPT ♦ 1 

4  CONTINUE 

WR I T  E ( 6 , 6 )  NPT 

6  FORMAT { T2, I  3, •  POINTS  WERE  DRAWN.*) 
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THE  REMINDER  OF  THE  CURVE  DRAWING  ORDER  IS  DRAWN 
ON  THE  GRAPH  PAPER  AFTER  COMPLETING  THE  THIRD  GRAPH. 

I F ( NG-3 )  8,7,8 

7  CALL  TPRINT  ( 1 ,0 . ,+8 .0, 0 . 1 ) 

CALL  SPRINTCIS  CRAWN  FIRST* ,14, 0.2, +.3. 0,0. 1,0.  ) 

CALL  TPRINT  ( 2 , 1 . 70 , *8 .0, 0 . 1) 

CALL  SPRINT ( *  SECOND  AND',  10, 1.90, +  8. 0,0. 1,0.) 

CALL  TPRINTC6, 3. 1, ♦3.0,0.21 

CALL  SPRINT  ( • TH IRD* , 5, 3. 3, +8.0, 0. 1 , 0 . ) 

8  RETURN 
ENU 


157 


o  n  o  o  nooo  n  o  o  o  onnn  oooooooooooooooooo 
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SUBROUTINE  SHPLlNS) 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SHFL  CHANGES  THE  ORDER  IN  WHICH  POINTS  ARE  DRAWN, 

TC  ELIMINATE  NEEDLESS  PEN  TRAVEL. 

THE  RESULTING  DRAWING  DIRECTIONS  ARE— 


D1 

D2 

CURVE 

1 

R-L 

L-R 

CURVE 

2 

L-R 

R-L 

CURVE 

3 

R-L 

L-R 

WITHOUT  THE  OPERATION  OF  SHFL,  ALL  CURVES  WOULD  BE 
DRAWN  FROM  RIGHT  TO  LEFT. 


ALLOCATIONS. 

COMMON  SI2) ,2Si2),SDI2),NPRI 2Jt XI 3, 25 ) , Y I  3, 25 ) ,2MAXI2> ,ZMlN 
$12) ,TITLEI2,9,5) ,NPTSI3) 

NP=  NPTSINS) 

ESTABLISH  THE  POINT  NUMBERS  TO  BE  INTERCHANGED,  THE 
M-TH  AND  THE  NHI-TH. 

DO  I  M= I , NP 
NHI*NP-M+1 

TEMPORARILY  RENAME  THE  VALUE  OF  THE  LOW- SEQUENCE- 
NUMBER  COORDINATES. 

TI  =  X I  NS  »  M) 

T2=YINS,M) 

PUT  THE  HIGH  SEQUENCE  NUMBER  VALUE  WHERE  THE  LOW  FOR¬ 
MERLY  WAS. 


XINS, M)=XINS,NHI ) 

Y(NS,M)=YINS,NHI  ) 

GIVE  THE  HIGH  END  LOCATION  THE  RENAMED  VALUE  OF  THE 
LOW  ENO  VARIABLE. 

XINS.NMI I *T l 
YINS.NHt )=T2 

IFUNP-MJ  .Gc.NP/2)  RETURN 
l  CONTINUE 
RETURN 
END 
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//D2  JOB  (XXXX,XX,035,06,1000, ), »S$$BUDRAYGER$*S»,CLASS=M 
/•PASSWORD 

//  EXEC  F4GCXM  FORT  G  COMPILE  (NODECK),  EXECUTE,  CLASS  M 

//FORT  *  SYS  I N  DO  • 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


c  c 

C  02  DRAWS  REPORT-READY  GRAPHS  OF  TYPES  4 — 9,  OF  C 
C  THE  PERMITTIVITIES  AND  HUMID-CAPACITANCE  CHANGE  C 
C  VS.  THE  FRACTIONAL  WEIGHT  INCREASE  BY  ADSORPTION  C 
C  OF  AMBIENT  WATER  VAPOR.  AXES  ARE  SCALED  FUR  EACH  C 
C  GRAPH.  C 
C  THE  MAIN  PROGRAM  READS  THE  DATA,  AND  CONTROLS  THE  C 
C  FLOW  THROUGH  THE  TEMPERATURE-MODE  COMBINATIONS,  C 
C  DRAWING  THE  GRAPH  OF  EACH  TYPE  IN  SUCCESION.  EVEN  C 
C  NUMBERED  GRAPH  TYPES  HAVE  CURVES  FOR  THREE  SPE-  C 
C  CIFIC  HIGH  FREQUENCIES,  ANO  ODD  TYPES  FOR  THREE  C 
C  LOW  FREQUENCIES.  C 
C  C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  ALLOCATIONS. 

C 

INTEGER  TITLE 

REAL  R(2),EP(25,20,2,5),EDP( 25,20,2,5) , DC ( 25 ,20, 2 , 5 ) ,0(2,201) , 
10RDNT (25, 20, 2, 5) 

COMMON  S(2),ZS(2)»SD(2)»NPRl2)»X(3»25)»Y(3,25)»ZMAX(2) ,ZMlN 
S(2) ,TITLE(2,9,5),N?TS(3i 
COMMON/SET49/WR(20,2,5),tSP(6f 20, 2 , 5 ) , KM< 2 , 5 ) 

NAMEL I  ST/TWO/NT, KM, LSP.WR, EP, E OP ,0, TITLE 
C 

C  READ  THE  DATA  FROM  THE  DATA  SET. 

C 

READ  (2)  NT,KM,0,TITLE,LSP»WR,EP»EOP,DC 
C 

C  IF  DESIRED,  WRITE  ITS  CONTENTS. 

C 

READ  (5,1)  QPRINT 
1  FORMAT  ( 79X , Fl.0  9 

I F ( QPRINT .NE.O . )  WRITE ( 6, T WO  ) 

C 

C  INITIALIZE  THE  CALCOMP  PACKAGE,  THE  GRAPH  SHEETS,  AND 

C  THE  CONTROL  INDICES. 

C 

CALL  PLSTRT ( 2222,22) 

CALL  PLOT (3*, 1.50, -3) 

NG=0 

NSP=0 

NGADVQ*! 


159 


ooooo  ooooo  ooo  o  <-i  o  noon 


/ 


C 

c  T-MC' S  ARE  I TERATIVELV  ESTABLISHED  HERE. 

DO  13  1*1, NT 
DU  13  J*l,2 

FOR  EACH  TEMPERATURE-MODE  COMBINATION*  SET  AND/OR  STEP 
THE  GRAPH  TYPE  INDEX. 

NT  P=3 

2  NT P*NTP*l 

SET  A  NEW  ORIGIN  WHEN  NECESSARY. 

IF(NGADVQ)  3*4*3 

3  NG=NG*l 

CALL  PL0T(O(l,NG),O(2,NG),-3) 

PRINT  A  MESSAGE  IDENTIFYING  EACH  GRAPH. 

4  WRITE  (6,5)  NTP*NG* I* J 

5  FORMAT ( 'O' » '  NTP  *'I2,'»  NG  *  '13, '»  I  =  '12,*,  J  *  '12) 

ELIMINATE  ALL  UNESSENTIAL  AND/OR  UNCONTROLLED  NUMBERS 
FROM  THE  BLUCK  OF  STORAGE  WHICH  WILL  BE  FILLED  WITH 
GRAPH  DRAWING  INFORMATION. 

DO  6  NS*l»3 
NPTS(NS)»0 
DO  6  N*I ,25 
X( NS»N)*0. 

6  Y( NS ,N ) *0. 

NGADVQ*! 

PROCEED  WITH  SETTING  UP  THE  GRAPH  OF  THE  PRESENT  TYPE, 
WITH  THE  PROPER  ORDINATE  ARRAY  (FIRST  ARGUMENT),  AND 
FRECUENCY-RANGE-STARTING-INDEX  (SECOND  ARGUMENT). 

GO  TO  (7, 7, 7, 7, 8, 9, 10, Li, 12), NTP 

7  CALL  SETL6(  EP, I,NG,NTP, I , J , NGADVQ »  13,  2) 

8  CALL  SETL6(  EP,4»NG»NTP, I, J, NGADVQ,  13,  2) 

9  CALL  SETL6(EDP, i»NG»NTP» I, J, NGADVQ,  13,  2) 

10  CALL  SETL6( EDP, 4, NG,NTP, I, J, NGADVQ,  13,  2) 

11  CALL  SETL6(  DC, l, NG, NTP, I, J, NGADVQ,  13,  2) 

12  CALL  SGTL6(  DC, 4, NG, NTP, I, J , NGADVQ,  13,  13) 

13  CONTINUE 
CALL  PLTEND 
STOP  222 


160 


/ 


c  •*■■■*«*•*******•*•******•**•*****•*••*•**•*•* 

C  •  * 

c  *****the  jgb  control  language  CARDS  DELOW*» •** 

c  *« ***8ETWE EN  "/•"  AND  "/*EOJ"  . 

C  «****ARE  PLACED  DEI-  IND  THE  DECK  FOR  SUB-  **••• 

C  •■•■•ROUTINE  SHFL  WHEN  ACTUALLY  RUNNING  *•**• 

C  . THE  PROGRAM.  . 

C  *  • 

c  . . . 

ENC 

/* 

//GO.  FT02F0Q I  CC  CSN  =  TFREE.CAY.FOR.D2,UNI T  =  2321 , 

//  VCL=SER=DCELLI,DISP=OLD 

/ /G0.FTI8F00I  CC  UNI T=TAPE7, DCB* I RECFM=U ,BLKSI 2E =1 440 ,0EN=1 , TR TCH=ET ) , 
//  DISP=OLD,LABEL=( ,BLP) . VOL=SER=OUTP T , D$N=CALCOMP 

/ /GO. SYS  I N  CC  • 

A  BLANK  IN  COLUMN  00  CELtTES  CATA  PRINTOUT. 

THE  1  IN  COLUMN  80  CAUSES  CATA  PRINTOUT.  I 

/•EOJ 
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SUBROUTINE  SETL6(0RCNT,NSPST,NG,NTP,I , J,M,# ,«) 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SETL6  ENGAGES  TEE  PROPER  TRIO  OF  FRCQUENC Y-PARA- 
PETERS,  AND  FILLS  THE  THREE  CURVES*  ORDINATE  AND 
ABSCISSA  ARRAYS. 


ALLOCATIONS. 

INTEGER  TITLE 

REAL  GRCNT ( 25f  20*  2 • 5 ) 

COPPCN  S(2),ZS(2),SC(2),NPR(2),X(3,25),Y(3,25),ZMAX(2),ZPIN 
S(2) .T1TLE(2,9,5),NPTS( 3) 

C0PPCN/SET49/WR(20,2, 5 ) , LSP ( 6, 20, 2 , 5 ) , KMI 2 , 5 ) 

INITIALIZE. 


Z  =  0 

KBAC=0 
NCVS  *3 

NSPENC=NSPST*2 

IF  THERE  IS  NO  CATA  AT  THIS  TEMPERATURE-MODE  COMBINA¬ 
TION,  RETURN  TO  MAIN  AND  ITERATE  TO  THE  NEXT. 

I F ( KP ( J  »  I  )  )  7,7,1 

FOR  EACH  OF  THE  SPECIAL  FRECliENC  I ES  IN  TURN,  ... 

1  CO  6  NSP=NSPST,N SPEND 
NS=NSP 

IF(NSP.Gc.4)  NS=NSP-3 
NRH=KP I J , I ) 

NPTSINS )=NRF 

...  AND  FOR  ALL  THE  RELATIVE  HUM  1 0 1  TIE  S — 


CO  5  K= 1 , NRF 
IF(NCVS)  2,2,4 
2  NSP1=NSP* 1 
NSP2=NSP*2 
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TEST  TO  ASSURE  THAT  CATA  HAS  BEEN  TAKEN  AT  AT  LEAST  ONE 
CF  THE  FREQUENCIES  IN  THE  PRESENT  GROUP.  IF  NOT,  WRITE 
A  MESSAGE ,  ONLY  CNCE  PER  T-MC,  THAT  THERE  IS  NO  DATA. 

IF(LSP(NSP,K,J, I ) • EC.O • ANO.LSP ( NSP I , K, J  » I ) .EQ. O.AN^.LSP ( 
$NSP2,K, J,  I)  .EQ.0.ANC.NTP.LE.5)  WR I TE ( 6, 3) I , J ,K ,NSPST 

3  FORMAT  (  •  THERE  IS  A  ZERO-TRIPLET  OF  LSP"S  FOR  I,J,K,», 

STARTING  NSP  =  *3(I3**,*)**  M2) 

IF  NO  RELATIVE  FuMICITY  FOR  THE  T-MC  HAS  OATA  FOR  THE 
FRECUENCY  GROUP,  RETURN  SO  AS  TO  CHANGE  ONLY  THE  GRAPH 
TYPE. 

IF  (LSP(NSPI,K,  J,  D.EQ.O)  KBAD=KBAD*l 
IF(KBAD.EQ.NRH)  RETURN  2 

4  I F ( LSP ( NSP, K, J, I )  .EQ.O)  GO  TO  5 
X(NS,K)=RR(K,J, I) 

Y(NS,K)=.ORDNT(  LSP(NSP,K,  J,  I  ),K,  J,  I  ) 

Z  =  X(NS,K)+Y(N'S»K)  +  Z 

5  CONTINUE 

IF  THE  SUM  OF  ALL  COORDINATE  VALUES  ON  A  CURVE  IS 
NOT  ZERO,  INCREASE  THE  EXPECTED  NUMBER  OF  CURVES 
eY  CNE. 

IF(Z.GT.O)  NCVS=NCVS*:i 

6  CONTINUE 

IF  THERE  ARE  POSSIBLY  SOME  CURVES  TO  DRAW,  PROCEED  TO 
SUBROUTINE  JUSTFY. 

IF(NCVS.GT.O)  CALL  JUSTFY  ( NG, NTP, NCVS ,  7,  8) 

7  M*0 
RETURN  l 

8  RETURN  2 
ENC 
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IV .  Discussion 

The  research  of  this  eighth  semi-annual  period  has  yielded  technical 
findings  that  bear  on  problems  of  interest  to  the  Department  of  Defense. 

Moreover  from  this  research  has  come  the  foundations  upon  which  our  future 
investigations  will  build. 

The  data  reported  here  taken  on  decectors  utilizing  the  PME  and  Dember 
effects  in  gold-doped  silicon  suggest  that  such  devices  may  find  use  in  special 
applications  for  photodetection  in  the  visible  and  near  infra-red  regimes. 

Greater  promise  yet  is  shown  by  new  detector  structures  to  be  described  in 
the  next  Scientific  Report.  For  example,  gold-silicon  Schottky  barrier  photo¬ 
diode  employing  a  new  mask  geometry  in  the  form  of  a  grating  shows  advantages 
in  comparison  with  conventional  silicon  Schottky  barrier  photodiodes  and  most 
pn  junction  photodiodes  now  available.  The  relevant  properties  include  the 
following:  fast  response  time,  high  responsivity ,  flat  spectral  reponse  between 
0.4  and  1.0  microns,  high  quantum  yield  and  low  noise  performance.  These  studies, 
together  with  those  reported  here  describing  fundamental  electronic  parameters  in 
Cr-doped  gallium  arsenide  and  in  n-type  indium  antiminide,  all  bear  on  the 
design  of  new  types  of  detectors  that  potentially  may  be  of  interest  in  many 
military  applications. 

Transistors  form  the  heart  of  the  electronic  systems  used  in  numerous 
applications  of  interest  to  the  Department  of  Defense:  detection,  telemetry, 
navigation,  control,  surveillance,  medicine,  etc.  In  this  report  we  describe 
two  studies  critical  to  the  design  and  optimum  application  of  transistors. 

The  first  deals  with  an  effect  that  can  lower  the  breakdown  voltage  of  a 
transistor,  constituting  a  failure  mode  in  some  circuit  applications.  This 
effect,  termed  premature  punch  through,  is  described  quantitatively  with 
instructions  given  as  to  how  it  can  be  controlled  or  essentially  eliminated. 

The  second  study  deals  with  the  incorporation  of  systematic  model  selection 
into  the  computer  simulation  of  semiconductor  circuit  behavior.  Compared  with 
other  procedures  now  available,  the  method  described  here  simplifies  computation, 
saves  computational  cost,  and  enables  the  simulation  of  larger  circuits.  It 
provoked  sufficient  interest  among  designers  of  circuits  for  radiation  environ¬ 
ments  that  the  essentials  of  the  method  constituted  one  of  five  invited  papers 
at  the  1971  IEEE  Nuclear  and  Space  Radiation  Effects  Conference  (Durham,  N.  H., 

July  1971). 
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Properties  exhibited  by  amorphous  semiconductors  make  them  potentially 
useful  in  a  variety  of  applications  that  may  involve  exposure  to  irradiation. 
From  the  research  reported  here  has  come  some  indications  of  the  ultimate 
limits  of  performance  to  be  expected  from  amorphous  semiconductors.  More 
comprehensive  indications  will  appear  in  future  Scientific  Reports. 
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